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Foreword 


Aircraft  engine  companies  in  the  U.S.  have  been  using  reacting  flow  QFD 
computer  codes  for  many  years  as  a  design  tool  for  new  combustor  designs  and 
for  the  analysis  of  flow  inside  existing  combustors  to  solve  design  problems.  In 
spite  of  this  record  of  successful  use  of  CFD  codes  for  combustors,  the  engine 
companies  all  agree  that  there  is  a  need  for  combustor  codes  with  capabilities 
much  greater  than  those  currently  in  use.  Generally,  current  codes  are  limited  in 
accuracy,  are  time-consuming  in  establishing  a  grid  system  which  represents  the 
complex  geometries  of  a  combustor,  present  difficulties  in  obtaining  a  stable  so¬ 
lution  (limited  robustness),  and  require  large  execution  times  on  mainframe  su¬ 
percomputers  for  practical  problems.  At  the  same  time,  the  down-sizing  that  has 
occurred  within  these  companies  has  left  them  with  very  limited  resources  to  in¬ 
dependently  pursue  an  improved  combustor  code.  They  envision  a  partnership 
with  the  U.S.  Government  (NASA,  Army,  Air  Force,  etc.)  to  be  a  more  cost-ef¬ 
fective  approach  for  producing  improved  combustor  code  capability  that  can  be 
used  by  all  within  the  U.S.  without  compromising  proprietary,  company-owned 
combustor  designs.  We  at  Lewis  Research  Center  are  pursuing  this  approach  in 
the  following  way. 

Efficient  and  accurate  CFD  algorithms  have  been  developed  for  compress¬ 
ible  flow  applications  over  the  last  several  years.  These  techniques  have  been 
modified  for  computational  efficiency  at  the  low  Mach  numbers  which  are  typ¬ 
ically  found  in  combustors.  Added  to  this  numerical  algorithm  are  state-of-the- 
art  physical  models  for  liquid  fuel  spray,  complex  turbulence,  and  finite  rate 
chemistry  for  combustion.  This  Workshop  marks  the  formal  release  of  ALL- 
SPD-3D  to  the  U.S.  combustion  community.  Concurrent  with  the  ALLSPD  ef¬ 
fort,  we  are  pursuing  even  more  promising  CFD  approaches  such  as 
unstructured  grids,  finite  element  flow  solvers,  pd  adaptive  re-meshing.  The 
hallmark  of  these  longer  term  efforts  is  the  initiation  of  the  National  Combustor 
Code  (NCC),  a  workhig  partnership  with  several  U.S.  companies,  NASA,  and 
other  agencies  to  jointly  develop  the  next  generation  code.  The  first  version  of 
the  NCC  is  expected  to  be  available  in  1998. 

We  believe  that  the  ALLSPD-3D  code  being  released  at  this  time  will  be  a 
major  step  forward  in  computational  capability  for  its  users.  It  should  provide  a 
user-friendly  platform  for  executing  the  analyses  of  complex  reacting  flow  prob¬ 
lems  with  state-of-the-art  physical  models.  Greater  physical  model  fidelity  and 
an  improved  gridding  methodology  are  being  pursued  through  the  NCC  partner- 
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ship  and  will  be  available  in  the  future.  In  the  meantime,  it  is  our  desire  that 
ALLSPD-3D  provide  computational  capability  to  the  U.S.  industry  until  the 
fixiits  of  the  NCC  ripen  and  are  harvested. 

The  ALLSPD  combustor  code  has  been  in  development  at  Lewis  since 
1991 .  An  initial  2-D  version  of  the  code  was  released  in  1993.  Some  of  the  major 
features  of  the  present  3-D  code  are  a  stochastic  fuel  spray  model,  a  k-£  turbu¬ 
lence  model,  a  library  of  finite  rate  combustion  models,  internal  obstacles,  mul¬ 
tiple  grid  zone  capability,  and  parallel  processor  execution.  ALLSPD-3D  has  a 
graphical  user  interface  to  increase  the  comfort  level  of  the  user  and  has  been 
run  on  supercomputers,  a  variety  of  scientific  workstations,  and  a  cluster  of 
workstations.  The  code  has  successfully  computed  a  wide  range  of  reacting  and 
non-reacting  flow  problems  of  increasing  geometric  complexity.  The  full  vali¬ 
dation  of  ALLPSD-3D  will  be  an  interactive  process  with  the  users  as  they  gam 
experience  with  the  code.  The  Lewis  Research  Center  will  work  with  the  ALL- 
SPD-3D  users  as  they  perform  calculations.  We  are  looking  forward  to  a  pro¬ 
ductive  working  relationship  with  you,  the  "customers,"  of  our  code 
development  effort! 


Dr.  Edward  J.  Mularz 

Chief,  Engine  Components  Division 

U.S.  Army  Research  Laboratory  -  Vehicle  Propulsion  Directorate 
NASA  Lewis  Research  Center 
November  14, 1995 


(Dr.  Mularz,  as  Chief  of  the  Aerothermochemistry  Branch  at  NASA  Lewis 
Research  Center,  initiated  the  ALLSPD  Code  effort  and  guided  its  development 
until  the  recent  move  to  his  new  position.) 
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Abstract 


This  report  describes  the  ALLSPD-3D  computer  program  for  the  cal¬ 
culation  of  three-dimensional,  chemically  reacting  flows  with  sprays.  ALL- 
SPD-3D  employs  a  coupled,  strongly  implicit  solution  procedure  for 
turbulent  spray  combustion  flows.  ALLSPD-3D  utilizes  generalized  cur¬ 
vilinear  coordinates  with  both  multi-block  and  flexible  internal  blocks^e 
capabilities  to  handle  complex  geometries.  The  ALLSPD-3D  algorithm, 
which  employs  preconditioning  and  eigenvalue  rescaling  techniques,  is  ca¬ 
pable  of  providing  efficient  solution  of  flows  over  a  wide  range  of  subsonic 
Mach  numbers.  The  chemistry  treatment  in  the  code  is  general  enough  that 
an  arbitrary  number  of  reactions  and  species  can  be  defined  by  the  users. 
In  addition,  for  general  industrial  combustion  applications,  the  code  pro¬ 
vides  both  dilution  and  transpiration  cooling  capabilities.  A  stochastic 
spray  model  and  an  efficient  method  for  treating  the  spray  source  terms  in 
the  gas-phase  equations  are  used  to  analyze  evaporating  spray  droplets.  Al¬ 
though  written  for  three-dimensional  flows,  the  code  can  be  used  for  two- 
dimensional  and  axisymmetric  flow  computations  as  well.  The  code  is  writ¬ 
ten  in  such  a  way  that  it  can  be  run  on  various  computer  platforms  (super¬ 
computers,  workstations,  and  parallel  processors)  while  the  GUI 
(Graphical  User  Interface)  should  provide  a  user-friendly  tool  for  setting 
up  and  running  the  code. 
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/ ;  Introduction 


It  has  been  about  two  decades^  since  the  U.S.  Government  sponsored  a  ma¬ 
jor  Computational  Fluid  Dynamics  (CFD)  code  development  effort  focused  on 
the  analysis  of  reacting  flows  inside  combustors  of  gas  turbine  engines  or  liquid 
rocket  engines.  During  the  ensuing  years,  many  advances  in  CFD  have  been  re¬ 
alized,  especially  in  the  area  of  compressible  flows.  Applying  these  advanced 
techniques  in  a  new  reacting  flow  code  has  the  potential  for  significant  improve¬ 
ments  in  accuracy,  efficiency,  and  robustness  over  the  current  design  codes  be¬ 
ing  used  in  the  industry.  For  the  past  few  years,  NASA  Lewis  Research  Center 
has  undertaken  the  development  of  a  new  combustor  analysis  code,  ALLSPD- 
3D,  which  employs  state-of-the-art  computational  technology  and  physical 
models.  During  this  development,  NASA  Lewis  will  provide  from  time  to  time 
a  current  version  of  the  code  to  the  aerospace  industry,  and  will  solicit  feedback 
from  users  on  the  code's  shortcomings  and  suggestions  for  further  improve¬ 
ments.  (Please  direct  any  questions  or  comments  regarding  the  ALLSPD-3D 
code  to  allspd@lerc.nasa.gov.) 

The  current  ALLSPD-3D  code  has  evolved  from  the  two-dimensional  ver¬ 
sion  released  in  June,  1993.  Besides  extension  to  three  dimensions,  the  new 
code  features  several  additional  improvements,  including  a  user-friendly  graph¬ 
ical  user  interface  (GUI),  multi-platform  capability  (supercomputers,  worksta¬ 
tions,  and  parallel  processors),  improved  turbulence  and  spray  models, 
generalized  chemistry,  a  new  solution  procedure,  and  more  general  boundary 
conditions.  In  this  report,  the  fundamental  structure  of  the  ALLSPD-3D  code, 
numerics,  physical  models,  and  I/O  procedure  are  described  in  detail. 

The  numerics  of  the  ALLSPD-3D  code  are  based  on  modem  compressible 
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algorithms  to  take  advantage  of  advances  made  in  CFD  during  the  past  two  de¬ 
cades.  It  is  well  known  that  numerical  methods  developed  for  compressible 
flows  are  often  ineffective  at  low  Mach  numbers.  There  are  two  reasons^’^  for 
this  difficulty.  First,  the  system's  eigenvalues  become  stiff  at  low  flow  veloci¬ 
ties.  In  theory,  this  can  be  circumvented  by  using  implicit  numerical  schemes. 
In  practice,  however,  a  large  approximate  factorization  error  in  multiple  dimen¬ 
sions  limits  an  optimal  CFL  number  beyond  which  convergence  slows  down.  At 
low  Mach  numbers,  there  are  large  disparities  among  CFL  numbers  and  eigen¬ 
values  in  three  directions  leads  to  further  uncertainties.  Therefore,  the  optimal 
CFL  numbers  can  be  difficult  to  determine.  Second,  the  pressure  term  in  the  mo¬ 
mentum  equations  becomes  singular  as  the  Mach  number  approaches  zero, 
yielding  a  large  round-off  error.  This  behavior  smears  the  pressure  variation  and 
often  produces  inaccurate  solutions.  Following  the  approach  in  Shuen  et.  al. 
and  Chen,^  these  difficulties  are  circumvented  in  the  ALLSPD-3D  code  by  us- 

.  The  re¬ 
sults  of  these  treatments  show  that  the  convergence  characteristics  of  the  code 
can  indeed  formulated  to  be  independent  of  the  flow  Mach  number.^ 

Chemically  reacting  flows  exhibit  another  class  of  numerical  difficulties 
because  of  the  wide  range  of  time  scales  involved  in  the  calculation.  To  avoid 
this  chemical  stiffness  problem,  the  chemical  source  terms  are  treated  fully  im¬ 
plicitly.  This  measure  is  equivalent  to  preconditioning  the  time  derivative  terms 
of  the  species  conservation  equations  so  that  all  chemical  and  convective  pro¬ 
cesses  proceed  at  approximately  the  same  numerical  rate. 

As  pointed  out  in  Shuen  et  al.,^  the  Navier-Stokes  computation  of  a  gen¬ 
eral  three-dimensional  turbulent  combustion  flow  has  been  a  tremendous  task  in 
the  past  and  even  with  today's  supercomputers  and  high  end  workstation  clus¬ 
ters,  it  is  still  very  difficult  to  accomplish  such  a  demanding  job.  To  understand 
the  magnitude  of  such  computations,  it  is  useful  to  analyze  a  typical  three-di¬ 
mensional  reacting  flow  problem  with  detailed  chemistry.  The  compressible 
Navier-Stokes  equations  are  usually  required  because  of  the  wide  density  vari- 


ations  resulting  from  the  combustion  process  despite  the  low  flow  speed  in  most 
practical  applications.  In  addition,  the  turbulent  flowfield  must  be  modeled  for 
most  cases  with  only  few  exceptions.  Thus,  with  the  standard  k-E  equations  to 
model  the  turbulence,  seven  equations  are  required  to  solve  for  turbulent  non¬ 
reacting  flows.  For  reacting  flows  with  N  species  considered,  the  total  number 
of  equations  to  be  solved  for  a  three-dimensional  turbulent  reacting  flow  is 
7+(N-l ),  and  the  number  of  species  (N)  to  completely  describe  a  typical  hydro¬ 
carbon  fuel  reaction  is  generally  greater  than  40.^®  Coupling  this  number  of 
equations  to  be  solved  with  a  grid  suitable  to  describe  a  typical  combustor  ge¬ 
ometry  (250,000-1-  cells)  quickly  exceeds  the  capabilities  of  existing  computer 
technology.  It  is  definitely  impractical  to  include  so  many  species  and  couple 
with  the  flow  equations  in  terms  of  today's  (and  in  the  near  future’s)  computer 
technology.  In  the  past,  ways  have  been  sought  to  reduce  the  number  of  species 
(reduced  mechanism)  for  some  particular  fuels  in  order  to  form  a  more  tractable 
CFD  problem.^  Unfortunately,  although  significantly  reduced,  the  number  of 
species  is  still  between  10  and  20  depending  on  the  particular  application.  Fur¬ 
ther  reduction  below  this  level  is  unlikely  because  the  reduction  will  also  lead  to 
inadequate  description  of  the  combustion  physics.^® 

The  numerical  issues  relating  to  solving  the  governing  equations  for  turbu¬ 
lent  reacting  flows  are  robustness,  CPU  time  requirements,  implicitness,  stabil¬ 
ity,  equation  coupling,  and  memory.  Some  of  these  issues  are  very  closely 
related,  such  as  robustness,  implicitness,  stability  and  equation  coupling.  During 
the  past  two  decades,  the  development  of  numerical  algorithms  and  the  usage  of 
various  solvers  in  CFD  have  been  more  or  less  dependent  on  the  availability  and 
accessibility  of  computer  resources.  The  trend  has  evolved  from  explicit,  decou¬ 
pled  schemes  to  strongly  implicit  and  fully  coupled  strategies.  The  pros  and  cons 
among  various  schemes  depend  on  the  problems  of  interest  and  sometimes  the 
vectorization  and/or  parallelization  on  computer  platforms  as  well.  Neverthe¬ 
less,  strongly  implicit,  fully  coupled  schemes  generally  give  better  stability,  ro¬ 
bustness  and  convergence,  at  least  for  steady  state  calculations;  however,  they 
usually  take  more  CPU  time  per  iteration  compared  to  explicit  schemes.  For 
combustion  flow  computations,  equation  coupling  (at  least  for  species  equa- 
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tions)  is  ineAdtable  due  to  the  strong  coupling  nature  of  the  chemical  source 
terms  in  a  detailed  chemistry  treatment. 

In  order  to  compromise  between  the  degree  of  coupling  and  the  total  mem¬ 
ory  requirement,  the  partially  decoupled  numerical  scheme  in  developed  by 
Chen  and  Shuen^^  has  been  adopted  for  ALLSPD-3D.  This  procedure  divides 
the  whole  set  of  equations  into  three  groups.  The  first  group  is  the  five  flow 
equations,  the  second  group  is  the  two  turbulent  equations  and  the  third  group  is 
the  species  equations.  Each  group  of  equations  is  solved  in  a  coupled  manner 
and  the  subsequent  sweep  through  the  three  groups  completes  one  iteration.  This 
arrangement  provides  greater  flexibility  in  various  types  of  flow  calculations 
while  maintaining  the  strongly  coupled  and  implicit  manner  within  the  existing 
combustion  code  structure.  The  numerical  issues  relating  to  this  partially  decou¬ 
pled  process  are  described  in  Chen  and  Shuen^^. 

The  solver  in  the  ALLSPD-3D  code  has  been  modified  to  include  the  LU 
scheme  which,  through  our  tests,^  is  thought  to  be  more  suitable  for  solving 
large  sets  of  coupled  three-dimensional  reacting  flow  equations. 

A  recently  developed  low  Reynolds  number  k-z  modeP"^  is  used  to  represent 
the  effects  of  turbulence.  In  this  model,  the  k-z  equations  are  solved  all  the  way 
to  the  wall  and  the  conventional  wall  function  treatment  of  near  wall  regions  is 
not  needed.  Some  variations  of  implementing  this  model,  especially  the  near 
wall  treatment,  are  also  available  as  an  option  to  alleviate  the  fine  grid  require¬ 
ments  near  walls  (which  is  typically  a  heavy  burden  associated  with  this  type  of 
turbulence  model  for  three-dimensional  calculations).  Detailed  treatment  of  the 
turbulent  model  is  discussed  later  in  this  report.  Although  turbulence  is  includ¬ 
ed,  the  mean  flow  quantities  are  still  used  in  the  chemistry  calculations  (laminar 
chemistry)  in  the  current  version  of  the  code.  A  simple  eddy  breakup  model 
(EBM)^^  is  available  as  an  option  to  model  the  turbulence  effect  on  combustion. 

The  liquid  phase  is  treated  by  solving  Lagrangian  equations  of  motion  and 
transport  for  the  life  histories  of  a  statistically  significant  sample  of  individual 
droplets.  This  involves  dividing  the  droplets  into  n  groups  (defined  by  position. 


velocity,  temperature,  and  diameter)  at  the  fuel  nozzle  exit  and  then  computing 
their  subsequent  trajectories  in  the  flow.  The  spray  model  used  in  this  study  is 
based  on  a  dilute  spray  assumption  which  is  valid  in  regions  where  the  droplet 
loading  is  low.^'^  The  liquid  fuel  is  assumed  to  enter  the  combustor  as  a  fully  at¬ 
omized  spray  comprised  of  spherical  droplets.  The  present  model  does  not  ac¬ 
count  for  the  effects  due  to  droplet  breakup  and  coalescence  processes  which 
might  be  significant  in  a  dense  spray  situation.  A  stochastic  process  is  included 
in  the  spray  calculation  to  account  for  the  effect  of  turbulence  dispersion  on  the 
spray  trajectories.  Also  the  temperature  within  each  droplet  is  non-uniform  and 
a  vortex  model^^  is  used  to  compute  the  internal  temperature  distribution. 

Most  combustors  in  aero-propulsive  engines  involve  very  complex  geome¬ 
tries.  To  enhance  geometrical  flexibility,  the  ALLSPD-3D  code  is  written  in 
generalized  curvilinear  coordinates  and  is  equipped  with  multi-block  grid  and 
internal  blockage  capabilities.  These  features  allow  the  users  to  handle  complex 
geometries.  Some  of  the  main  features  of  ALLSPD-3D  are  listed  in  Table  1 
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TABLE  1.  ALLSPD-3D  Main  Features 


•  3D/2D/  axisymmetric  geometries _ _ _ 

•  Finite-difference,  compressible  flow  formulation  with  low  Mach  number 

preconditioning  (current  version  is  intended  only  for  subsonic  flow  simu¬ 
lations  and  uses  central-differencing  for  convective  and  viscous  terms  on 
right  and  left-hand  sides) _ _ _ _ 

•  Strongly  coupled,  implicit  steady-state  algorithm _ 

•  Generalized  finite  rate  chemistry _ 

•  Generalized  property  and  chemical  reactions  database _ 

•  Eddy  breakup  model  for  turbulence-chemistry  closure _ 

•  Low-Reynolds  number  k-e  turbulence  model _ _ 

•  Spray  combustion  (a  stochastic,  separated  flow  model  with  vortex  model 

for  droplet  internal  temperature  distribution) _ 

•  Dilution  and  transpiration  cooling  capabilities _ _ 

•  Serial  &  parallel  versions  (using  PVM  for  message  passing  language)  on 

various  platforms  (BBM,  SGI,  SUN,  HP  workstations,  CRAY)  _ _ 

•  GUI  for  problem  setup  and  post-processing _ _ _ 

•  Various  boundary  conditions  (multiple  inlets/outlets,  dilution  holes, 

transpiration  holes,  periodic,  symmetry,  far-field,  adiabatic  or  isothermal 
walls,  centerline  singularity) _  _ 
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II :  Mathematical 

Formulations 


2.1  Gas-Phase  Equations 

The  three-dimensional,  unsteady,  compressible,  density- weighted,  time-av¬ 
eraged  Navier-Stokes  equations  and  species  transport  equations  for  a  chemical¬ 
ly  reacting  gas  of  N  species  written  in  generalized  non-orthogonal  coordinates 
can  be  expressed  as 

where  the  vectors  Q ,  E,F ,G,E^,F and  Hi  are  defined  as 

q=]q 

£  =  1  (4,(2  +  %yF  +  4,G)  i  +  \yF  +  Ifi) 

f  J  +  +  %F + \G)  Fy  =  \  + 11/ + n/) 

G  =  i  (5,(2  +  q/  +  ?/  +  ?/)  =  i  (5/  +  (;,£  +  ;/) 

In  the  above  expressions,  x,  il  and  4  are  the  time  and  generalized  spatial 
coordinates,  and  and  4,  are  the  grid  speed  terms.  The  variables 
q  ,  q  ,  q  ,  C  *  C  and  4  are  the  metric  terms  and  J  is  the  transformation  Jacobian. 
The  vectors  Q,E,F,G,E  ,F  ,G  ,H  and  H,  in  the  above  definitions  are 

^  ^  V  V  V  c  / 


Q  =  l^p,  pM,  pv,  pw,  p£^p^,pe,  pFi...pFjv-iJ 

E  =  ^p«,  pM^+/7,  pMV,  PMW,  ^PjE^  +  ^jw,  P^M,  pEu,  pMF,...pMyjv_i 
F  =  ^pv,  pMv,  pv^+p,  pvVV,  ^pE^+p^v,  pkv,  pEv,  pvyi...pvy;v_i 
G  =  ^pw,  pMw,  pvw,  pw^  +  p,  ^p£^+pjw,  pfcw,  pew,  pwyi-.pwy^. 
Ey=  (0,  T,,,  «T^  +  VX^  +  WX„  +  q^,  X^,  X^,  ,)  ^ 

=  (O’  V  V  V  “V +  «y’  V  V  ^yv‘iyN- ^ 

<5v  =  (O’ V  ''zz’  ''V ^ 

and  the  source  term  vectors  H^.  and  Hi  are 


=  ^0, 0, 0, 0, 0, 0,  Sg  +  A.  00  j ...  (0^ 


^/  = 

/  .  o  dw 

/  X-'  X""  .  47C  3  P  V , 

y  n  m  ,y  nmu^--:rpr  ”„j:  ’  Z^' 
^  p  p’  jLi  p  p  p  3  P  P  Pat  ^ 

\  p  P  P 


An  3 


H-vi.  y 

p  p  p  3*^P  P  Pdt 


A  -x  dw  2 

y„  ^  -^p  r  n  ^^’0,0, 

^  P  P  P  3  *^P  P  pdt  ■ALi  p  p  fs  p  p 
P  P 

_  v 

y  n„m  ,  0, 0, . .  0,  0 

^  P  P 

P  ^ 

where  p,  p,  u,  v,  w,  k,  e  and  y.represent  the  density,  pressure,  Cartesian  velocity 
components,  turbulent  kinetic  energy,  dissipation  rate  of  turbulent  kinetic 
energy  and  species  mass  fraction,  respectively;  E^e+l  l2{t^W'+vP‘)  is  the  total 
internal  energy  with  6  being  the  thermodynamic  internal  energy ,  and  Oj  is  the 
rate  of  change  of  species  i  due  to  chemical  reactions.  The  co.  in  the  chemical 
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source  term  will  be  described  later  in  Section  2.4.  The  quantities  related  to  the 
source  term  in  the  turbulent  equations  are  described  in  Section  2.6.  The  normal 
and  shear  stresses,  energy,  species,  and  turbulent  diffusion  fluxes  are  given  by 


^  du  2  f 3m  .  dv  ,  dw\ 


dx^ dy^dz  ■ 


^  9v  2..  f  3m  .  3v 

V  -  3^'Ux  a)-  3zJ 

-  dw  2..  r3M.3v  3w^ 

■'zz  -  21^*02  3^«Ux  3y  3zJ 


(6) 


where  T,  |X  ,  |X  and  are  the  temperature,  effective  viscosity,  turbulent  vis¬ 
cosity,  molecular  viscosity  and  the  effective  thermal  conductivity,  respectively. 
The  vector  H,in  Equation  (1)  represents  the  source  term  that  accounts  for  the 
interactions  between  the  gas  and  liquid  phases.  The  is  the  number  of  droplets 
in  thep^  characteristic  group  of  droplets.  The  is  the  evaporation  rate  of 
each  particle  group;  p  is  the  liquid  density.  The  is  the  droplet  radius.  The  u^, 
V  and  w  are  the  particle  velocities  for  the  particle  group.  The  is  the 
enthalpy  of  fuel  vapor  at  the  droplet  surface,  and  A  A  T  is  the  convective  heat 
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transfer  between  two  phases.  Detailed  discussions  about  the  liquid  phase  equa¬ 
tions  will  be  described  in  Section  2.3.  In  reacting  flow  calculations,  the  evalua¬ 
tion  of  thermo-physical  properties  is  of  vital  importance.  In  this  code,  the 
values  of  C  .,  k.,  and  \i.  for  each  species  are  determined  by  fourth-order  polyno¬ 
mials  in  temperature,  as  described  in  Section  2.7  and  Section  2.8. 


2.2  All-Mach-Number  Formulation 

As  noted  earlier,  the  two  main  difficulties  that  render  the  compressible  flow 
algorithms  ineffective  at  low  Mach  numbers  are  the  round-off  error  caused  by 
the  singular  pressure  gradient  term  in  the  momentum  equations  (the  pressure 
term  is  of  order  1/M^  while  the  convective  term  is  of  order  unity  in  the  non-di¬ 
mensional  momentum  equations)  and  the  stiffness  caused  by  the  wide  dispari¬ 
ties  in  eigenvalues.  To  circumvent  these  two  problems  regarding  the  low  Mach 
number  ceilculations,  a  time  preconditioning  term  has  been  applied  to  rescale  the 
system  eigenvalues  and  the  pressure  variable  is  decomposed  into  a  constant  ref¬ 
erence  pressure  (Pjef)  ^  E^nge  pressure  (Pg)-  This  approach  has  been  suc¬ 
cessfully  applied  to  solve  chemical  non-equilibrium  viscous  flows  with  and 
without  sprays  by  Shuen  et.  al.^  and  by  Chen  and  Shuen.”^  After  including  the 
all-Mach-number  formulation,  the  resulting  Navier-Stokes  equations  in  a  con¬ 
servative  form  are 


where  the  primitive  variable  vector  Q  and  the  preconditioned  matrix  are  given 
as 


Q  =  i^Pg>  V,  w,h,k,z,  Yy..Y^_  j 


(8) 
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0  0  0000.... 
p  0  0000.... 

0  p  0000.... 

0  0  p  0  0  0  ...  . 

pii  pv  pw  p  0  0  .  .  .  . 

0  0  0  .  p  .  .  .  .  0 

0  0  0  .  .  p  .  .  .  0 

0  0  0  0  .  .  p  .  .  . 

0  0  0  0  0  .  ’.  p  .  . 

. 0  .  .  p  . 

. 0  .  .  p 

0  0  0  0  0  .  .  0  .  . 

where  x*  is  the  pseudo-time,  h=e+plp  is  the  specific  enthalpy  of  the  gas  mix¬ 
ture.  The  eigenvalues  (in  the  ^-direction)  in  pseudo-time  can  be  obtained  firom 
the  matrix  F‘A,  where  A  is  the  Jacobian  .  The  eigenvalues  for  a  real  gas 

are 


where  c  is  the  speed  of  sound,  =  ^,a^  =  4^,  and  U  is  the  contravari- 

ant  velocity  component  defined  as  U  =  u  +  ol^v+ a^w.  To  obtain  well-condi- 
*  tioned  eigenvalues,  the  scaling  factor  P  is  taken  to  be  P  =  u^+v^+w^;  such  that 

all  the  eigenvalues  will  have  the  same  order  of  magnitude.  The  definition  of 


11 


vectors  in  Equation  (7)  is  identical  to  those  in  Equation  (1)  except  that  the  pres¬ 
sure  terms  in  the  momentum  equations  are  replaced  by  the  gauge  pressure. 

2.3  Liquid-Phase  Equations 

2.3.1  Droplet  Motion  Equations 

The  liquid  phase  is  treated  by  solving  Lagrangian  equations  of  motion  and 
transport  for  the  life  histories  of  a  statistically  significant  sample  of  individual 
droplets.  This  involves  dividing  the  droplets  into  n  groups  (defined  by  position, 
velocity,  temperature  and  diameter)  at  the  fuel  nozzle  exit  and  then  computing 
their  subsequent  trajectories  in  the  flow.  The  spray  model  used  in  this  study  is 
based  on  a  dilute  spray  assumption  which  is  valid  in  regions  where  the  droplet 
loading  is  low.  The  liquid  fuel  is  assumed  to  enter  the  combustor  as  a  fully  at¬ 
omized  spray  comprised  of  spherical  droplets.  The  present  model  does  not  ac¬ 
count  for  the  effects  due  to  droplet  breakup  and  coalescence  processes,  which 
might  be  significant  in  a  dense  spray  situation.  The  Lagrangian  equations  gov¬ 
erning  the  droplet  motion  are 


II 

II 

16  2  ^ 

0  r 

Vp  p 

dyp 

— =  V 

dt  P 

dv^ 

and  = 

dt 

16  2  «  V 

p/p 

(11) 
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16  2  K  V 

P/;, 

where  the  particle  Reynolds  number,  Re^,  and  the  drag  coefficient,  C^,  are 

defined  as  ♦ 
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V-. 
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Re„ 

1  +  -4 


2/3  \ 


(12) 


0.44 


for  Re^  <  1000 
for  Re^SlOOO 


The  subscript  g  represents  the  gas-phase  quantities  andp  represents  the  liquid- 
phase  (or  “particle”)  quantities.  Equation  (11)  and  Equation  (12)  are  used  to 
calculate  the  new  droplet  positions,  and  the  new  droplet  velocities  at  the  new 
locations.  A  second-order  Runge-Kutta  scheme  was  used  to  integrate  Equation 
(11). 


23.2  Droplet  Mass  and  Heat  Transfer  Equations 

As  described  by  Faeth,^"^  the  following  correlations  were  used  to  approxi¬ 
mate  the  mass  and  heat  transfer  coefficients  for  a  single  isolated  spherical  drop¬ 
let. 


> 


=  NA\+B) 


pDf 


(13) 


^  (14) 

where  h  is  the  heat  transfer  coefficient,  is  the  droplet  diameter,  D^is  the  fuel 
mass  diffiisivity,  k  is  the  thermal  conductivity  of  fuel  v^or  and  m"  is  the  fuel 
mass  evaporation  rate  per  unit  area.  The  and  are  defined  as. 


1  + 


0.276^.>-^_  ^ 


Q21(>Re\^'^Sc'^ 


1  + 


1.232 


Re/r 


3/4 


1/2 


1  + 


1.232 


1 1/2 


(15) 


RCpSc 


3/4 


13 


where  Sc  and  Le  are  the  Schmidt  and  Lewis  numbers,  respectively.  The  Spald¬ 
ing  number,  B,  in  Equation  (13)  and  Equation  (14)  is  defined  as. 


B  = 


1-L 


fgP 


(16) 


where  Y  is  the  fuel  vapor  mass  fraction  at  the  surface  of  the  droplet  and  is 
the  mean  fuel  mass  fraction  of  the  ambient  gas.  In  the  present  study,  Y^^^  is 
obtained  from  the  following  equation. 


(1 


(17) 


where  W  is  the  molecular  weight  of  gas  nuxture  excluding  fuel  vapor,  W^is  the 
molecule  weight  of  fuel  and  is  the  mole  fraction  of  fuel.  The  is 
obtained  from  the  assumption  of  Raoult’s  law.  Based  on  this  assumption,  the 
mole  fraction  at  the  droplet  surface  is  equal  to  the  ratio  of  the  partial  pressure  of 
fuel  vapor  to  the  total  pressure.  For  the  present  spray  calculation,  the  partial 
pressure  of  fuel  viq)or  is  computed  using  the  empirical  correlation  in  Reid.^^ 


+  ^v2'* 


1.5 


+  /*, 


v3- 


jc^  +  F, 


v4-^ 


(18) 


where  x=l-T/r,  T  and  p  are  the  critical  temperature  and  critical  pressure  of 
the  fuel  vapor,  respectively.  The  constants,  P^^  are  obtained 

from  Reid. 


2.33  Droplet  Internal  Temperature  Equations 

As  a  single  droplet  enters  a  hot  environment,  the  small  portion  of  the  droplet 
near  the  surface  will  heat  up  quickly  while  the  center  core  of  the  droplet  remains 
“cold.”  The  heat  will  be  conducted  and  convected  to  the  droplet’s  interior  as  the 
droplet  penetrates  further  into  the  hot  ambient  gas.  Eventually,  the  temperature 
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within  the  droplet  will  become  nearly  uniform  before  the  end  of  its  lifetime.  To 
solve  this  transient  process  occurring  within  the  droplet  is  not  trivial.  In  the  past, 
certain  approximations  were  usually  made  to  alleviate  this  computational  bur¬ 
den  while  retaining  the  solution  integrity.  The  simplest  approximation  is  the  uni¬ 
form  temperature  noodel.  This  model  assumes  that  the  thermal  conductivity  of 
the  fuel  is  infinite.  Of  course,  this  is  not  valid  at  the  beginning  of  the  heating-up 
process  of  the  droplet.  Another  model  considers  the  heat  diffusion  inside  the 
droplet,  yielding  a  conduction  model.  The  temperature  distribution  within  the 
droplet  is  obtained  by  solving  the  one-dimensional  heat  conduction  equation 
subject  to  the  convective  boundary  conditions  at  the  droplet  surface.  The  con¬ 
duction  model  completely  neglects  the  convective  phenomena  within  the  drop¬ 
let  which  mighf  occur  due  to  significant  shear  forces  at  the  surface  induced  by 
high  slip  velocities.  Tong  and  Sirignano  developed  a  vortex  model  for  the  in¬ 
ternal  temperature  of  a  single  droplet  which  accounts  for  the  convective  effect 
of  the  Hill’s  vortex  formed  inside  the  droplet.  Tong  and  Sirignano’s  vortex  mod¬ 
el  is  used  in  ALLSPD-3D  code  to  obtain  the  internal  temperature  distribution  of 
the  droplet.  The  equations  governing  the  internal  temperature  distribution  based 
on  this  model  are. 


ar,  i7fc, 

where  the  value  of  a  varies  between  0  and  1  with  0=0  referring  to  the  vortex 
center  and  0=1  referring  to  the  droplet  surface.  The  initial  and  boundary  condi¬ 
tions  for  Equation  (19)  are 
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inj 


a  =  0 


a  =  1 


da  m  ki  Jypdt  J 


da 


(20) 
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where  di^/da  is  obtained  from  the  energy  balance  at  the  droplet  surface, 


(21) 


where  m],  and  h  are  calculated  from  Equation  (13)  and  Equation  (14),  respec¬ 
tively,  is  the  latent  heat  of  the  fuel  and  A  T  =  where  is  droplet  sur¬ 

face  temperature  and  fg  is  the  mean  gas  temperature  evaluated  as. 


f  = -T  +-T 

*  3  «  3  p* 


(22) 


An  implicit  scheme  is  used  to  solve  Equation  (19)  subject  to  the  initial  and 
boundary  conditions.  Second-order  central  difference  is  used  for  the  spatial  dif¬ 
ferential  terms  and  a  first-order  difference  is  used  for  the  time  term.  These  treat¬ 
ments  render  a  scalar  tridiagonal  algebraic  system  and  is  solved  by  the  Thomas 
algorithm.^* 


2.4  Chemical  Source  Term 

For  a  set  of  elementary  reactions  involving  N  species,  the  reaction  equations 
can  be  written  in  the  following  general  form. 
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v’.,C5.  ^ 

ik  i  -> 
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i  =  1 


,k=  1,2,3  ...Nr 


(23) 


where  and  V”^  are  the  stoichiometric  coefficients  for  species  i  appearing  as 
a  reactant  in  the  ife*  forward  and  backward  reactions,  respectively.  The  law  of 
mass  action  states  that  the  rate  of  change  of  concentration  of  the  species  is 
defined  by  the  movement  of  each  species  summed  over  all  of  the  contributing 
reactions.  Therefore,  the  rate  of  change  in  species  concentration  is  expressed  as 
follows. 
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The  forward  and  backward  reaction  rate  constants  {Kji  K^,)  are  computed 
from  Arrhenius’  relationship  with  phenomenologically  developed  forms.  The 
forward  and  backward  reaction  rate  constants  of  the  reactions  are  defined  as. 


K^(k)=Af^(k)T  exp 


,  FT", 

y  u  y 


ctfik) 


'EAk)y 

Kj,(k)  =  Ajrik)T  '  exp^-^^ 


(25) 


“it  = 


(26) 


The  influence  of  other  bodies  on  a  particular  chemical  reaction  is  modeled 
using  a  typical  third  body  efficiency  correction, 


N 

K^{k)  =  K^{k)  Y,  Third.  (fc)C. 

'  =  1  (27) 

N 

K^ik)  =  Kj^ik)  Y  Third.  WC. 
i  =  1 

It  is  noted  that,  because  of  the  vastly  different  time  scales  that  may  be  involved 
in  the  elementary  reactions  and  the  flow,  and  the  strong  dependence  of  the 
source  terms  on  temperature  (exponentially)  and  density  (p^  or  p^,  depending 
on  the  order  of  reaction),  the  set  of  equations  may  become  very  stiff  for  most 
flow  conditions  of  interest.  To  mitigate  this  problem,  chemical  source  terms  are 
treated  in  an  implicit  fashion  in  the  ALLSPD-3D  code  to  improve  numerical 
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stability.  Details  for  the  constraction  of  the  source  term  Jacobian  are  described 
in  Appendix  A. 


2.5  Equilibrium  Rate  Kinetics  Model 

There  is  a  large  database  of  finite-rate  chemistry  data  where  only  one  of  the 
two  required  reaction  rate  constants  is  reported.  The  chemical  rate  information 
in  the  secondary  direction,  either  forward  or  backward  direction,  can  be  deter¬ 
mined  from  an  equilibriiun  assumption.  Typically,  the  equilibrium  constant  (Xp) 
for  a  reaction  can  be  found  in  many  standard  chemistry  references  as  a  function 
of  temperature.  However,  it  is  convenient  to  deternune  Xp  as  a  part  of  the  solu¬ 
tion.  This  can  be  done  by  taking  advantage  of  the  large  thermodynamic  infor¬ 
mation  available  in  the  NASA  Lewis  database.  The  equilibrium  rate 
information  can  be  determined  for  a  large  number  of  hydrocarbon  related  spe¬ 
cies  using  the  following  procedure. 

In  an  equilibrium  state,  it  is  known  that  the  partial  pressure  ratio  of  the  spe¬ 
cies  component  of  an  ideal  mixture  is  equal  to  the  equilibrium  constant  (Xp).^^ 
In  equilibrium  dOJi- p  is  zero  and  therefore  only  the  Gibbs  free  energy  of  forma¬ 
tion  for  each  species  is  needed  to  determine  the  equilibrium  constant.  For  a  sim¬ 
ple  reaction,  a  CS(a)+  b  CS(b)  <5  c  CS(c)+d  CS(d),  the  equiUbrium  constant  in 
terms  of  partial  pressures  (nondimensionalized  by  atmospheric  pressure)  is  de- 
fmed  as. 


d  b 

pg]  ra 

c  d 

Pc] 


exp 


R  T 
u 


(28) 


where  AG®  is  defined  as  standard  free  energy  change.  This  relationship  can  be 
extended  to  any  general  chemical  reaction  using  the  standard  free  energy 
change  expressed  in  terms  of  molar  Gibbs  free  energy  [kcal/kg-kmole].  The 
molar  Gibbs  free  energy  is  defined  by  the  following  expression. 
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n=l 


n=  1 


where  the  partial  molar  Gibbs  free  energy  of  formulation  of  n*  species,  » 
consistent  with  the  polynomial  representations  of  Gordon  and  McBride  is. 


2  2 

AgliT)  =h-Ts  =  C^  (0)il-ln(T))T-Cp 
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(32) 


The  equilibrium  constant  in  terms  of  partial  pressures  is  converted  to  the  mo¬ 
lar  based  constant  required  for  Equation  (24)  using  the  following  expression. 


=  SA 


(v) 


Sf= 


CR. 

1  ir 


(33) 


Once  the  forward  reaction  rate  constants  are  known,  the  backward  reaction 
rate  constant  is. 


2.6  Two  Equation  Model  of  Turbulence 

In  general,  there  are  several  types  of  low-Reynolds  number  k-£  models  and 
the  major  distinguishing  feature  of  these  models  is  the  handling  of  the  no-slip 
wall  singularity  that  exists  in  the  destruction  term  of  the  specific  dissipation 
equation.  The  most  common  practice  is  to  use  a  modified  form  of  the  transport 
equations,  where  a  fictitious  dissipation(i)  rather  than  the  real  dissipation(E)  is 
used.  The  fictitious  dissipation  is  usually  defined  as  £  =  £  ■  Swall*  There  are  three 
types  of  low-Reynolds  number  k-£  models.  The  first  type  uses  flowfield  charac¬ 
teristics,  such  as  curvature  of  the  turbulent  kinetic  energy(£vyaii  =  2v(d  l^Jdy  )  ), 
to  satisfy  the  near  wall  asymptotic  behavior,  while  at  same  time  avoiding  the 
singularity  caused  by  the  no-slip  boundary  condition.  The  early  k-£  model  for¬ 
mulation  proposed  by  Launder  and  Sharma^^  is  typical  of  this  type  of  formula¬ 
tion.  The  second  type  of  low-Reynolds  number  k-£  model,  proposed  by  Chien, 
uses  the  wall  normal  distance(£^all  “  2vk/yn^)  and  Van-Driest  t5q)e  damping  to 
develop  the  necessary  near  wall  characteristics.  The  third  type  of  low-Reynolds 
number  formulation  makes  no  attempt  to  modify  the  transport  equation;  rather, 
finite  near-wall  characteristics  are  used  as  boundary  conditions.  The  low  Rey¬ 
nolds  number  model  proposed  by  Shih^^is  typical  of  the  models  that  avoid  the 
zero  boundary  condition  and  therefore  the  singularity.  All  three  employ  addi¬ 
tional  wall  damping  to  satisfy  the  near  wall  physics. 

In  practice,  the  low-Reynolds  number  k-£  models  are  easier  to  implement 
numerically  than  high  Reynolds  number  models  with  wall  function  boundary 
conditions  for  complex  geometries  because  of  the  simple  boundary  conditions. 
However,  all  of  the  turbulence  models  have  limitations  imposed  by  the  numer¬ 
ical  implementation  and  physics.^*  Therefore,  the  two-equation  model  solver  in 
ALLSPD-3D  has  been  designed  to  be  generic  and  modular. 

Presentiy,  the  low-Reynolds  number  model  developed  by  Shih,  et  al.  of 
CMOTT  has  been  incorporated  into  ALLSPD-3D  as  the  default  turbulence 
model.  However,  the  modular  design  of  the  numerical  implementation  allows 
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the  exact  form  of  the  low  Reynolds  number  characteristics  to  be  changed  with 
only  minor  modifications.  Furthermore,  several  additional  numerical  options  are 
incorporated  to  address  the  numerical  stiffness  caused  by  the  high  near-wall  grid 
resolution.  The  typical  form  of  the  source  term  required  for  the  Shih  et.  al.  model 
is  given  by  the  following, 


2 

C  -P  —C  f 
_zlk  k  ^2^8  k  . 

where  the  exact  turbulent  production,  P^,  is. 


(35) 
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Here  a  linear  stress-strain  relationship  is  used  to  model  the  Reynolds’  stress 
components.  The  low-Reynolds  number  corrections  required  for  the  source 
term  are  given  by  the  following  expression. 


A  = 


0 

-2pp,(V2m)^ 


(38) 


The  modeling  constants  fe»  (^el»  and  Ce2  are  summarized  in  Table  2  . 


Model 

Cel 

0^2 

Ok 

Oe 

Shih  et  al. 

1.44 

1.92 

0.09 

1.00 

1.30 

l-0.22exp(-Rt2/36) 

TABLE  2.  Turbulence  Model  Parameters 
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The  eddy  viscosity  is  defined  by 

2 

3  9 

and  the  modeled  damping  coefficients  are  a^  =  -1.7x10'  ,  a3  =  -1x10  ,  and  a5= 
-5x10'^®. 

The  eddy  viscosity  constant  can  also  be  varied  as  a  function  of  the  mean 
shear  as  suggested  by  Shih.  The  formulation  of  a  variable  eddy  viscosity  con¬ 
stant  has  been  developed  to  satisfy  reahzabihty  criteria  and  has  been  shown  to 
be  more  accurate  for  non-homogeneous  flows  and  wall  bounded  flow  situa- 
tions,^'^ 


Cjj  =  min 
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2.7  Transport  Property  Model 

All  of  the  transport  properties  used  in  ALLSPD-3D  have  been  derived  from 
the  database  of  the  physical  property  information  developed  by  Gordon  and 
McBride.  However,  in  order  to  reduce  complexity,  property  curves  for  the  dif¬ 
fusion  and  conductivity  coefficients  are  converted  into  single  forth  order  poly¬ 
nomial  curves.  The  typical  diffusion  property  curves  for  the  i*  species  are 
defined  by  the  following  expressions, 

4  4 

;fc=o  '  fc=o 

The  C  (ifc)  constants  have  the  dimensional  property  of  [kg/(m-s)/K*],  and 
C  (k)  constants  have  the  dimensional  property  of  [Joules/(m-s-K)/K  ].  Once 
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the  individual  species  viscosities  are  found,  the  mixture  viscosity  and  conduc- 
tivity  are  computed  using  Wilke’s  mixing  rule, 


, •  =  1  V  j  =  1 


'  I  i  I  I  i:  I  % 

-  t  .  -f  ,  *  •  1 


i  =  1  V  j  =  1 


Here,  the  inter-collisional  parameter((t))  between  species  i  and  j  is  defined  as 


The  inter-species  diffusion,  between  the  species  i  and  j,  is  modeled  using  the 
Chapman-Enskog^  formula  for  dilute  gases. 


where  the  effective  collision  integral  factor,  Qjj,  and  the  effective  temperature, 
Tj),  are  approximated  by 


7^145 + 


(45) 


Once,  the  individual  species  diffusivity  properties  have  been  determined,  the 
overall  diffusion  coefficients  are  computed  by  lumping  all  of  the  species  togeth¬ 
er  into  a  mixture  and  then  computing  the  binary  diffusion  between  the  i**^  species 
and  the  mixture.  An  approximate  form  of  the  diffusion  coefficient  is  as  follows. 


(46) 


Once,  the  binary  diffusion  coefficients,  Dij,  for  all  of  the  species  combinations 
are  known,  the  diffusion  velocity  of  each  species  is  found  using  Pick’s  law.  The 
species  diffusion  velocity  components,  7. J  ,  are  defined  as 


(47) 


The  turbulent  diffusion  of  momentum,  mass  and  heat  are  modeled  based  on 
the  assumption  that  the  eddy  diffusivity  is  directly  proportional  to  the  turbulent 
diffusivity  by  a  constant,  which  is  the  turbulent  Schmidt  number  (Sc^)  for  mass 
diffusion  and  the  turbulent  Prandtl  number  (ft-,)  for  thermal  diffusion.  Once 
combined  with  the  mean  properties,  effective  diffusion  coefficients  can  be  ex¬ 
tracted  and  defined  as  follows, 


|X 

’^eff  "  Prt 

The  turbulent  Schmidt  number  and  the  turbulent  Prandtl  number  are  assumed 
to  be  equal  to  0.90. 


2.8  Thermodynamic  Model 

The  mass-averaged  constant  pressure  specific  heat  constant  ( )  is  as¬ 
sumed  to  be  a  fourth  order  polynomial  function  of  temperature.  The  curve  fitted 
constants  have  been  adopted  from  Gordon  and  McBride  Once  the  individual 
species  specific  heat  ( )  have  been  found,  the  mixture  specific  heat  (C^^ )  can 


be  found  by  weighting  over  all  of  the  species. 


c„=icp.W7* 

jfc  =  0  .  =  1 

The  constants  have  the  dimensional  property  of  [Joules/kg/K*].  The  internal 
energy  of  the  gas  mixture,  found  by  using  mixture  properties  such  as  pressure, 
density  and  enthalpy,  is  defined  as 
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=  1 

1  =  1 


h.  =  hr 

I  h 


j  Cpx 
Tref 


where  ftg  is  an  integration  constant  to  satisfy  the  reference  temperature  condi¬ 
tion  of  Gordon  and  McBride.^^  Since  the  specific  heats,  Cp,,  are  formulated  as  a 
function  of  temperature  only,  the  enthalpy  relationship  can  be  easily  integrated 
to  yield  a  polynomial  relationship  between  enthalpy  and  temperature  for  the 
gas  mixture.  Newton’s  iteration  is  used  to  find  the  temperature  from  the  pre¬ 
dicted  internal  energy  of  the  gas  mixture.  The  pressure  of  the  gas  mixture  is 
determined  from  Dalton’s  law  of  partial  pressure. 


N  Y. 

P  =  X  A^. 

j  =  1  ' 


2.9  Tbrbulent  Chemistry  Model 

The  influence  of  turbulent  intermittence  on  the  finite  rate  chemistry  is  mod¬ 
eled  using  a  simple  eddy  dissipation  concept.  The  simple  formulation  of  Mag- 
nussen  and  Hjertager^^  is  employed  where  the  intermittency  of  the  turbulent 
flame  is  modeled  in  terms  of  rates  of  dissipation  of  eddies.  Here  the  rate  of  eddy 
dissipation  is  related  to  the  mean  concentrations  of  fuel  species  and  the  turbulent 


time  scales,  where  the  turbulent  time  scale  is  defined  as  the  ratio  of  the  turbulent 
kinetic  energy  to  specific  dissipation.  The  eddy  dissipation  concept  offers  one 
distinct  advantage  over  the  eddy  break  up  model  concept  of  Spalding^^  in  that 
the  modeled  term  is  directly  related  to  the  mean  species  quantities  and  does  not 
require  additional  modeling  of  species  variance  terms(^  ).  This  turbulent 

chemistry  model  is  a  first  order  correction  that  has  been  shown  to  be  accurate  in 
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combustor  applications. 

The  species  source  term  modeling  the  mass  transfer  between  the  species 
shown  in  Equation  (24)  is  replaced  with  the  following. 
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(52) 
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Ill :  Numerical  Methods 


Equation(7)  is  the  final  gas-phase  governing  equation  to  be  solved  numeri¬ 
cally.  As  mentioned  in  the  introduction,  solving  the  coupled  three-dimensional 
equations  numerically  with  “many”  species  using  a  reasonable  number  of  grid 
points  (say,  a  quarter  million)  in  a  fully  coupled  manner  is  impractical.  The  main 
problem  with  using  the  fully  coupled  approach,  like  the  one  in  Chen,et.  al.^  is 
that  it  requires  a  large  amount  of  storage  for  the  left-hand-side  (LHS)  coefficient 
matrix.  Although  great  effort  has  been  expended  to  take  advantage  of  several 
storage  saving  tactics,  this  problem  has  not  been  fully  solved.  As  the  number  of 
species  and  the  grid  density  increase,  this  problem  becomes  more  difficult  to  re¬ 
solve. 

Although  it  is  believed  that  the  fiiUy  coupled,  strongly  implicit  algorithm  is 
the  best  choice  (as  far  as  the  stability,  robustness  and  convergence  are  con¬ 
cerned)  for  stiff  problems  typical  of  combustion  flows,  it  is  not  very  economical 
or  affordable  for  the  practical  complex  problems  found  in  industry.  Other  par¬ 
tially  impUcit  and  decoupled  (sequential)  schemes,  such  as  TEACH,  can  expe¬ 
rience  difficulties  in  handling  reacting  flows,  especially  with  sprays.  There  is 
currently  a  wide  spectrum  of  numerical  schemes,  whose  applicability  depends 
on  the  specific  problem  of  interest.  The  advantages  of  a  particular  scheme  may 
be  the  disadvantages  of  others.  In  the  area  of  three-dimensional  combustion 
flow  simulation,  however,  the  freedom  of  choice  is  very  narrow. 

The  approach  in  ALLSPD-3D  employs  a  fully  coupled,  strongly  implicit  al¬ 
gorithm  with  a  slight  modification  to  allow  more  efficient  computation  for 
three-dimensional  reacting  flows.  The  partially  decoupled  procedure  described 
here  is  designed  to  accommodate  more  species  and  more  grid  points  in  three- 
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Flow  equations  •  5  Eqs. 

< — 1 

K-e  equations  -  2  Eqs. 

_ \ 

1 

Species  equations  -  N-1  Eqs. 

FIGURE  1.  Partially  Decoupled  Procedure 

dimensional  combustion  calculations.  It  is  believed  that  this  approach  represents 
a  reasonable  trade-off  between  various  numerical  approaches  available  for  the 
computation  of  finite-rate  combustion  flows.  As  the  growth  of  computer  tech¬ 
nology  in  terms  of  memory  and  speed  progresses,  the  advantages  of  the  current 
approach  will  become  more  pronounced. 

In  ALLSPD-3D,  the  physical  time  (second  term  in  Equation(7))  term  has 
been  eliminated  and  the  solution  marches  forward  in  pseudo-time  to  steady 
state.  The  fully  coupled  governing  equations,  Equation(7),  are  partially  decou¬ 
pled  into  three  sets  of  equations,  i.e.,  five  flow  equations,  the  K-e  turbulent  equa¬ 
tions  and  (N-l)  species  equations.  The  numerical  scheme  used  in  this  code 
solves  each  set  of  equations  in  a  coupled  manner  and  iteratively.  This  process  is 
depicted  in  Rgure  1  . 

The  governing  equations  for  each  set  are  similar  and  have  the  following  gen¬ 
eral  form: 


.(m)  d  .  df  ^(m)  ~  (m) 


(53) 
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where  m=l  for  flow  equations,  m=2  for  turbulent  equations,  and  /n=3  for  spe¬ 
cies  equations.  For  example, 

2^*^  =  ijiPg’  V,  w,  h)  (54) 

and 
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The  term  that  appears  in  Equation(53)  represents  a  source  term  arising 
from  the  partially  decoupled  procedure  and  is  expressed  as. 


1  d  k 
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3.1  LU  Scheme 


After  linearization  and  applying  a  first-order  time  differencing  to  the  pseudo¬ 
time  term,  Equation(53)  can  be  expressed  in  the  following  form. 


9  (m)  ^  d('")  ^  V  f  ^ 


3^  3^/  v3ti  3ti  titi  ^T| 


where. 


where  the  superscript  p  denotes  the  previous  iteration  level,  and  5^"*^  and 
are  the  Jacobian  matrices  for  chemical  and  turbulent  source  terms,  respectively. 
^  (m)  ^(m)  ^m)  inviscid  term  Jacobians  and  and 

are  the  viscous  Jacobians  and  ‘  =  A^  +  Aq'”^  .  The  expressions  for  these 
Jacobians  are  given  in  Appendix  A.  is  the  fourth  order  dissipation  term 
that  will  be  defined  later.  The  lower-upper  symmetric  Gauss-Seidel  (LU-SGS) 
scheme  described  here  was  originally  developed  by  Yoon  and  Jameson  for 
Euler  and  Navier-Stokes  calculations  without  preconditioning.  With  the  current 
timp.  preconditioning,  the  LU-SGS  scheme  has  to  be  modified  to  incorporate 
the  preconditioning  matrix.  Consider  Equation(53),  for  example.  The  LHS 
inviscid  terms,  which  include  the  preconditioning  matrix,  need  to  be  modified. 
For  example,  in  the  ^  direction. 


where  A^'"^].  A^"^  is  then  split  into  two  parts. 


(59) 


where  I  X  is  the  absolute  value  of  the  maximum  eigenvalue  in  the  ^ 

direction  at  each  grid  point,  a  is  a  parameter  greater  than  1,  and  I  is  the  umt 
matrix.  The  term  ^  in  Equation(59)  is  then  approximated  by  the  follow¬ 

ing  one-sided  operator, 


where 


=  l(A<->±a|X4«r<">) 

The  operators,  d'/d  ^  and  d*/d  ^  in  Equation(60)  are  the  first  order  backward 
and  forward  differencing  operators,  respectively.  Similar  procedure  can  be 
applied  to  r\  and  dCf'^Vd  C-  It  can  be  easily  shown  that  the  effect  of  the 

above  eigenvalue  augmentation  is  equivalent  to  adding  a  second  order  central- 
differencing  dissipation  term  (scaled  by  the  local  maximum  eigenvalues)  to  the 
implicit  operator.  Therefore,  the  discretized  governing  equations  for  LU-SGS 
scheme  can  be  written  as  follows. 


(m) 


d  (m) 

5^; 


OT]  OTl' 


(m) 

'2 


where  the  second  and  fourth  order  dissipation  terms  are  defined  as. 
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where  IX and  IX^^"*^!  are  the  absolute  values  of  the  maximum  eigen- 
values  in  the  Ti  and  ^  directions,  respectively.  The  parameter  a  and  a  are 
used  to  define  the  strength  of  the  dissipation.  Here,  the  eigenvalues  for  the  flow 
equations  are  identical  to  those  previously  defined.  The  eigenvalues  for  turbu¬ 
lent  and  species  equations  are  simply  the  contravariant  velocity  components. 
The  LU-SGS  scheme  utilizes  the  following  factorization  procedure  to  obtain 
the  solution. 


MAQ  =  b 

(61) 

iL  +  D)D'\u  +  D)AQ  =  b 

(62) 

Equation(62)  can  be  solved  in  the  following  steps. 

Step  1  (L  +  £>)  AQ  =  b 
Step  2  {U  +  D)^Q  =  d(  Ag*) 

Step  3  A^**  =  A^  +AQ 

where  M  and  b  are  the  LHS  and  RHS  of  Equation(53),  respectively.  The  L,  U 
and  D  are  the  non-zero,  off-diagonal  lower,  upper  and  diagonal  part  of  matrix 
M,  respectively.  The  LU-SGS  scheme  can  be  vectorized  along  a  constant 
(i+j-»-k)  oblique  planes  (see  Figure  2  ).  Therefore,  the  nx  n  LHS  coefficient 
blocks  are  stored  only  for  the  two-dimensional  oblique  plane.  Here,  n  is  the 
number  of  the  equations  to  be  solved. 


33 


k 


FIGURE  2.  Vectorized  Oblique  Plane 


3.2  Mass  Conservation 

The  first  set  of  equations  in  the  partially  decoupled  procedure  is  straightfor¬ 
ward  to  solve.  There  are  no  numerical  differences  between  this  set  of  equations 
and  those  in  the  fully  coupled  approach^  except  that  only  five  equations  are 
solved.  However,  when  the  turbulent  equations  or  species  equations  are  solved, 
there  is  an  extra  source  term,  i.e.,  bf  .  These  terms  consist  of  a  common 
factor,  i.e,  [d  /Bt*  J/p ,  which  has  to  be  evaluated  from  the  continuity 
equation  in  order  to  conserve  mass  in  a  numerical  sense.  Since  the  continuity 
equation  can  be  written  as. 

The,  [d  /3x’J/p  term  can  be  evaluated  as: 
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•(p^)+X"- 


p  p 


(63) 


3.3  Spray  Source  Terms 

In  Equation(53),  the  liquid-phase  interaction  with  the  Navier-Stokes  equa¬ 
tions  is  modeled  as  a  source  term  which  accounts  for  mass,  momentum  and  en¬ 
ergy  exchanges  between  the  two  phases.  In  order  to  compute  this  source  term, 
the  liquid-phase  governing  equations  described  in  Section  2.3  have  to  be  inte¬ 
grated  in  time  along  each  particle  trajectory  from  a  specified  initial  condition. 
As  a  liquid  droplet  begins  its  journey  from  the  injection  nozzle,  its  position,  ve¬ 
locity,  temperature,  and  size  will  change  according  to  its  interaction  with  the 
gas-phase  flowfield.  In  contrast  to  the  Eulerian  approach  for  the  gas-phase  equa¬ 
tions,  the  Lagrangian  treatment  for  the  liquid-phase  equations  requires  interpo¬ 
lation  of  the  flow  quantities  from  the  Eulerian  grid  to  the  particle  positions  and 
re-distribution  of  the  spray  source  terms  from  the  computed  particle  positions  to 
the  Eulerian  grid  for  the  gas-phase  calculation.  The  overall  computational  time 
step  is  carefully  chosen  since  the  time  step  required  for  the  liquid  phase  calcula¬ 
tions  is  much  smaller  than  that  of  the  flow  solver  and  because  the  present  LU 
scheme  allows  the  use  of  a  large  CFL  number.  A  large  number  of  particle  groups 
are  required  to  statistically  represent  the  spray  behavior.  Therefore,  it  is  very  im¬ 
portant  to  keep  the  numerical  efficiency  of  the  present  implicit  scheme  from  be¬ 
ing  severely  degraded  by  the  coupling  (or  interaction)  between  the  gas  and 
liquid  phases.  A  strategy  used  in  ALLSPD-3D  code  to  overcome  this  prob- 
lem,^®in  steady  state  flows,  is  described  in  the  following  subsections. 

33.1  Stochastic  Process 

In  the  present  code,  the  stochastic  separated  flow  (SSF)  model  is  applied  to 
account  for  the  dispersion  effects  on  spray  characteristics.  In  this  model,  the  gas- 

phase  velocity  fluctuations,  u  ’  v  and  w  ’,are  generated  by  randomly  sampling 

s  s  s 
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a  Gaussian  probability  density  distribution  having  a  standard  deviation  of 
V(2ik)/3 .  The  instantaneous  velocities  are  then  used  to  evaluate  Equation(l  1).  A 
fixed  number  of  samplings  is  conducted  for  each  group  of  particles.  The  final 
spray  quantities  are  obtained  by  averaging  the  results  (source  terms,  trajecto¬ 
ries,..., etc.)  among  the  total  samples.  Details  of  the  SSF  model  can  be  found  in 
Faeth^"^  and  Shuen.^^ 

33,2  Determination  of  Spray  Time  Step 

For  spray  combustion  calculations,  several  time  scales  exist  in  the  flow  field 
which  can  differ  by  several  orders  of  magnitude.  The  chemical  reaction  time 
scale  is  usually  very  small  compared  to  the  rate  of  evolution  of  the  gas  flow.  This 
is  also  true  for  spray  evaporation.  To  accurately  calculate  particle  trajectories, 
size  and  temperature,  the  integration  time  step  has  to  be  small.  This  restriction 
is  especially  severe  as  the  droplet  becomes  smaller  and  smaller  toward  the  end 
of  its  lifetime.  For  an  unsteady  problem,  the  time  step  for  the  entire  system  will 
be  controlled  by  the  smallest  time  scale.  For  a  steady  state  calculation,  however, 
the  time  step  for  chemical  reaction  is  usually  not  a  problem  in  the  present  for¬ 
mulation  due  to  the  implicit  treatment  of  the  chemical  source  term.  Experiences 
indicates  ihaf  for  single-phase  combustion  calculations,  the  same  CFL  number 
usually  can  be  used  for  both  non-reacting  and  reacting  calculations  with  the 
present  LU  method.  This  ensures  that  the  convergence  properties  for  reacting 
flow  calculations  do  not  deteriorate  using  the  present  numerical  algorithm. 
However,  the  time  step  for  spray  equations  still  remains  small  and  has  to  be  se¬ 
lected  (computed)  very  carefully  in  order  to  obtain  accurate  spray  results  and  to 
Tpaintain  stabihty.  Here  we  describe  the  determination  of  the  spray  time  step. 
The  spray  interactions  with  the  gas  flow  are  described  in  the  next  section. 

The  spray  time  step  at  any  instant  along  its  trajectory  is  determined  from  the 
following  time  step  constraints:  (1)  droplet  velocity  relaxation  time  (t^,  (2) 
droplet  life  time  (t^),  (3)  droplet  surface  temperature  constraint  time  (r^,  (4)  local 
grid  time  scale  (O  and  (5)  turbulent  eddy-droplet  interaction  time  (t  ).  The  final 
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spray  time  step  (A  ^  is  determined  by  taking  the  minimum  of  t^,  and  t. 

to  ensure  the  accuracy  and  stability  of  the  spray  calculations.  A  factor  between 
0.1  and  0.5  is  used  to  multiply  the  selected  time  step  in  the  current  spray  calcu¬ 
lation.  These  time  steps  are  described  in  more  detail  below. 

Droplet  Velocity  Relaxation  time  (t  ) 

The  locally  linearized  droplet  equations  of  motion,  Equation(l  1),  have  an 
exact  solution  in  terms  of  the  local  slip  velocity  with  an  exponential  decay  form. 
The  rimft  constant  for  the  exact  solution  can  be  expressed  as. 


Droplet  Life  Time  (t^ 

To  ensure  that  the  drop  size  remains  positive  for  the  purpose  of  practical 
computations,  the  droplet  lifetime  at  any  instant  of  time  is  estimated  by  the  fol¬ 
lowing  equation. 


t 


(65) 


Droplet  Surface  Temperature  Constraint  Timeit  1 

When  the  governing  equation  for  the  droplet  internal  temperature  distribu¬ 
tion,  Equation(19),  is  solved,  the  temperature  solution  can  become  completely 
incorrect  due  to  the  use  of  an  inappropriately  large  time  step.  This  is  particularly 
important  for  fuel  with  a  low  boiling  temperature  (close  to  room  temperature). 
In  order  to  ensure  the  success  of  the  temperature  calculation  using  Equation(19), 
the  exact  solution  of  the  infinite  conductivity  model  (also  called  the  uniform 
temperature  model)  is  used  to  estimate  the  time  step  for  the  present  vortex  mod¬ 
el.  The  temperature  equation  for  the  droplet  based  on  this  model  is 


37 


dT^ 


-[^r-rp 


dt  Pfydp 

Equation(66)  has  a  exact  solution  (after  local  linearization)  of  the  form. 


(66) 


ATp^  (67) 

where  A  T  is  the  droplet  temperature  change  within  the  integration  time  step, 
and  A,  A’  and  B’  are  defined  as, 


A'  = 
B’  = 


Pl^v^p 
6 


— -x—rdiZhf^ 
PlCydp  P 


(68) 


Therefore,  a  time  scale  can  be  obtained  based  on  the  above  solution  if  a  desired 
A  is  specified.  This  time  scale  is  expressed  as. 


\ 

1 — ^ 

t  =  _i _ dJ.  (69) 

»  A' 

In  ALLSPD-3D  spray  calculations,  A  T  =3  K  is  specified  for  particles  just  leav¬ 
ing  the  injector,  where  they  experience  a  sudden  temperature  jump,  and 
AT  =  3.0ir  for  the  rest  of  the  spray  calculation  (toward  the  end  of  their  life¬ 
time). 

T.ocal  Grid  Time  Scale  ft  > 

A  particle  can  travel  across  several  grid  cells  and  may  experience  sudden 
changes  in  local  gas  properties  if  the  time  step  is  too  large.  This  not  only  causes 


1 


AT„ 

1 _ E 

A  J 


A' 
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inaccuracy  in  the  integration,  but  also  increases  the  difficulty  of  locating  the  par¬ 
ticle  positions.  Therefore,  a  time  scale,  is  computed  to  ensure  that  the  distance 
a  particle  moves  is  less  than  one  local  cell  size  in  one  time  step. 


Turbulent  Eddy-Droplet  Interaction  Time  (t\ 

According  to  Shuen  et  al,^^  a  particle  is  assumed  to  interact  with  an  eddy  for 
a  time  which  is  the  minimum  of  either  the  eddy  lifetime  or  the  transit  time  re¬ 
quired  for  the  particle  to  cross  the  eddy.  These  times  are  estimated  by  assuming 
that  the  characteristic  size  of  an  eddy  is  the  dissipation  length  scale  given  by. 


Le 


and  the  eddy  lifetime  is  estimated  as. 


(70) 


t 


e 


(71) 


The  transit  time  of  a  particle  was  found  using  the  linearized  equation  of  motion 
for  a  particle  in  a  uniform  flow. 


t.  =  -X  In 


1- 


where 


X  = 


t  tip 


(72) 


and  |jl’ -  Vl  relative  velocity  at  the  start  of  the  interaction.  When 

>  x|d”  -  ,  the  linearized  stopping  distance  of  the  particle  is  smaller  than  the 

characteristic  length  scale  of  the  eddy,  and  then  Equation(72)  has  no  solution. 
In  this  case,  the  eddy  has  captured  the  particle  and  the  interaction  time  is  the 
eddy  lifetime.  Therefore, 


h 


=  if  L,>x|r-Vl 

min  t,)  if  ^  x|il”  -  ii^’j 


39 


333  Interaction  Between  Two  Phases 

For  the  gas-phase  equations,  Equation(53),  the  presence  of  the  spray  appears 
in  the  form  of  a  source  term,  This  source  term  represents  the  interchange 

of  the  mass,  momentum  and  energy  between  the  two  phases.  As  particles  are  in¬ 
jected  into  the  flow  domain,  their  subsequent  behavior  (positions,  velocities, 
size  and  temperature)  is  affected  continuously  by  the  neighboring  gas  properties 
and  vice  versa.  Usually  a  large  number  of  spray  particles  is  desirable  to  accu¬ 
rately  predict  the  spray  behavior.  However,  this  requires  tremendous  computa¬ 
tional  effort  for  the  spray  calculation  alone.  To  minimize  the  computational  time 
for  the  spray,  the  spray  source  term,  need  not  be  updated  at  every  gas- 
phase  iteration  (pseudo-time  marching).  Usually  the  spray  source  term  is  updat¬ 
ed  every  10  to  20  iterations  in  ALLSPD-3D  spray  calculations.  When  the  spray 
soiurce  term  is  updated,  each  group  of  particles  is  integrated  either  to  the  end  of 
its  lifetime  or  imtil  it  leaves  the  computational  domain.  It  should  be  noted  that, 
for  the  present  steady  state  spray  combustion  calculation,  the  spray  time  step  de¬ 
termined  previously  is  independent  of  the  pseudo-time  used  for  the  gas-phase 
equations,  which  is  determined  from  the  maximum  CFL  number  obtainable  ac¬ 
cording  to  the  local  eigenvalues.  Therefore,  this  strategy  maintains  both  the  ef¬ 
ficiency  of  the  flow  solver  and  the  accuracy  of  the  spray  Lagrangian  integration. 
This  treatment  of  the  gas-liquid  interaction  is  different  from  those  of  Raju  and 
Sirignano^^  where  time-accurate  solutions  are  of  their  primary  concern.  In  the 
present  computation,  it  is  assumed  that  when  the  droplets  impinge  on  the  walls 
they  evaporate  completely. 

3.4  Pseudo-Time  Step 

In  Equation(58),  A  r  is  the  pseudo-time  which  is  determined  from  the  eigen¬ 
value  analysis.  For  the  flow  set  of  equations,  the  eigenvalues  are  (defined  previ¬ 
ously), 
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where, 


F- 

+  4p[^/  +  ^/  +  ^/) 

2 

^  c  y 

+4p(ii/+ii/+n/) 

(.4' 

i-C) 

l/„  =  T1  «+Tl  V+T|  W, 
B  'j:  7  'z 

1/ -  =  c «  +  C  v+  C  w, 
C  rx  7  7 


P  =  maj:[p,  C/^,J 

The  pseudo-time,  A  r  is  then  calculated  by  the  following  equation, 

AT=CFU(K^^+X^^+X^^)^'^  (74) 

where  CFL  is  a  user-input  CFL  number  and  Uref  is  a  reference  velocity. 

3.5  Boundary  Conditions 

Boundary  conditions  in  ALLSPD-3D  are  applied  implicidy.  Certain  bound¬ 
ary  conditions  make  use  of  ghost  cells  which  are  invisible  to  the  user.  In  partic¬ 
ular,  grid  interface  boundaries  acquire  information  for  their  ghost  cells  from  the 
neighboring  blocks.  The  grid  connection  between  blocks  is  assumed  to  be  con¬ 
tiguous  and  the  neighboring  volumes  should  be  approximately  the  same  size. 


41 


3.5.1  Wall  Boundary  Condition 

3.5.1.1  Mean  Equations 

The  boundary  condition  applied  to  the  mean  equations  assumes  that  the 
pressure  gradient  is  zero  at  the  walls.  Since  the  grid  may  be  non-orthogonal,  the 
projected  normal  pressure  gradient  is  found  from  the  grid  metrices.  For  exam¬ 
ple,  on  an  i-face 


dp  _  d^dp  dr]dp  dtjdp 
dn  dnd^  dndl]  dnd^ 

A  no-slip  wall  condition  with  either  an  adiabatic  or  an  isothermal  condition 
may  be  chosen.  The  adiabatic  condition  is  imposed  by  setting  the  enthalpy  gra¬ 
dient  to  zero;  similar  to  the  pressure  condition  described  above.  To  specify  a 
wall  temperature,  the  static  enthalpy  is  determined  for  the  desired  wall  temper¬ 
ature.  The  species  mass  fractions  used  in  this  calculation  are  the  lagged  values 
resulting  from  the  previous  species  iteration. 

3.5. 1.2.1  Dilution  Holes 

A  special  boundary  condition  has  been  defined  for  problems  containing  in¬ 
jection  ports.  These  areas  are  assumed  to  be  relatively  small,  and  to  have  a 
known  mass  flow  rate  and  static  temperature.  The  pressure  is  found  from  a  first 
order  extrapolation  from  the  interior  domain.  The  first  step  in  applying  this 
boundary  condition  is  to  determine  the  actual  injection  area.  This  is  done  by  tak¬ 
ing  the  average  of  the  magnitude  of  the  area  normal  to  the  four  comers  for  every 
grid  cell  within  the  prescribed  bounds. 

Next,  the  change  in  velocity  is  found  by  calculating  the  change  in  mass  flow 
rate  required  to  meet  the  imposed  conditions  assuming  that  the  density  is 
unchanged  from  the  previous  iteration. 
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The  velocity  direction  comes  from  the  direction  cosines  prescribed  by  the  user. 
The  boundary  condition  for  the  energy  equation  is  set  by  holding  the  injected 
fluid  to  a  constant  stagnation  enthalpy.  This  is  accomphshed  by  calculating  the 
static  enthalpy  for  the  prescribed  injected  temperature  using  the  thermal  prop¬ 
erties  from  the  previous  iteration. 


£,h  =  MV  + 


This  differs  from  the  upstream  boundary  condition  in  that  the  mass  flow  rate  is 
prescribed  rather  than  the  velocity. 


3.5. 1.3.2  Transpiration  Cooling  (or  Bleed) 

Walls  containing  transpiration  cooling  or  bleed  have  been  handled  in  a  man¬ 
ner  similar  to  the  injection  boundary  condition.  However,  for  this  boundary  con¬ 
dition,  only  the  mass  flow  is  prescribed  as  described  above.  The  conditions  on 
pressure  and  temperature  are  the  same  as  for  a  wall.  In  other  words,  the  energy 
equation  can  be  bounded  by  either  an  adiabatic  wall  or  by  an  isothermal  wall  op¬ 
tion. 


3.5.1.4  liirbulence  Model  Equations 

At  the  wall,  the  turbulent  kinetic  energy  and  turbulent  dissipation  are  bound¬ 
ed  by  non-zero  values  determined  by  Shih,  eLal.  to  maintain  realizability, 

2 

*bc  =  =bc  =  O'^’T 

The  above  constraint  on  e  is  accurate  for  very  near  wall  flows  but  may  be  diffi¬ 
cult  to  satisfy  without  grid  clustering  near  the  walls.  For  this  reason,  an  optional 


boundary  condition  has  been  included  to  relax  the  grid  resolution  requirement  , 

near  walls.  For  this  condition,  the  gradient  of  turbulent  dissipation  given  by 

Equation(75)  is  replaced  with  3e/3Ti  =0.  » 

3.5.1.5  Dilution  Holes 

The  values  for  the  injected  turbulence  field  are  found  from  the  user  specified 
turbulent  intensity 


Ak  = 


VimT.. 


2pA, 


inj 


where  r„  is  the  injected  turbulent  intensity,  and  Ainj  is  the  total  injection  area. 
The  assumed  injected  dissipation  is  then  calculated  from  the  turbulent  kinetic 
energy  and  a  specified  length  scale.  For  a  dilution  hole,  the  assumed  length 
scale  is  the  square  root  of  the  injection  area. 


At  = 


f 


riiT.. 


3 

^2pA 


inj  J 


j^inj 


3.5.1.6  Transpiration  Cooling  (or  Bleed) 

For  this  type  of  boundary  condition,  the  turbulence  quantities  are  taken  from 
the  closest  interior  grid  point 

3.5.2  Species  Continuity  Equations 

At  a  wall  the  gradient  of  the  species  mass  fractions  is  set  to  zero. 
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3  J.3  Dilution  Holes  or  Transpiration  Cooling  (or  Bleed) 

For  these  boundary  conditions  the  composition  of  the  injected  fluid  is  pre¬ 
scribed  by  the  user  and  remains  unchanged. 

3.5.4  Upstream 

At  an  upstream  boundary,  the  stagnation  enthalpy  is  held  constant.  This  is 
accomplished  by  retaining  the  initial  values  for  the  three  velocity  components, 
species  mass  fractions,  and  by  holding  the  inlet  temperature  constant.  Pressure 
is  extrapolated  from  the  interior  points.  The  default  extrapolation  method  is  a 
second  order  difference.  The  turbulence  quantities  are  held  at  their  initial  defi¬ 
nitions  as  are  the  species  mass  flections. 

The  default  upstream  conditions  for  the  turbulent  flow  scalar  quantities  are 
generated  using  following  expressions.  The  velocity  scales  and  turbulent  kinetic 
energy  are  computed  based  on  local  or  reference  velocity  levels  and  turbulence 
intensity  levels. 


k  = 


(76) 


The  specific  dissipation  is  modeled  using  the  following  three  formulations 
for  attached  flow  situations.  The  first  expression  uses  that  the  most  fundamental 
form  of  the  specific  dissipation  where  it  is  expressed  as  a  function  of  velocity 
and  length  scales.  Here,  the  length  scale  is  assumed  to  be  the  Karman’s  mixing 
length  scale  formulation. 


,3/2 

ref 

'i'l 

The  second  formulation  uses  the  experimental  correlation  of  Mellor.  Here, 
the  experimental  data  is  expressed  as  a  function  of  eddy  viscosity  profiles. 
Therefore,  the  following  expression  is  used  to  find  the  specific  dissipation. 
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The  third  method  is  similar  to  Equation(77).  Here,  the  length  scale  definition 
has  been  modified  to  yield  the  following, 


e  = 


(79) 


3.5.5  Downstream 

The  flow  through  the  downstream  boundary  is  assumed  to  remain  subsonic. 
An  exit  pressure  is  imposed  at  the  grid  cell  past  the  exit  plane.  For  multiple  block 
problems,  different  exit  pressures  may  be  set  for  each  block,  however  care  must 
be  taken  if  the  exit  pressures  vary  greatly  from  the  initial  pressure(which  is  set 
to  user  inputted  value  of  Pref  over  the  exit  face).  If  an  exit  pressure  does  require 
a  significant  departure  from  the  initial  value,  a  large  pressure  gradient  at  the  exit 
can  be  avoided  by  slowly  decreasing  the  exit  pressure  over  a  number  of  iteration. 

The  remaining  flow  quantities  are  extrapolated  from  the  interior  flow  do¬ 
main.  A  zero  order  extrapolation  method  is  the  default.  An  expert  user  flag  is 
included  in  subroutine  qplusl.f  to  change  the  order  of  this  extrapolation  if  need¬ 
ed. 


3.5.6  Periodic 

For  the  periodic  boundary  condition,  flow  quantities  are  copied  from  the  op¬ 
posite  face  into  the  ghost  cells,  e.g.,  the  information  for  the  j=l  plane  would 
come  from  the  j=jmax  plane.  This  boundary  condition  requires  only  two  grid 
points  in  the  periodic  direction  and  is  therefore  the  best  choice  for  a  pseudo-two- 
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dimensional  problem. 

3.5.7  Symmetry 

The  symmetry  boundary  condition  copies  the  flow  field  from  the  neighbor¬ 
ing  interior  point  to  the  “ghost  cell”,  e.g.,  the  information  for  the  j=0  plane 
would  come  from  the  j=2  plane.  Since  ALLSPD-3D  requires  a  minimum  of  five 
grid  points  (including  the  ghost  cells),  this  option  requires  the  user  to  supply  a 
minimum  of  three  grid  points. 

3.5.8  Far  Field 

The  same  condition  is  used  as  for  the  downstream  boundary  condition. 

3.5.9  Singularity 

For  axisymmetric  grids,  where  a  plane  collapses  to  a  single  line,  a  special 
boundary  condition  has  been  implemented.  This  boundary  must  be  applied  for 
the  k=l  boundary.  The  solution  domain  will  not  include  the  points  along  the  sin¬ 
gularity,  so  that  the  centerline  values  must  be  found  by  extrapolation  from  the 
interior  points.  A  check  is  made  for  internal  blockages  at  the  centerline  so  that 
the  appropriate  wall  boundary  conditions  (see  above)  can  be  applied  to  the  tur¬ 
bulence  equations. 
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IV :  Input/Output  Files 


4.1  Input/Output 

The  input  and  output  file  structure  is  explored  in  this  section.  Four  user  or 
GUI  developed  input  files  are  required  to  execute  the  ALLSPD-3D  code.  To  run 
the  ALLSPD-3D  code,  users  need  the  grid  file,  grid.dat,  block  link  file, 
gridlintdat,  property  data  file,  prop.dat,  and  the  flow  data  file,  allspd.dat.  A 
typical  grid.dat  read  format  is  shown  in  Figure  3.  For  cases  where  users  would 
like  to  provide  the  initial  conditions  for  the  flow  variables  within  the  whole  do¬ 
main,  the  user-defined  initial  condition  file,  uvwtdat,  has  to  be  generated.  The 
structure  of  these  input  files  will  be  discussed  along  with  several  examples.  The 
gridlink.dat  file  contains  the  grid  size,  connectivity  and  boundary  condition  in¬ 
formation.  A  typical  gridlink  file  is  shown  in  Figure  4.  A  detailed  example  is 
presented  in  Table  3 .  A  typical  allspddat  file  is  shown  in  Figure  5 .  Various  op¬ 
tions  are  summarized  in  Table  3  .  The  ALLSPD-3D  property  input  file, 
prop.dat,  is  arranged  in  three  major  groups:  Thermal  Property  Information,  Fi¬ 
nite  Rate  Chemistry  Information  and  Spray  Property  Information.  They  are  ar¬ 
ranged  in  top-down  groups  as  shown  in  Figure  6.  The  requirements  of  each 
group’s  input  information  is  summarized  later.  For  liquid  spray  problems,  users 
can  provide  their  own  spray  information  data  file,  spray.dat,  instead  of  using  the 
default  spray  particle  distribution  generator. 

For  the  output  files,  the  most  important  files  are  the  PLOT3D  format  solu¬ 
tion  files  which  can  be  examined  using  the  NASA  PLOT3D  plotting  package*^ 
or  other  post-processors  which  accept  this  file  format.  If  a  liquid  spray  is  initi¬ 
ated,  the  particle  trajectories,  temperatures,  sizes  are  also  written  to  a  single  file 
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for  post-processing.  A  restart  file  for  subsequent  runs  may  also  be  generated  and 
is  usually  needed  for  running  large  combustion  problems. 


4.2  Input  Data  Files 


4.2.1  Grid  Data:  griddat 

The  grid  is  generated  separately  by  the  user  and  read  into  the  code  by  the 
subroutine  MESH.  The  default  grid  data  file  name  is  griddat.  The  format  for  the 
grid  file  is  consistent  with  PLOT3D’s  multi-block  format.  It  is  generated  as  a 
fi-ee  formatted  file  and  therefore  is  portable.  There  are  two  formats  available  for 
grid.dat.  One  is  a  non-blanked  file,  and  the  other  is  an  Iblanked.  The  Iblanked 
file  requires  an  extra  set  of  Iblanked  data  (0  or  1  for  each  gnd  point)  in  addition 
to  the  usual  x,  y,  and  z  data.  For  most  users,  a  non-Iblanked  format  is  preferred. 
The  input  flag  for  either  Iblanked  or  non-Iblanked  grid.dat  file  is  in  gridlintdat 
that  will  be  described  later.  The  format  for  griddat  is  listed  as  follows. 

For  non-Iblanked  format; 

readJ{  iunit,  *)nblk 

read(  iunit,  *)(  imax(  n),jmcix(n),kmax(n),n=l,nblk) 
do  n=l,nblk 
read(  iunit,*) 

(((x(i,j,k),i=l,imax(n))J=lJmax(n)),k=],kmajdn)), 
(((y(i,j,k),i=l,imajdn))J=ljrnax(n)),k=l,kmax(n)), 
(((z(i,J,k),i=l,iTnax(n))j=l,jrnax(n)),k=l,knuix(n)) 
end  do 

For  Iblanked  format; 

read(  iunit,  *)nblk 

read(iunit,*)(imax(n),jmax(n),kmax(n),n=l,nblk) 

do  n=l,nblk 
read(  iunit,*) 

(((x(i,j,k),i=l,imax(n))j=ljmax(n)),k=l,kmax(n)), 

(((y(i,j,k),i=J,imax(n))J=l,jmax(n)),k=l,kmax(n)), 

(((z(i,j,k),i=l,imax(n))J=l,jniax(n)),k=l,hnax(n)), 

(((iblk(i,j,k),i=l,irnax(n)),j=l,jniax(n)),k=l,hnax(n)) 

end  do 

FIGURE  3.  Read  file  format  for  grid  Jot 
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4.2.2  Multi-Block  Data:  gridlinkJat 

The  multi-block  connection  information  is  defined  in  gridlinLdat  and  is  read 
by  the  subroutine  MESH.  The  input  variables  are  described  in  Table  3. 

iblank 

^  Number  of  block  information 

33  49  43 

-Q  43  Number  of  cells  per  block 


33 

49 

43 

34 

49 

43 

30 

49 

43 

0.0  0.0 

0.0  0.0 


101325.0  0.0  0.0  0.0  0.0 


.001058e-l 

.002116e-l 


Boundary  Conditions 


Connectivity 


1  39  43  .001058e-l 


20  24  49  49  1  5 
20  24  49  49  20  24 
20  24  49  49  39  43 
15  16  1  1  14: 
15  16  49  49  1  4 


5  .001058e-l 
24  .002116e~l 
43  .001058e-l 
1  43  .00345e-2 
1  43  .00345e-2 


.258  .965 
.258  .965 
.258  .965 
.258  -.965 
.258  -.965 
.258  -.965 


0.0  .233  0.0 
0.0  .233  0.0 
0.0  .233  0.0 
0.0  .233  0.0 
0.0  .233  0.0 
0.0  .233  0.0 


0.0  .76 
0.0  .76 
0.0  .76 
0.0  .76 
0.0  .76 
0.0  .76 


4  12  1  1  8  16  .00826e-l 

4  12  1  1  28  36  .00826e-l 

4  12  49  49  8  16  .00826e-l 

4  12  49  49  28  36  .00826e-l 

23  27  1  1  1  5  .002077e-l 

23  27  1  1  20  24  .003939e-l 

23  27  1  1  39  43  .002077e-l 

23  27  49  49  1  5  .002077e-l 

23  27  49  49  20  24  .003939e-l 

23  27  49  49  39  43  .002077e-l 


325.0 

.0010 

1.00 

.0 

0.0 

0.0 

i  .233  0.0  0.0 

.76 

325.0 

.0010 

1.00 

.0 

0.0 

0.0 

1  .233  0.0  0.0 

.76 

325.0 

.001 

.139 

.990 

0.0 

0.0  . 

233  0.0  0.0  . 

767 

325.0 

.001 

.139 

.990 

0.0 

0.0  . 

233  0.0  0.0  . 

767 

325.0 

.001 

.139  - 

.990 

0.0 

0.0  . 

233  0.0  0.0  . 

767 

325.0 

.001 

.139  - 

.990 

0.0 

0,0  . 

233  0.0  0.0  . 

767 

325.0 

.001 

.087 

,996 

0.0 

0.0 

.233  0.0  0.0 

.767 

325.0 

.001 

.087 

.996 

0.0 

0.0 

,233  0.0  0.0 

.767 

325.0 

.001 

.087 

.996 

0.0 

0.0 

.233  0.0  0.0 

.767 

325.0 

.001 

.087 

-.996 

0.0 

0.0 

.233  0.0  0.0 

.767 

325.0 

.001 

.087 

-.996 

0.0 

0.0 

.233  0.0  0.0 

.767 

325.0 

.001 

.087 

-.996 

0.0 

o 

o 

.233  0.0  0.0 

.767 

Dilution  Information 


FIGURE  4.  Typical  gridlink,dat61e 
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TABLE  3.  Detailed  Information  about  gridUnk.dat  File 


Input  Variable 

Comments 

Iblank 

The  flag  for  non-Iblanked  or  Iblanked  gruLdat 

0  =  non-Iblanked  griddat 

1  =  Iblanked  gruLdat 

(imaxbl(nb)jinaxbl(nb), 

kmaxbl(nb)4ib=l^blk) 

Number  of  grid  points  in  i,  j,  and  k  direction,  respectively,  in  the  vtth  block 

nblk 

Total  number  of  blocks 

least(nb)4west(nb), 

inorth(nb)4south(nb), 

ltop(nbybot(nb),bcin(nb,l) 

bciii(nb^),bdn(nb^), 
bdn(nb,4),bdn(nb^, 
bdn(nb,6),  nb=l^blk 

Block  boundary  control  indices,  specifying  the  boundary  type  in  the  east,  west, 
north,  south,  top,  and  bottom  faces  of  the  nbrii  block. 

The  boundary  types, 

least(nb)4west(nb)4north(nb)4south(nb),ltop(nb)4bot(nb), 

=  1  upstream  (inflow)  (all  six  faces),  =  2  downstream  (outflow)  (all  six  faces) 

=  3  solid  wall  (all  six  faces),  =  4  block  interface  (all  six  faces) 

=  5  symmetry  (only  for  north,  south,  top,  and  bottom  faces) 

=  6  periodic  boundary  (only  for  nortih/south,  and  top/bottom  faces) 

=  7  singularity  (only  for  bottom  face,  i.e.,  k=l) 

=  8  far-field  (for  external  flows)  (all  six  faces) 

bdn(nb,l),bcin(nb^),bciii(nb^),bciii(nb,4),bcm(nb^,bcin(nb,6) 

are  used  to  assign  either  pressure  or  temperature  at  each  faces  depending  on  the 
boundary  types.  Zero  values  are  usually  used  except  for  those  fac^  that  are 
either  wall  or  exit  conditions.  For  wall  boundary  faces,  zero  (0.0)  is  used  to 
assign  an  adiabatic  wall  and  desired  temperature  (®K)  is  used  to  assign  an 
isothermal  wall.  For  exit  faces,  these  are  used  to  specify  the  ambient  pressure  (iW 

Ti?)  that  is  to  remain  constant  The  corresponding  faces  are: 

bcin(nb4)  east  face,  bdn(nb^)  west  face,  bcin(nb^)  north  face,  bcm(nb,4) 
south  face,  bdn(nb^  top  face,  bdn(nb,6)  bottom  face 

lee(nb)4ww(nb)4nn(nb), 
lss(nb)4tt(nb)^bb(nb) , 
nb=l^blk 

Block  boundary  neighboring  control  indices,  indicating  the  identity  of  the 
neighboring  block,  such  as: 

lee(nb)^m  -  nbth  block’s  east  side  borders  mth  block;  zero  if  nbth  block’s  east 
side  is  not  a  block  interface. 

liin(nb)=n  -  nbth  block’s  north  side  borders  nth  block;  zero  if  nbth  block’s  north 
side  is  not  a  block  interface,  etc ... 

intblk 

Internal  blockage  flag 

rs  0  no  internal  blockage  defined;  skip  the  following  input  data  and  move  directly 
to  the  next  one  for  dilution  holes  (nthl). 

:=1  with  internal  blockage(s);  provide  the  following  input  data  for  internal 
blockaee^s)  in  each  block. 

Input  Variable 

Comments 

TABLE  3.  Detailed  Information  dboxA  gridUnkJat  File 


Input  Variable 


Comments 


4.23  Flow  Data:  attspd.dat 

The  flow  data  and  major  control  of  the  program  execution  is  done  in  the  all- 
spddat  and  is  read  by  the  subroutine  INPUT.  Each  input  is  described  in  Table 
4. 


0 

10 

10 

10 

10  0 

1 

1 

8.000000e-01  3.000000e-01 
l.OOOOOOe+6  l.OOOOOOe+6  l.OOOOOOe+6 
1.050000e+00  1.050000e+00  1.050000e+00 
1000 
0.1 
le-13 
1 
1 
3 

62.5  0.0  0.0 

307.0 

1.013250e+05 

.00254 

-1 

0 

150 

50000 

1 

I 

5  31  3  43  3  43 

1.600000e+03 

5.0e+07 

0 

1 

1 

3  3  22  28  19  25 

3  45.0  65.0 

5.5  7.S 
100 

20 

0 

300.0 

4.67e-4 

50.0 

1  2  3  4  5  15  16  17  18  19  11  12  13  14  8 
5 

'CH30H'  O.OOOOOOe+00 
'02'  2.330000€-01 

'C02'  O.OOOOOOe+OO 
'H20’  O.OOOOOOe+OO 
'N2'  7.670000e-01 


FIGURE  5.  Typical  Listing  of  the  alkpd.dat  Input  FUe 
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TABLE  4.  Detailed  Layout  of  allspd.dat  File 


Input  Data 

Input  Variable 

Comments 

0 

iread 

=  0  a  fresh  run 

=  1  a  restart  run 

10 

ntmax 

Total  number  of  iterations  (psuedo-time  steps) 

10 

nplot 

Every  nplot  number  of  iterations  the  plot3d  files  are  generated, 
and  diese  files  are  overwritten  every  next  nplot  number  of 
iterations 

10 

nsave 

Every  nsave  number  of  iterations  the  restart  file  is  generated,  and 
this  file  is  overwritten  every  next  nsave  number  of  iterations 

1,0,0 

iturb,  icmu, 
idisbc 

turbulence  model  switch, 

=  0  laminar  flow 

=  1  turbulent  flow, 

eddy  viscosity  coefficient 

=  0  constant  Cp  in  turbulence  model  (^.09) 

=  1  variable  Cp  in  turbulence  model 

idisbc 

=  0  Dirichlet  wall  boundary  condition  for  e  equation. 

=  1  Neumann  wall  boundary  condition  for  £  equation. 

1 

icomb 

diemistry  switch 
=  0  non-reacting  flow 
=  1  reacting  flow 

1 

iartv 

Explicit  artificial  viscosity  switch 
=  0  no  artificial  damping 
s  1  with  artificial  danq>ing 

8.000e-0l  3.000e- 

yis4, 

vis4  =  4th  order  danq)mg  coefficient 

01 

vis2 

vis2  =  2nd  order  dancing  coefficient 

le6,  Ie6,le6 

CFL 

cfl(l )  =  CFL  number  for  the  flow  equations 
cfl(2)  =  CFL  number  for  the  K-e  equations 
cfl(3)  =  CFL  number  for  the  species  equations 

1.05, 1.05, 1.05 

precfac 

precfac(l)  =  pre-conditioning  factor  for  LU  scheme  for  the  flow 
equations 

precfac(2)  =  pre-conditioning  factor  for  LU  scheme  for  the  K-e 
equations 

precfac(3)  =  pre-conditioning  factor  for  LU  scheme  for  the  species 
equation 

TABLE  4.  Detailed  Layout  of  allspdJat  File 


Input  Data 


Input  Variable  Comments 


For  some  difficult  flows,  a  smaller  CFL  number  initially  is  helpful 

to  start  the  calculation.  The  CFL  number  can  then  be  ramped  up  as 

the  flow  develops.  The  CFL  number  is  calculated  using  the 
following  fonnuJa: 

CFL=CFLxfactau 

where  the  factor  factau  has  initial  value  of  facini  (defined  next) 
and  is  linearly  increased  to  1  over  nswitch  iterations. _ 

The  initial  value  used  in  conjunction  with  the  above  nswitch. 
facini  is  usually  set  to  0.1 .  _ _ 


The  convergence  criterion,  say,  l.Oe-?  for  exanyle _ 

The  turbulence  equations  are  updated  (called)  every  nkecaU  flow 
iterations.  Say,  if  nkecall=10,  the  turbulence  equations  are  updated 
every  10  flow  iterations 

The  species  equations  are  updated  (called)  every  nscallO  flow 
iterations.  Say,  if  nkecall=10,  the  turbulence  equations  are  updated 
every  10  flow  iterations 


Initial  condition  controlling  flag 

=  -1  only  inlet  faces  have  non-zero  flow  specified  by  uin,vin,  and 
win  described  next. 

=  0  uniform  flow  initialization  for  the  entire  domain  except  for  the 
walls.  The  flow  is  defined  by  uin,vin,  and  win  described  next 

=  1  user  provides  the  initial  flow  condition  for  the  entire  domain. 
In  this  case  the  flow  data  defined  here  (uin,vin,  win)  will  be  us^ 
only  as  reference.  The  description  for  tiie  data  file,  uvwt.dat^  is  in 
another  section. 

s=  2  this  is  used  only  for  easy  specification  of  the  swirling  inlet 
condition.  The  file  name  is  the  same  as  in  iuvwt=l  case  th^  is 
uvwtdax.  The  format  of  this  file  for  this  case  will  be  described 
later.  _ _ 

Initial  velocities.  These  are  used  to  assign  initial  flow  velocity  for 
the  entire  domain  or  inlet  only,  depending  on  die  value  of  iuvwt 
above.  They  are  used  to  calculate  the  inlet  flow  Reynolds  number 
also. 

uin  =  X  conqionent  of  velocity  (m/s) 
vin  =  y  coii5)onent  of  velocity  (m/s) 
win  =  z  component  of  velocity  (m/s) _ 

Reference  tenq)erature  (K)  for  the  flow  domain.  TWs  is  used  to 
assign  initial  flow  teix^rature  for  the  entire  domain.  It  is  used  to 
calculate  the  inlet  flow  Reynolds  number  also.  However,  the 
initial  temperature  will  be  overwritten  by  the  temperature  defined 
in  uvwtdat  if  iuvwt=l,or2. _ 

Reference  pressure  (N/m^)  for  the  flow  domain.  This  is  used  to 
assign  initial  flow  pressure  for  the  entire  domain  except  along 
outlet  faces.  It  is  used  to  calculate  the  inlet  flow  Reynolds  number 
also.  However,  the  initial  pressure  will  be  overwritten  by  the 
pressure  defin^  in  uvwt.dat  if  iuvwt=l ,  or  2. 
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TABLE  4.  Detailed  Layout  of  aUspd.dat  File 


Input  Data 

Input  Variable 

Comments 

.00254 

Iref 

A  reference  length  scale(m)  used  to  calculate  the  inlet  Reynolds 
number  and  turbulence  quantities. 

-1 

iebm 

Eddy  break-up  model  for  turbulence-combustion  closure. 

=  0  no  eddy  break-up  model 
=  -1  with  eddy  break-up  model 

0 

iequate 

A  flag  to  compute  the  reverse/forward  reaction  constant  in 
reacting  flow  calculation. 

=  0  skip  the  reverse/forward  reaction  constant  calculation 
=  1  perform  the  reverse  reaction  constant  calculation 
=  -1  perform  the  forward  reaction  constant  calculation 

150 

ntligt 

Ignition  energy  is  added  to  the  gas  phase  energy  equation  starting 
at  the  time  step  of  ntlig;t. 

50000 

iit2igt 

Ignition  energy  is  turned  off  from  the  gas  phase  energy  equation 
starting  after  the  time  step  of  nt2igt. 

1 

noigt 

Number  of  ignition  sources,  (if = 0,  skips  to  next  input  of  toff) 

1,531343343 

iblkigt 

iligt42igt^j^^ 

ibkigt(iio)41igt(no)42igt(no)jligt(no)j2igt(no), 
kligt(no)4^igt(iio)4io=l^oigt:  Specification  of  the  position  of 
the  ignition  spark(s). 

ibkigt(iio):  The  noth  ignition  is  located  in  ibkigt(no)  block. 
iligt(no):  The  starting  /-index  of  the  ignition  spark  location. 
i2igt(no):  The  ending  /-index  of  the  ignition  spark  location. 
jligt(no):  The  starting  ;-index  of  the  ignition  spark  location. 
j2igt(no):  The  ending  y-index  of  the  ignition  spark  location. 
kligt(no):  The  starting  k-index  of  the  ignition  spark  location. 
k2igt(no):  The  ending  Jk-index  of  the  ignition  spark  location. 

1.600000e+03 

tofr 

The  cut-off  temperature  within  the  ignition  region  to  turn  off  the 
energy  sparks. 

5.0C+07 

sparkuser 

The  strength  of  ignition  spark  (J/(m^-s)).  If  zcto,  then  the  default 
strength  in  the  code  is  used. 

0 

nsdata 

=  0  The  spray  d?ta  are  to  be  generated  in  the  following  inputs. 

-1  The  spray  data  are  to  be  provided  through  user-generated  file, 
spray.dat.  The  format  for  spray.dat  is  described  in  another  section. 

1 

ispray 

=  0  Single-phase  flow  calculation,  (no  liquid  phase  equations  are 
solved) 

=  1  Two-phase  flow  calculation. 

If  ispray=0  or  nsdata=l,  skip  to  the  input  for  inplot(i). 

1 

izospr 

Nozzle  zone 

3  3  22  28  1925 

]sprayl4spray2, 

jsprayljspray2, 

ksprayl,kspray2 

The  spray  particles  will  be  issuing  from  die  (i,j,k)  points  between 
isprayl^pray2  points  in  /  direction,  between  jsprayljspray2 
points  in>  direction,  and  between  ksprayl^pray2  points  in  k 
direction. 
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TABLE  4.  Detailed  Layout  of  aUspd.dat  File 


Input  Data 

Input  Variable 

Comments 

3  45.0  65.0 

Dangles,  angiel, 
angle2 

nan^e  =  Number  of  angles  between  (anglel,angle2)  from  which 
the  spray  particles  are  issuing. 

angiel  =  The  lower  angle  (say  10  degree)  for  the  spray  particles  to 
shoot  out.  (The  angle  currently  is  only  in  the  x-z  plane). 

angie2  =  The  upper  angle  (say  45  degree)  for  the  spray  particles  to 
shoot  out.  (The  angle  currently  is  only  in  the  x-z  plane). 

5.5  7,5 

sizel,  size2 

sizel  =  The  smallest  size  of  spray  particles  (in  terms  of  p  m  in 
radius). 

size2  =  The  largest  size  of  spray  particles  (in  terms  of  p  m  in 
radius). 

100 

ntspr 

The  time  step  to  initiate  the  spray  subroutine. 

20 

ncalspr 

The  spray  subroutine  is  called  every  ncalspr  flow  iterations 
(pseudo  time  steps) 

0 

imspr 

Always  set  to  0  (reserved  for  later  version  use) 

300.0 

tspr 

Liquid  spray  temperature(K) 

4.67e-4 

flowispry 

Liquid  spray  mass  flow  rate  (kg/s) 

50.0 

vispr 

Liquid  spray  velocity 

In  spray  calculations,  the  number  of  groups  of  particles  is 
specified  through  the  parameter  statement(alls^.h)  when 
compiling.  From  the  above  specifications,  a  complete  spray  table 
is  established.  A  spray  particle  is  characterized  by  its  position 
(i  j,k),  angle,  size  and  temperature.  The  mass  of  each  spray 
particle  is  flien  calculated  from  the  flowispr  divided  by  the  tot^ 
number  of  particles  (defined  in  the  allspd.h  file).  Based  on  the  size 
information  for  each  particle,  a  finite  number  of  droplets  can  ht 
calculated  in  each  spray  particle.  Therefore,  when  a  spray  particle 
is  mentioned,  it  actually  represents  a  group  of  small  droplets  that 
have  the  same  size,  velocity,  position,  and  temperature. 

1234515  1617 

18  1911  12  13  148 

inpiot 

Plot  file  flags  specify  the  flow  variables  to  be  written  to  the 
plot3d  files.  There  are  three  files  to  be  written  out  for  plotting  and/ 
or  post-processing  purposes.  The  first  file  contains  vector 
quantities,  allspd J^FV,  and  the  second  and  the  third  files  contain 
all  scalar  quantities,  allsp<LPFSl  and  al]spd^FS2.  In  plot3d,  five 
variables  are  ne^ed  in  each  file  for  three-dimensional  problems. 
AIXSPD-3D  users  are  asked  to  specify  15  flow  variables  for  these 
three  output  files.  The  inplot(i)  has  the  follovdng  valid  integer 
values. 

=1  density,  =  2  u  momentum,  =  3  v  momentum,  =  4  w  momentum 
=  5  static  ten5>crature,  =  6  gage  pressure,  =  7  static  pressure 
=  8  turbulent  Cp,  =  9  static  enthalpy,  =  10  total  enthalpy 
=  11  turbulent  kinetic  energy,  =  12  turbulent  dissipation 
=  13  turbulent  eddy  viscosity/laminar  viscosity 
=  14  distance  to  the  walls,  =  i+14  mass  fraction  of  species  i 

5 

nspr 

Number  species  to  read  in  Initial  Condition 
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TABLE  4.  Detailed  Layout  of  allspdJat  File 


Input  Data 

Input  Variable 

Comments 

‘CH30H’ 

O.OOOOOOe+0 

name,  Ysi(l) 

Species  name  and  species  mass  fraction 

‘02’  2.330000e-01 

name,  Ysi(2) 

Species  name  and  species  mass  fraction 

‘C02’ 

O.OOOOOOe+00 

name,  Ysi(3) 

Species  name  and  species  mass  fraction 

‘H20’ 

0.0000006400 

name,  Ysi(4) 

Species  name  and  species  mass  fraction 

‘N2’  7.6700006-01 

name,  Ysi(ns) 

Species  name  and  species  mass  fraction 

4.2.4  User  Defined  Initial  Conditions:  uvwt.dat 

When  iuvwt=l  or  2  in  allspd.dat  is  specified,  the  code  needs  the  uvwt.dat  data 
file.  This  file  is  read  by  subroutine  INIT. 

For  iuvwt=l:  uvwt.dat  file  must  be  an  unformatted  file. 


do  n=l,nblk 

ii^imaxbl(n) 

jj-jmaxbl(n) 

kk-kmaxbl(n) 

do  k=l,kk 

doj=l,jj 

do 

read(9)  pp,  u,  v,  w,  t,  tke,  tep,(ys(nn),  nn=l,  ns) 
end  do 


where 

nblk  =  number  of  blocks 

imaxbl(n)  =  grid  points  in  /  direction  in  nt/i  block 

jmaxbl(n)  =  grid  points  in  j  direction  in  uth  block 

kmaxbl(n)  =  grid  points  in  k  direction  in  nth  block 

pp  =  gauge  pressure  (p-pref) 

u  =  X  component  of  velocity 

V  =  y  component  of  velocity 

w  =  z  component  of  velocity 

t  =  static  temperature 

tke  =  turbulent  kinetic  energy 

tep  =  turbulent  dissipation  rate 

(ys(mi),  nn=l,  ns)  =  mass  fraction  for  each  species 
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For  iuvwt=2:  uvwt.dat  must  be  a  free  formatted  file,  only  the  inlet  face  data  are 
defined  in  this  file  (say  i=l). 

This  option  is  designed  for  the  users  who  want  to  specify  an  inlet  swirling  pro¬ 
file  that  can  be  processed  and  fitted  into  the  grid  points  at  the  inlet  face.  Right 
now,  a  flat  inlet  face  in  the  west  boundary  (i=l  and  x=constant)  is  assumed  to 
apply  this  option.  Users  have  only  to  provide  a  radial  profile  of  the  swirl  cup, 
and  an  axi  symmetric  condition  is  assumed  for  the  swirling  flow  at  the  inlet.  If  a 
user’s  geometry  (and  grid)  does  not  fit  into  the  above  swirling  condition,  the 
user  can  always  use  the  iuvwt=l  option  to  define  the  variables  for  the  entire 
domain  including,  of  course,  the  inlet  swirling  profiles. 


do  n—l,nblk 
read(iunit,*)  irmax 
readi iunit,  rmax,yOyZO 
read(iunit,*)  (yscup(nx),nx=l,ns) 
do  ii-1, irmax 

read(iunit*)  rcup(ii),ucup(ii),vcup(ii),wcup(ii),  tcup(ii),tkcup(ii) 
end  do 
end  do 


where 

nblk  =  number  of  blocks 

irmax  =  number  of  radial  points  used  to  define  a  swirl  profile 
rmax  =  radius  of  the  swirler  cup  (m). 
y0  =  y  coordinate  in  the  center  of  the  swirler  cup  (m) 
zO  =  z  coordinate  in  the  center  of  the  swirler  cup  (m) 

(yscup(iix)^=l^)  =  Mass  fractions  of  the  species 

^gup(ii)  =  r/rmax,  where  r  is  the  radial  distance  from  the  center  of  the  swirler 

ucup(ii)  =  u-velocity  (m/s) 
vcup(ii)  =  v-velocity  (m/s) 
wcup(ii)  =  w-velocity  (m/s) 
tcup(ii)  =  temperature  (K) 
tkcup(ii)  =  turbulence  intensity 
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4.2.5  Property  Data:  prop.dat 

Property  data,  chemistry  kinetics  and  liquid  fuel  property  are  read  by  subrou¬ 
tine  PROPIN  and  the  data  file  name  is  prop.dat. 


c  propane  cbaaiacry 
5 

apaciM  naaa  wn(aa) 

'cShS'  44.09620 

'02*  31.99880 

'co2'  44.00980 

'h2o'  18.01528 

'ii2'  28.01348 

apeeia*  nama  CpO(os) 

*c3bB'  -5.ei6733S4e-01 

'o2'  2.e3608646E*00 

'C02'  2.91174042E«00 

»b20'  3.726394031*00 

'n2'  3.14666076E*00 

apaelaa  aaoa  Bu0(&a] 

'c3he'  -3.54277999E-02 

'o2'  6.5394S961E-06 

'C02'  -1.58795067E-06 

'h2o'  -5.03891166E-06 

'a2'  6.19327062E-06 

apaciaa  Dane  CdOlaa} 

'e3b8'  -3.54277999E-02 

'o2'  3.31925134E-03 

’co2'  -1.I5515278E-02 

'h2o*  -1.087S1854E-02 

♦n2*  4.45461595E-03 

apaciaa  nano  aiGBaCna) 

'c3b8'  5.061 

•o2'  3-467 

'co2'  3.941 

'h20'  2.641 

'a2*  3.798 

apaciaa  naae  aUlcCna] 

'e3h8'  254.0 

'o2'  106.7 

•co2'  195.2 

'h2o*  809.1 

'n2'  71.4 

c  alnpla  atap  propaaa  chemlatry 


cpi(na) 

3.788764998-02 
2. 548823 60E-03 
6.79806423E-03 
5.062 60913E-04 
8.082974328-04 
Eul (na I 

1.70883581E-04 

5.33912193E'08 

5.57379337E-08 

4.76SS0638E-0e 

4.47125S44E-08 

CdKaa) 

1.70883581E-04 

B.27054720E-05 

9.84856560E-0S 

7.44375760E-05 

7.35313003E-05 


Cp2  (Of) 

-2.223S9422E-05 
-1.69585753E-06 
-4.39571263E-06 
1.31410709E-06 
1.12738268B-07 
Eu2 (na) 

-1.75S46ai4E-08 

-1.33802792E-11 

-1.60280202E-11 

-5.41357490E-12 

-1.16368441E-11 

Cd2<na) 

-1.7S546814E-0B 

-1.5206945SE-Oa 

-2.1S630454E-08 

4.30071180E-08 

-1.18328259E-08 


Cp3taa) 

6.63236633E-09 

5.65797716E-10 

1.37312098E-09 

-6.76347697E-10 

-1.74168650E-10 

£u3(na) 

5.580440788-13 

2.64S71434E-15 

2.9S798155E-15 

1.43424321E-16 

2.38277601E-15 

CdJCna) 

S.5S044076E-13 

2.738:0631E-12 

3.5530214BE-12 

-9.10953801E-12 

1.62392e41E-12 


Cp4(na) 

-7.8833e24SE-13 

-7.0e06S610E-14 

-1.6S574806E-13 

9.806783298-14 

3.20583B31E-14 

EU4<ne> 

S.69022174E-17 

-2.07706250E-19 

-2.11109S68E-19 

3.08036420E-20 

-1.91204822E-19 

Cd4(na) 

5. 6902217 4E-17 
-1.93110391E-16 
-2.34533614E-16 
6.56923107E-16 
-9.40907174E-n 


Bl(na) 

-1.38295270E*04 

-9.428713S8C*02 

-4.846919006*04 

-3.022385668*04 

-9.68316219E*02 


B2(&a) 

2.54S79774E*01 

7.83483S21E*00 

7.260358378*00 

1.27805551E*00 

4.SS961323E*00 


c  1  c3hB  ♦  5o2  ->  3eo2  - 
e  nr  ar<nr] 


»  4b2o 
aar(nr) 


c  nx 


8.4a09  -15098.1 

brtnrl  abr(nr) 

0.0  0.0 

x)caqp(nr.BS)  aa-> 

1.0  5.0  0.0  0.0  0.0 


Thermal  Property  Information 


Finite  Rate  Chemistry  Information 


c  ar  xkatpptnr.nal 

1  0.0  0.0  3.0  4.0  0.0 

e  nr  thirdOtar)  chlrd<nr.nx]  ija-> 

1  0.0  0.1  1.6  0.0  0.0  0.0 

1  0.0  0.0  0.0  0.0  0.0  0.0 

c  liquid  pbaaa  oC  propana  (eonditiona  ae  300 

42.64a5  369.64  230.93 

426. 2a3  2.231a3  0.0185 

-6.72219  1.33236  -2.13868 


85.31 

582.0 

-1.38551 


Liquid  Phase  Information 


HGURE6.  Typical />rop.diar  File 


TABLE  5.  Detailed  Information  about  prop.dat  File 


Input  Variable 

Comments 

heading: 

title  of  this  file  (up  to  80  characters) 

hsf: 

number  of  species 

heading: 

title  of  the  property  (up  to  80  characters) 

specnm(nx),wm(nx),nx=l^: 

specnm(nx)  =  species  name  (up  to  80  characters). 
wm(nx)  =  molecular  weight  for  each  species. 

heading: 

title  of  the  property  (up  to  80  characters) 
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TABLE  5.  Detailed  Information  about  prop.dat  File 


Input  Variable 

Comments 

speciim(iix),q)0(iix),cpl(nx),q) 

2(nx),cp3(nx),cp4(nx), 

q)5(nx),cp6(iix) 

coefficients  for  4th  order  specific  heat  polynomial  curve  as  function  ot 
temperature 

The  Cp  curve  is  generated  by  a  forth  order  polynomial  defined  in 

McBride^^.  The  specific  heat,  enthalpy  and  entropy  are  defined  in  Equation 
(30)  and  Equation  (49)  where  cp5(nx)  is  a  constant  whidi  combines  the 
heat  of  formation  of  each  species  at  298.15  K  and  an  integration  constant 
for  the  enthalpy  and  cp6(iix)  is  an  integration  constant  for  the  entropy. 

specnm(nx)=species  name  (up  to  80  characters). 

heading: 

title  of  the  property  (up  to  80  characters) 

specnm(nx),euO(nx),eul(iix),eu 

2(nx),eu3(iix),eu4(nx) 

coefficients  for  4th  order  dynamic  viscosity  polynomial  curve  function  of 
temperature 

beading 

title  of  the  property  (up  to  80  characters) 

specnm(nx),cdO(nx),cdl(nx),cd 

2W)»cd3(nx)>cd^*“) 

coefficients  for  4th  order  thermal  conductivity  polynomial  curve  function 
of  temperature,  Equation  (41) 

heading 

title  of  the  property  (up  to  80  characters) 

specnm(nx)^gm  (nx) 

species  name,  effective  molecular  collision  area.  Equation  (45) 

heading 

title  of  the  property  (up  to  80  characters) 

specnni(nx),edk(nx)  i 

species  name,  effective  tenq)erature(normali2ed  depth  of  potential  energy 
for  a  species).  Equation  (45) 

♦♦Kinetic  Information^^^ 

See  Section  4.4 

heading 

title  of  the  property  (up  to  80  characters) 

pcri,tcri,tboil,tmelt 

Liquid  phase  critical  property  information 
pcri-liquid  phase  critical  pressxire  (N/m^) 
ten  -  liquid  phase  critical  temperature  (K) 
tboil>  liquid  phase  boiling  point  temperature  (K) 
tmelt-liquid  phase  melting  point  temperature  (K) 

hfg,cppar4*kpar^hoipar 

hfg-  latent  heat  of  vj^rization  (J/kg) 
eppar-  liquid  phase  specific  heat  (J/kg-K) 
rkpar-liquid  phase  conductivity  (J/m-s-K) 
riiolpar-liquid  phase  density  (kg/m^) 

vpl,vp2,vp3,vp4 

Critical  pressure  curve  fit  coefficients.  Equation  (1 8) 

4.2.6  User  Defined  Spray  Data;  spray.dat 

When  nsdata=l  and  ispray=l  are  specified,  the  code  expects  users  to  supply 
the  spray  information  from  the  file,  spray.dat.  The  format  is  described  as  fol¬ 
lows. 
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read(iunit,  *)  npgrp 
do  np-1, npgrp 
read(iunit*) 

xpari)( np,  l)pq)aTO( np,2),xpar0(np,3), 
velparO(np,l),velparO(np,2),velparO(np,3), 
tsprO( np ),  radpaWi np),numpar(np) 
end  do 


where 

npgrp  =  number  of  group  of  particles  for  spray 
xpar0(np4)  =  x  coordinate  of  the  initial  spray  location  (m) 
xparO(np^)  =  y  coordinate  of  the  initial  spray  location  (m) 
xparO(np^)  =  z  coordinate  of  the  initial  spray  location  (in) 
velparb(np,l)  =  x  component  of  the  spray  injection  velocity  (m/s) 
yelpar0(np,2)  =  y  component  of  the  spray  injection  velocity  (m/s) 
velparO(np^)  =  z  component  of  the  spray  injection  velocity  (m/s) 
tsprO(np)  =  liquid  temperature  (K)  for  each  group  of  spray 
radparO(np)  =  radius  of  each  drop  size  (m) 
numpar(np)  =  number  of  droplets  per  spray  group  (#/s) 


4.3  Output  Data  Files 

The  ALLSPD-3D  code  generates  up  to  eleven  output  data  files.  The  input  data 
is  echoed  out  to  out.dat.  The  files,  allspd.PFV,  allspd.PFSl,  andallspd.PFS2 
contain  solutions  in  PLOT3D  format.  If  spray  was  used  in  the  calculation,  the 
spray  data  files  loss.dat,  allspdSPY  are  also  generated.  The  restart  information 
is  located  in  the  save.dat  file.  The  convergence  histories  are  written  out  to  the 
res.dat  and  reske.dat  files.  In  addition,  the  turbulence  model  restart  information 
is  written  to  theface.dat  file  and  the  pre-processed  geometry  information  is 
written  to  the  ynbLdat  file. 


43.1  Plot3d  Data  FUes 

There  are  three  plot3d  files.  One  is  the  vector  file,  allspd.PFV,  and  the  other  two 
are  scalar  files,  allspd.PFSl  and  allspd.PFSl.  The  contents  for  each  file 
depends  on  the  user  inputs  specified  in  allspd.dat 


63 


read(iunit,*)nblk 

read(  iunit,  *)(  imax(  n)jTnax(  n),kniax(  n),n=l,nblk) 
do  n-hnblk 

read(  iunit,  '^)finO, angle,  reno,  time 
read(iunit,*) 

(((ql(iJMi=Urruxx(n))J=ljmaMn)),k=l,knuiMn)), 

(((q2(iJMi=l.irnax(n))J=]Jmax(n)),k=l,kmax(n)), 

(((q3(iJMi=l.irnax(n))J=lJnuix(n)),k=l,kniax(n^ 

(((q4(iJMi=l,imax(n))j=lJmax(n)),k=l,kmax(n)). 

(((q5(ijMi-lirf><vc(n))J=lJm(vc(n))Ml,knui^n)) 

end^ 


where 

fmO  =  the  free  stream  Mach  number 
angle  =  flow  angle 
reno  =  flow  Reynolds  number 
time  =  time 

ql  to  q5  =  are  the  flow  variables  defined  in  the  allspd.dat 
4.3^  Spray  Data  Files 

All  spray  files  are  generated  and  read  within  the  subroutine  SPRAY.  The 
spray  results  are  written  into  two  files.  The  file  countdat  contains  a  brief  sum¬ 
mary  of  the  initial  distribution  of  the  spray  particles  groups.  The  file  loss.dat 
contains  a  summary  of  the  spray  group  time  integration  history.  This  file  also 
contains  the  number  of  time  integrations  required  for  a  fuel  droplet  to  vaporize 
or  to  hit  a  particular  boundary.  The  ALLSPD-3D  spray  model  assumes  that 
complete  vaporization  occurs  when  a  drop  hits  a  wall.  It  is  also  assumed  that, 
when  a  droplet  passes  through  the  exit  boundary,  it  is  no  longer  part  of  the  cal- 
culation(i.e.  particle  mass  is  lost).  The  current  code  does  not  continue  to  track 
the  particle  trajectory  once  it  passes  through  the  interface  between  blocks  (if 
more  than  one  block  is  used  in  the  calculation).  A  true  multi-block  spray  code 
will  be  released  in  the  future.  Fortunately,  in  many  practical  spray  applications, 
the  liquid  particles  evaporate  rather  quickly  before  hitting  any  interface  bound¬ 
ary.  The  maximum  time  integration  limit  written  out  to  the  loss.dat  file  should 
be  used  to  set  the  integration  limit  in  the  parameter  statement,  allspdh,  when 
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compiling.  Users  can  also  try  to  increase  npmax  to  a  large  number.  The  spray 
data  are  written  into  allspd.SPY  and  its  format  is  as  follows. 


write(iumt,  npgrp 

do  np=l,npgrp 

write( iunit,  *)  np.npmarkl np) 

write(iunit,*) 

xplot(  npidx,  1  ),xplot( npidx,2),xplot( npidx,  3), 
rplot(npidx),tplot(npidx),dtplot(npidx),npidx=l,npmark(np) 
end  do 


where 

npgrp  =  number  of  group  of  particles  for  spray 

npmark(np)  =  the  last  number  of  time  integration  for  each  group  of  particle 

xpIot(npidx,l)  =  x-coordinate  of  spray  trajectory  (m) 

xplot(npidx^)  =  y-coordinate  of  spray  trajectory  (m) 

xplot(npidx^)  =  z-coordinate  of  spray  trajectory  (m) 

rplot(npidx)  =  radius  of  the  particle  adong  its  trajectory 

tplot(npidx)  =  droplet  surface  temperature  along  its  trajectory  (K) 

dtplot(npidx)  =  integration  time  step  along  its  trajectory  (s) 

43.3  Restart  File 

The  primary  restart  file  is  called  save.dat.  This  file,  and  the  filcsface.dat  and 
ynbl.dat,  are  the  only  files  required  to  restart  a  previous  ALLSPD-3D  run.  They 
contain  all  necessary  data  to  restart  the  program.  It  should  be  noted  that  the 
save.dat  files  generated  by  serial  mode  and  parallel  mode  execution  are  not  com¬ 
patible.  Therefore,  users  should  not  try  to  restart  a  run  in  serial  mode  with  the 
restart  file  that  is  generated  by  a  parallel  version,  and  vice  versa.  The  file 
save.dat  is  generated  in  subroutine  SAVE  and  is  read  by  subroutine  RESTART. 
The  file/flce.ddt  is  generated  in  subroutine  SAVE  and  is  read  in  by  subroutine 
INTT.  The  file  ynbldat  is  generated  in  subroutine  PREGRID  and  is  read  in  by 
subroutine  MESH. 


43.4  Residual  File 

The  code  writes  the  convergence  history  to  both  the  computer  screen  and  files 
named  res.dat  and  reske.dat.  The  residual  shown  on  the  computer  screen  is  the 
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L2  norm  of  the  total  residual  for  flow  variables  and  the  residuals  shown  in  the 
res.dat  contain  residual  for  different  flow  variables.  The  write  statement  is 
write(imit*)  itertot,  rescon,  res_p,  res_u,  res_v,  res_w,  resji  where 

itertot:  iteration  number 

rescon:  L2  norm  of  residual  for  flow  variables 

ies_p:  residual  for  pressure 

res_u:  residual  for  u  velocity 

res_v:  residual  for  v  velocity 

ies_w:  residual  for  w  velocity 

res_h:  residual  for  energy 

The  reskeJat  contains  residual  for  turbulent  variables,  i.e.,  K  and  8. 


4.4  Typical  Chemistry  Input  File  Structure 

The  most  important  rule  to  chemistry  inputs  is  that  N2  must  be  the  last  spe¬ 
cies  in  the  input  file.  For  turbulent  chemistry  model  calculations,  the  first  spe¬ 
cies  and  second  species  must  be  fuel  and  air,  respectively. 
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4.4.1  One-Step  Global  Chemistry  Mode 


TABLE  6.  Typical  Single-step  Chemistry  Input  File 


Input  file  content 

Conunents 

c  methanol  chemistry 

Title 

1 

Number  of  reactions 

c  nr  ar(nr)  ear(nr)  tfrr 

Heading  Title  for  Forward  reaction  rates 

1  4.6ell  -15098.1  0.0 

Reaction  number.  Pre-exponent,  Activation  Ten^ieraturelX], 
Temperature  Correction  Factor 

c  nr  br(nr)  ebr(nr)  tbrr 

Heading  Title  for  Backward  reactions 

c  1  ch3oh  +  1.5o2  ->  co2  +  2h2o 

Title,  helpful  to  write  the  reaction  mechanism  (blank  lines  arc 
acceptable) 

10.0  0.0 

Reaction  number,  Prc-oqponent,  Activation  Ten^rature[K], 
Temperature  Correction  Factor  (note  this  is  forward  only  reac¬ 
tion  and  GUI  will  not  goierate  this  line!) 

c  nr  xkstp(nr4is)  ns-> 

Heading,  Forward  Stoichiometric  Coefficients 

1  1  1.5  00  0 

Reaction  number.  Forward  Stoichiometric  Coefficients  (There 
should  be  nspe  of  them) 

c  nr  xkstpp(nr^)  ns-> 

Heading,  Backward  Stoichiometric  Coefficients 

100120 

Reaction  number.  Backward  Stoichiometric  Coefficients 
(There  should  be  nspe  of  them) 

c  nr  thirdO(nr)  third(nr^)  ns-> 

Heading  Title  for  Correction  Mechanism  (Note,  this  input  is  to 
used  to  change  the  power  of  the  reaction  constant) 

1-1.01.01.00.00.00.0 

Reaction  number.  Reaction  type.  Power  of  the  reaction  mech¬ 
anism  in  forward  direction  (for  the  nspe  species) 

1  -1.00.0  0.00.00.00.0 

Reaction  number.  Reaction  type.  Power  of  the  reaction  mech¬ 
anism  in  Reverse  direction  (for  the  nspe  species) 

Note 

ThirdCKl)  =  -1 , Global  type  kinetics 

ThirdCKl)  =  0,  Modified  kinetics  with  fractional  power 
Third0(l)  >0  detailed  kinetics  with  3-body  corrections 

Kf  =  A  exp(Ea/RT)  [chSohJ*  (02]’ 

Kb  =  0.0 

Mechanism  of  Westbrook  and  Dryer^^ 
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4.42  Simple  Two-Step  Global  Chemistry  Model 


TABLE  7.  Typical  Two-step  Chemistry  Input  File 


Input  file  content 

Comments 

cTwo-step  propane  chemestry 

Title 

2 

Number  of  reactions 

cl  c3h8  +  7/2o2->3co  +  4h2o 

c2  co  +  .5o2<->co2 

Title,  helpful  to  write  the  reaction  mechanism  (blank  lines  are  acceptable) 

c  nr  ar(nr)  ear(nr)  tfrr 

Heading  Title  for  Forward  reaction  rates 

1  5.62e09  -15098.1 0.0 

2  2.24el2  -20140. 0.0 

Reaction  number,  Pre-exponent,  Activation  Ten^ieraturelK],  Temperature  Cor¬ 
rection 

c  nr  br(nr)  ebr(nr)  tbrr 

Heading  Title  for  Backward  reactions 

1  O.OeOOO.0  0.0 

2  5.00e08  -20140  0.0. 

Reaction  number.  Pre-exponent,  Activation  Ten^rature[K],  Temperature  Cor¬ 
rection  factor  (reaction  one  is  a  global  reaction  and  therefore  does  not  have  any 
reverse  reactions) 

c  nr  xkstp(nr4is)  ns-> 

Heading  for  Forward  Stoichiometric  Coefficients 

1  1  3.5  0  0  0  0 

2  0  0.5  1.0  000 

Reaction  number.  Forward  Stoichiometric  Coefficients  (There  should  be  nspe 
of  them) 

c  nr  xkstpp(nr4is)  ns-> 

Heading  for  Backward  Stoichiometric  Coefficients 

1  0  0  3  04  0 

2  0  0  0  1.00  0 

Reaction  number,  Backward  Stoichiometric  Coefficients  (There  should  be 
nspe  of  them) 

c  nr  third0(nr)  third(nr^s)  ns-> 

Heading  Title  for  Correction  Mechanism  (This  input  is  used  to  change  the 
power  of  the  reaction  constant) 

TABLE  7.  Typical  Two-step  Chemistry  Input  File 


Input  file  content 

Comments 

1  0.0  0.1  1.65  0.0  0.0  0.00 

10.0  0.0  0.0  0.0  0.0  0.00 

20.0  0.0  0.25  1.0  0.0  0.5  0 

20.0  0.0  0.0  0.0  1.0  0.00 

Reaction  number,  Reaction  type.  Power  of  the  reaction  mechanism  in  forward 
and  backward  direction  (Note  that  this  medianism  required  modified  power  for 
both  reactions) 

NOTE: 

Third0(l)  =  -1  ,Global  type  kinetics 

ThirdCKl)  =  0,  Modified  kinetics  with  fractional  power 

Third0(l)  >0  detailed  kinetics  with  3-body  corrections 

c3h8  +  7/2  o2  ->  3co  +  4h2o 

Kf  =  A  exp(Ea/RT)  [cshg]®  !  [Oj'  ® 

Kb  =  0.0 

CO  +  .5o2  <->  co2 

Kf  =  A  exp(Eafll.T)  [CO]*  * 

Kb  =  Aexp(Eb/RT)1®[C02]' 

Mechanism  of  Westbrook  and  Dryer^^ 

4.43  Typical  Detailed  Chemistry  Input  file 


TABLE  8.  Detailed  Chemistry  Input  File 


Input  file  content 

Comments 

c  Large  complex  chemistry 

Title 

18 

Number  of  reactions 

c1H2  +  0H<->H20  +  H 
c2  02  +  H<->0H  +  0 
c3H2  +  0<->0H  +  H 
c4  O2  +  H  <->  HO2 
c5H2  +  M<->H  +  H 
c6  02  +  M<->0  +  0 
c7H  +  H02<->H2  +  0 
c8H20  +  M<->OH  +  H 
c9H20  +  0<->20H 
clO  0  +  HO2  <->  C)2  +  OH 
c11H  +  H02<->20H 
cl2  OH  +  HO2  <  ->  0  +  H2O 
cl3  H2  +  H2O2  <->  H2O  +  HO2 
cl4  OH  +  H2O2  <->  H2O  +  H2O2 
cl5  2HO2  +  M  <->  02  +  H2O2 
cl6  H  +  H2O2  <->  OH  +  H2O 
c17H202+M<->20H 

cl8H  +  0<->0H+M 

Title,  helpful  to  write  the  reaction  medianism  (blank  lines  are 
acceptable) 

c  nr  ar(nr)  ear(nr)  tfrr 

Heading  Title  for  Forward  reaction  rates 

TABLE  8.  Detailed  Chemistry  Input  File 


Input  file  content 

Comments 

1  4.74E+10 -3068.95  0.0 

Reaction  number,  Pre-exponent,  Activation  Temperature[K], 

2  1.85E+11  -8253.65  0.0 

Temperature  Correction 

3  4.2E+n  -6919.98  0.0 

4  1.35E+9  503.27  0.0 

5  2.2E+11  -48314.040.0 

6  1.8E+15  -59396.074  1.0 

7  1.3E+11  -1509.8140.0 

8  1.3E+12  -52913.940.0 

9  6.8E+10  -9242.577  0.0 

10  5.0E+10  -503.27  0.0 

11  2.04E+11  -538.5  0.0 

12  8.0E+9  -1499.75  0.0 

13  7.91E+10 -12581.78  0.0 

14  6.1E+9  -719.678  0.0 

15  1.8E+9  0.00.0 

16  7.8E+8  0.00.0 

17  1.44E+14  -22903.875  0.0 

18  7.1E+12  0.00.0 

c  nr  br(nr)  ebr(nr)  tbrr 

Heading  Backward  reactions 

1  2.03138E+12  -10985.23  0.274934 

Reaction  number.  Pre-exponent,  Activation  Temperature[K], 

2  4.523I5E+08  412.94  0.410439 

Temperature  Correction 

3  2.2825E+11  -5998.66  0.0246801 

4  5.1015E+13  -25765.13  0.423663 

5  4.3848E+7  4166.86  0.0320618 

6  1.51713E+9  840.326  0.525515 

7  1.275E+10  -27672.10.4273321 

8  5.73655E+6  7491.96  0.313192 

9  8.61 3E+8  -404.7870.250384 

10  2.6504E+9  -25743.37  0.403293 

11  2.768E+7  -17119.5  0.8083150 

12  4.05365E+10  -35609.15  0.130335 

13  5.07334E+12  -3160.89  0.411882 

14  3.761E+09  -18043.58  0.14644 

15  1.0845E+10  -16734.10.0205195 

16  1.04212E+5  -33982.3  0.89952 

17  8.487E44  3519.1250.0 

18  2.1057E+16  -51573.5  1.0666587 

c  nr  xkstp(nr,ns)  ns-> 

Heading 
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TABLE  8.  Detailed  Chemistry  Input  File 


Input  file  content 

Comments 

1  1.0  0.0  1.0  0.0  0.0  0.0  0.0  0.0  0.0 

2  0.0  1.0  0.0  0.0  1.0  0.0  0.0  0.0  0.0 

3  1.0  0.0  0.0  0.0  0.0  1.0  0.0  0.0  0.0 

4  0.0  1.0  0.0  0.0  1.0  0.0  0.0  0.0  0.0 

5  1.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 

6  0.0  1.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 

7  0.0  0.0  0.0  0.0  1.0  0.0  1.0  0.0  0.0 

8  0.0  0.0  0.0  1.0  0.0  0.0  0.0  0.0  0.0 

9  0.0  0.0  0.0  1.0  0.0  1.0  0.0  0.0  0.0 

10  0.0  0.0  0.0  0.0  0.0  1.0  1.0  0.0  0.0 

11  0.0  0.0  0.0  0.0  1.0  0.0  1.0  0.0  0.0 

12  0.0  0.0  1.0  0.0  0.0  0.0  1.0  0.0  0.0 

13  1.0  0.0  0.0  0.0  0.0  0.0  1.0  0.0  0.0 

14  0.0  0.0  1.0  0.0  0.0  0.0  0.0  1.0  0.0 

15  0.0  0.0  0.0  0.0  0.0  0.0  2.0  0.0  0.0 

16  0.0  0.0  0.0  0.0  1.0  0.0  0.0  1.0  0.0 

17  0.0  0.0  0.0  0.0  0.0  0.0  0.0  1.0  0.0 

18  0.0  0.0  0.0  0.0  1.0  1.0  0.0  0.0  0.0 

Reaction  number.  Forward  Stoichiometric  Coefficients  (There 
should  be  nspe  of  them) 

c  nr  xkstpp(nr4is)  ns-> 

Heading 

1  0.0  0.0  0.0  1.0  1.0  0.0  0.0  0.0  0.0 

2  0.0  0.0  1.0  0.0  0.0  1.0  0.0  0.0  0.0 

3  0.0  0.0  1.0  0.0  1.0  0.0  0.0  0.0  0.0 

4  0.0  0.0  0.0  0.0  0.0  0.0  1.0  0.0  0.0 

5  0.0  0.0  0.0  0.0  2.0  0.0  0.0  0.0  0.0 

6  0.0  0.0  0.0  0.0  0.0  2.0  0.0  0.0  0.0 

7  1.0  1.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 

8  0.0  0.0  1.0  0.0  1.0  0.0  0.0  0.0  0.0 

9  0.0  0.0  2.0  0.0  0.0  0.0  0.0  0.0  0.0 

10  0.0  1.0  1.0  0.0  0.0  0.0  0.0  0.0  0.0 

11  0.0  0.0  2.0  0.0  0.0  0.0  0.0  0.0  0.0 

12  0.0  1.0  0.0  1.0  0.0  0.0  0.0  0.0  0.0 

13  0.0  0.0  0.0  0.0  1.0  0.0  0.0  1.0  0.0 

14  0.0  0.0  0.0  1.0  0.0  0.0  1.0  0.0  0.0 

15  0.0  1.0  0.0  0.0  0.0  0.0  0.0  1.0  0.0 

16  0.0  0.0  1.0  1.0  0.0  0.0  0.0  0.0  0.0 

17  0.0  0.0  2.0  0.0  0.0  0.0  0.0  0.0  0.0 

18  0.0  0.0  1.0  0.0  0.0  0.0  0.0  0.00.0 

Reaction  number.  Backward  Stoichiometric  Coefficients  (There 
should  be  nspe  of  them) 
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TABLE  8.  Detailed  Chemistry  Input  File 


Input  file  content _ Comments _ _ 

c  nr  thirdCKnr)  thirdCnr^ix)  ns->  Heading _ 

I  1.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  Reaction  number.  Reaction  type.  Power  of  the  reaction  mecha- 

21.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  nism  in  forward  direction  (for  the  nspe  species) 

3  1.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 

42.0  2.5  0.3  0,0  20.3  0.0  0.0  0,0  0.0  0.0 

5  1.0  0.0  0,0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 

6  1.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 

7  1.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 

8  2.0  3.2  3.2  15.7  0.0  0.0  0.0  0.0  0.0  0.0 

9  1.0  0,0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 

10  1.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0,0 

II  1.0  0.0  0.0  0.0  0,0  0.0  0.0  0.0  0.0  0.0 

12  1.0  0.0  0.0  0,0  0.0  0.0  0.0  0.0  0.0  0.0 

13  1.0  0.0  0.0  0,0  0.0  0.0  0.0  0.0  0.0  0.0 

14  1.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 

15  1.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 

16  1.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 

17  1.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 
18  2.0  4,0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 

Taken  from  Lee^^ 
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V :  Graphical  User 

Interfacei  GUI) 


This  is  a  preliminary  manual  for  alpha  version  0.1  of  the  ALLSPD-3D 
Graphical  User  Interface  (GUI).  This  GUI  corresponds  to  version  1.0  of  the 
ALLSPD-3D  Code. 


5.1  Installation 

The  ALLSPD-3D  GUI  was  written  in  C  and  requires  X/Motif  Release  1.1. 
It  has  been  built  successfully  on  the  following  platforms:  IBM  RS6000,  SGI, 
Also  beta  ported  to  Solaris  and  HPUX. 

To  install  the  GUI,  first  create  a  separate  directory  and  untar  the  GUI  tar  file. 
This  will  create  a  directory  tree  with  the  following  structure: 

CHANGES  chemistry/  help_files/  species/ 

Makefile  common/  include/  view/ 

Manual/  file/  model/  xmallspd.c 

bin/  grid/  solve/ 

Set  the  XMALLSPD_ROOT  environment  variable  to  point  to  the  directory 
where  the  GUI  will  be  built.  For  instance: 

setenv  XMALLSPD_ROOT  /u/quealy/gui 

Also,  modify  your  path  to  include  $XMALLSPD_ROOT  and 
$XMALLSPD_ROOT/bin.  For  example: 

set  path=(  $path  $XMALLSPD_ROOT  $XMALLSPD_ROOT/bin  ) 

These  can  be  included  in  your  .cshrc  file  so  they  are  issued  every  time  a  new 
shell  is  created. 
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A  top  level  Makefile  will  compile  the  entire  directory  structure  and  create  a 
GUI  executable  named  ‘xmallspd’.  At  this  point,  you  can  type  make ,  and  you 
should  see  the  code  compile.  Sample  output  from  this  stage  is  included  below. 


>make  grid  model  solve  view  ccsnmon*  ;  \ 

for  dir  in  $dir_list  ;  do  \ 

echo  making  GUI  in  directory  Sdir  ;  \ 

{  cd  Sdir  ;  make  all  )  ;  \ 

done 

making  GDI  in  directory  file  „  ^ 

cc  -g  -l/usr/lerc/fsang/gui/gui/include  -c  file.c 

echo  'FILE  module  built. . . ' 

FILE  module  built. .. 

making  GDI  in  directory  grid  ^  j  ^ 

cc  -g  -l/usr/lerc/fsang/gui/gui/xnclude  -c  grid.c 
echo  'GRID  module  built. . . ' 

GRID  module  built . . . 

making  GDI  in  directory  model 

cc  -g  -i/usr/lerc/fsang/gui/gui/include  -c  model. c 

cc  -g  -l/usr/lerc/fsang/gui/gui/include  -c  geom_info.c 

cc  -g  -l/usr/lerc/fsang/gui/gui/include  -c  geonuscale.c 

cc  -g  -l/usr/lerc/fsang/gui/gui/include  -c  geoirL.inplk.c 

cc  -g  -l/usr/lerc/fsang/gui/gui/include  -c  geom_tile.c 

cc  -g  -l/usr/lerc/fsang/gui/gui/include  -c  geom_Pfile.c 

cc  -g  -I/usr/lerc/fsang/gui/gui/include  -c  geom_Cfile.c 

cc  -g  -l/usr/lerc/fsang/gui/gui/include  -c  geoi^util.c 

cc  -g  -I/usr/lerc/fsang/gui/gui/include  -c  probl^.c 

cc  -g  -l/usr/lerc/fsang/gui/gui/include  -c  prop_fluid.c 

cc  -g  -l/usr/lerc/fsang/gui/gui/ include  -c  prop^spray.c 

cc  -g  -l/usr/lerc/fsang/gui/gui/include  *c  prop_chem.c 

cc  -g  -l/usr/lcrc/fsang/gui/gui/include  -c  i)rop_gas_species . c 

cc  -g  -l/usr/lerc/fsang/gui/gui/include  -c  propjreactxon.c 

cc  -g  -l/usr/lerc/fsang/gui/gui/include  -c  prop^util.c 

cc  -g  -i/usr/lerc/£sang/gui/gui/include  -c  prop_ig.c 

cc  -g  -l/usr/lerc/fsang/gui/gui/include  -c  bc.c 

echo  'MODEL  module  built-  -  - ' 

MODEL  module  built — 

making  GDI  in  directory  solve  .  ,  ^  , 

cc  -0  -l/usr/lerc/fsang/gui/gui/xnclude  -c  solve. c 
cc  -g  -I/usr/lerc/fsang/gui/gui/include  -c  control. c 
cc  -g  -l/usr/lerc/fsang/gui/gui/ include  -c  output. c 
cc  -g  -i/usr/lerc/fsang/gui/gui/include  -c  solution. c 
echo  'SOLVE  module  built. .. ' 

SOLVE  module  built... 

making  GDI  in  directory  view  . 

cc  -g  -l/usr/lerc/fsang/gui/gui/include  -c  view.c 
echo  'VIEW  module  built. . . ' 

VIEW  module  built... 

making  GDI  in  directory  common 

cc  -g  -l/usr/lerc/fsang/gui/gui/include  -c  util.c 
echo  'COMMON  module  built- . . ' 

COMMON  module  built. . . 

make  xmallspd  ,  .  ,  .  ti  « 

cc  -g  -l/usr/lerc/fsang/gui/gui/include  7=  ^  ,  . 

cc  -o  xmallspd  xmallspd. o  file/file. o  srid/grid.o  model/model. o 
model/geoaufile.o  model/geon^Pf ile.o 

mode l/geoRuinf 0.0  model/geom_scale.o  model/geo^intblk.o 
mode l/oec«i_uti  1.0  model/problem. o  model/prop_...luid.o 
model /prop_spray.o  model/prop_cbem.o  inodel/prop^as_species-o 
model/prop_reaction.o  model/prop_util.o  model/prop_ig.o 
model/bc.o  solve/solve. o  solve/ con trol .o  11  -IPW 

solve/solution. o  view/view.o  common/util.o  -ixm  -IXt  ixil  ifw 


At  the  time  of  writing  this  manual,  several  assumptions  have  been  made,  espe¬ 
cially  pertaining  to  running  the  parallel  version  of  the  code. 


Assumptions: 

If  PVMis  used,  it  is  installed  and  the  proper  PVM  environment 

variables  have  been  set. _ 

Both  the  GUI  server  and  aU  compute  nodes  share  a  common  file 

system.  _ 

User  ids  are  identical  across  all  platforms. 

If  LSF  is  used,  it  is  available  on  the  GUI  host  as  well  as  the  com¬ 
pute  nodes.  _ _ 

Pre  and  post  processors  exist  on  the  host  where  the  GUI  is  run¬ 
ning,  and  all  executables  are  in  the  user’s  path. _ 

If  no  reactions  are  specified  in  a  pwp.dat  file,  it  is  assumed  that 
there  are  no  liquid  properties. 

The  GUI  sets  default  values  for  most  variables  used  by  the  ALLSPD-3D  code. 
These  default  values  are  set  in  two  locations  in  the  GUI  code: 

model/modelc:  MODELJnvokeModelSimDefaults  () 
solve/solve.c:  SOlVEJnvokeSolveSiTnDefaults  () 

The  user  may  edit  these  routines  and  adjust  the  default  values  if  desired,  and 
then  recompile  the  GUI  to  generate  a  new  ‘xmallspd’  executable. 


5.2  General  Overview 

The  ALLSPD-3D  GUI  allows  the  user  to  create  the  following  ALLSPD-3D 
input  files: 


alispd.dat 

(Input  file) 

prop.dat 

(Property  file) 

giidlink.dat 

(Grid  linkage  file) 
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The  GUI  allows  a  user  to  start  either  the  serial  or  parallel  code  on  local  or 
remote  hosts,  using  several  startup  options.  This  is  subject  to  the  assumptions 
listed  in  the  installation  section. 

The  ALLSPD-3D  GUI  also  gives  a  user  access  to  a  number  of  pre-  and  post¬ 
processing  tools,  again  subject  to  the  assumptions  in  the  installation  section. 

The  following  ‘Scorecard’  is  included  for  people  who  wish  to  work  with  the 
GUI  source  code  itself.  This  is  a  cross  reference  between  variable  names  used 
by  the  ALLSPD-3D  code,  and  their  corresponding  GUI  variable  names. 


All*p<a3<S 


GUI 


nplot 

nsave 

iturb 


idiabc 

icomb 

iartv 

via4 

vis2 

c£l(l> 

cfll2) 

c£l(3> 

pr«c£ac(l) 

precfac(2] 

pr«ef«c(3) 

nawitch 

facini 

rUcacall 

nacallO 

iuwt 

vin 

win 

tref 

pzaf 

rlref 

iabm 

iaqrate 

otligc 

iit2igc 

noigt 

ibteigt 

iligt 

i2igt 

jiigt 

j2igt 

Jtligt 

k2igt 

toff 

apar)cuaar 

nsdata 

iapray 

izospr 

i spray 1 

ispray2 

jsprayl 

jspray2 

icsprayl 

]cspray2 

nanglas 

anglel 

angla2 

sizel 

siza2 

ntspr 

ncalspr 

iinpspr 

tispr 

f lowispr 

vispri 

inplot 

iblanle 


_solve . control .use_raBtart 

Zaolva. control . iteration. iterations 

“solve . control . iteration . iter jplot 

“solve. control. iterat ion. iter_seve 

_;nodel .  prob.  turbulence 

_Aiodel .  prob .  variable_cmu 

_;BOdel .  prob .  epsi  1  on.>c 

_podel  .prob.  eesibustion 

_sol ve. control. nisc.artificial^vi sc 

^solva. control.aisc. coeff 4 

“solve . contro 1 . mi sc . coef  £2 

_aolve . control .ai sc. c£l . eql 

Zsolve. control . aisc. c£l .  oq2 

_solve . control . aisc . cfl .  eq3 

Zsolve. control .aisc. precondition. egl 

“so 1 ve . c ontr o 1 . ai sc . precondi t i on . eg2 

3solva, control. aisc. precondition. eg3 

“solve .  control .  iteration,  iter^ewton 

Isolve. control . iteration. iter^factor 

_solve .  control .  i  teration .  icer_3ee 

_solve . control . i teration . itor_species 

jDodei.prop. fluid. input.src 

^odel.prop.  fluid,  initjvel.^ 

jnodel.prop. fluid. init„vol_y 

_aodel .prop. fluid. init_vel_z 

_podel  .prop,  fluid.  init_toap 

jnodel .prop. fluid. init ^pressure 

jnodel .prop. fluid. ref_length 

_podel .  prob .  e<ldy_br«atojp 

_;aodol .  prop .  chea .  react  ion .  egui  libr  iua_rote 

.Aiodel  .prop,  ignition,  ig^start 

jjodel  .prop .  ignition.  ig_end 

_podol . prop . ignition . nua_srcs ; 

jBodel  .prop,  ignition.  ignition„src->’bloc)?_no 

_podel  .prop,  ignition.  ignition„src->i8tart 

_j»odol  .prop,  ignition.  ignition_src->iend 

_poda I .prop. ignition. ignition_sre->3 start 

jnodeX .prop. ignition. ignition^src->3end 

jaodel  .prop,  ignition.  ignitioa_src->Jcstart 

_podel .prop. ignition. ignition._srC">kand 

_Aiodel.prop.  ignition.  toap_liait 

_aodel .  prop,  ignition.  spar)i;_streagth 

jBodel . prop. spray. input_src 

^odal.prob.  spray 

_podol.prop.  spray. nor  tle_bloe)c 

_jnodel .prop. Spray. location->i start 

_aodel .prop. spray. location->iend 

_aodel .prop. apray . location- >3  start 

_aodel .prop. spray. location- >jend 

_Aiodel .  prop .  spray,  locati  on->Jcstart 

^odel  .prop,  spray,  locati  on->Jcond 

_aodel  .prop .  spray.  nuaL.anglos 

_aodol .prop. spray. angle.start 

_aodel .prop- spray. angle^end 

_podel  .prop  -  spray. particlo_si  redstart 

_podel  .prop.  spray.particlo_sire_ond 

_;do<3o1  .prop .  spray.  sproy_start 

jBodel -prop. spray. spray_£req 

_podol -prop. spray. souree^tera 

jBodel  .prop .  spray,  temp 

jBodel  .prop .  spray,  initial^rop.  input_data.  aass_f  low^rate 
_aodal  .prop,  spray.  initial_prop.  input_data.  vol_aag 
_solve . output .plot3d 
jnodel . geoa. iblank 
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The  main  menu  of  the  ALLSPD-3D  GUI  is  illustrated  in  Figure  7  .  The  remain¬ 
ing  sections  of  this  chapter  highlight  each  of  the  options  from  this  main  menu 
briefly  describes  each  function. 


FIGURE  7.  ALLSPD-3D  GUI  Main  Menu  Window 


5.3  FUe 

The  File  section  of  the  main  menu  allows  you  to  read  and  write  the  three 
main  ALLSPD-3D  files  (input,  property  and  grid  linkage  files),  and  also  to  exit 
the  GUI. 

For  clarification  in  this  manual,  the  files  will  be  referred  to  by  their  default 
names:  allspddat  (input  file),  gridlinLdat  (grid  linkage  file),  and  prop.dat  (property 
file).  Note  that  the  GUI  accepts  and  will  generate  filenames  which  differ  from 
the  default.  The  default  names  are  used  in  this  manual  for  simplicity. 

The  order  of  reading  the  three  input  files  can  be  important.  Since  species  infor¬ 
mation  may  be  included  in  each  of  these  three  files,  it  is  wise  to  read  an  existing 
prop.dat  file  first.  If  prop.dat  is  not  read  first,  all  species  information  at  the  end  of 
the  allspd.dat  file  and  the  gridlintdat  file  is  ignored,  and  a  warning  message  is  dis¬ 
played. 

When  reading  the  gridlink.dat  file,  the  grid.dat  filename  is  also  requested.  The 
GUI  does  not  read  the  entire  grid.dat  file.  It  simply  reads  the  initial  information 
in  the  file  to  verify  that  the  gridlink.dat  file  matches  the  gruLdat  file. 
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5.4  Grid 

The  Grid  section  of  the  main  menu  gives  the  user  access  to  grid  generation 
packages.  Currently  the  GUI  allows  starting  the  CFD-GEOM  package  devel¬ 
oped  by  CFD  Research.  Gridgen  access  has  not  yet  been  implemented. 

The  user  may  set  the  XMALLSPD_CFDGEOM_EXEC  environment  variable 
to  point  to  the  location  of  the  CFD-GEOM  executable.  If  this  environment  vari¬ 
able  is  not  set  it  is  assumed  the  default  name  ‘CFD-GEOM’  is  used,  and  that 
the  executable  is  in  the  user’s  path.  It  is  assumed  that  licenses  to  access  CFD- 
GEOM  are  in  place,  and  that  CFD-GEOM  will  run  on  the  same  host  as  the 
GUI. 


5.5  Model 

The  model  menu  is  illustrated  in  Figure  8  .  Each  of  the  options  from  this 
model  menu  are  briefly  described  below 


5.5.1  Model/File 

The  File  option  from  the  Model  Module  simply  allows  you  to  close  the  Mod¬ 
el  Module  window. 
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5  J.2  Model/Manipulate  Grid 

The  ‘Manipulate  Grid’  option  from  the  Model  Module  window  allows  the 
user  to  import  a  grid  developed  externally,  display  grid  information,  perform 
grid  scaling  if  necessary,  and  define  internal  blockages. 

The  GUI  currently  supports  PLOT3D  and  CFD-GEOM  grid  files.  The  files  can 
be  2D  or  3D  files.  All  2D  files  are  converted  by  the  GUI  to  a  3D  format  usable 
by  ALLSPD-3D.  The  option  to  import  a  Gridgen  grid  file  has  not  been  imple¬ 
mented  yet. 

The  ‘Grid  Information’  option  displays  grid  parameter  information  for  each 
grid  block.  If  a  grid  file  has  not  been  specified,  and  error  message  is  displayed. 

The  ‘Grid  Scaling’  option  allows  the  user  to  scale  the  current  grid  from  either 
inches,  feet  or  centimeters  to  meters  (which  is  what  AIiSPD-3D  expects).  The 
user  is  allowed  the  option  to  select  a  new  file  name.  If  a  new  file  name  is  not 
selected,  the  original  grid  is  overwritten. 

WARNING:  Since  this  is  an  alpha  version  of  the  GUI,  it  is  STRONGLY 
ADVISED  that  you  do  not  overwrite  the  original  grid.  Any  errors  occurring 
during  the  scaling  process  could  destroy  the  original  file.  To  prevent  this  from 
happening,  keep  a  copy  of  the  original  files  in  a  protected  place  with  read  only 
access. 

The  ‘Define  Internal  Blockages’  option  allows  the  user  to  define  any  number  of 
internal  blockages  for  each  block. 

5.5.3  Model/Problem  Description 

The  ‘Problem  Description’  section  of  the  Model  Module  allows  the  user  to 
select  several  problem  options.  The  options  to  select  are: 
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If  the  application  is  Turbulent,  the  user  has  the  additional  option  of  using  a 
variable  C\i  in  the  turbulent  eddy  viscosity  formulation  or  the  CMOTT  e 
boundary  condition  in  place  of  the  default  £  boundary  condition. 

5.5.4  Model/Properties 

The  ‘Properties’  section  of  the  Model  Module  is  the  most  complex  section 
of  the  ALLSPD-3D  GUI.  Under  this  option,  the  user  may  specify  Fluid  proper¬ 
ties,  Spray  properties.  Chemistry  properties  (Thermal  &  Kinetics),  and  Ignition 
properties.  Each  of  these  areas  is  described  in  detail  below. 

5.5.4.1  Model/Properties/Fluid 

Under  the  Fluid  Properties  window,  the  user  is  able  to  specify  various  initial 
fluid  properties,  including: 

Initial  Velocity  (x.y,z) _ 

Initial  Temperature 
Initial  Pressure 

Reference  length _ 

Also,  the  user  is  given  four  options  for  setting  initial  flow  conditions: 

Initial  velocity  plane  is  defined  ONLY  _ 

Uniform  Initial  conditions  _ 

User  Defined  Initial  conditions  (expects  uvwtfile) _ 

Initial  plane  firom  Swirl  Cups  (expects  uvwt  file) 

At  the  time  of  writing  this  manual  the  uvwt.dat  filename  must  be  entered  via  the 
Model/Properties/Fluid  window  if  the  allspd.dat  file  indicates  that  such  a  file  is 
to  be  read.  If  the  ALLSPD-3D  code  expects  the  filename  but  it  has  not  been 
entered  through  the  GUI,  an  eiror  message  is  displayed  when  the  user  attempts 
to  start  the  code. 


5  J.4.2  General  Information  on  Species  Database 

The  GUI  relies  on  a  species  database  to  obtain  property  information  for  the 
pwp.dat  file.  The  species  database  that  is  distributed  with  the  GUI  source  code  is 
located  in  $XMALLSPD_ROOT/species/speciesJist.  There  are  currently  31  species  in 
this  database. 

The  user  is  able  to  customize  this  species  database  by  entering  information  via 
the  GUI  and  storing  the  newly  created  species  information  into  the  database. 
This  allows  the  user  to  add  to  this  database  as  desired. 

WARNING:  Since  this  is  an  alpha  version  of  the  GUI,  it  is  STRONGLY 
ADVISED  that  you  keep  a  backup  copy  of  any  species  database  that  is  built  by 
the  user.  Any  errors  occurring  during  the  writing  of  new  information  to  the 
existing  species  database  could  destroy  the  existing  database.  To  prevent  this 
from  happening,  keep  a  copy  of  a  customized  species  database  in  a  protected 
place  with  read  only  access.  Also,  if  you  are  in  the  process  of  creating  a  species 
database,  save  it  frequently  and  check  for  errors.  And  do  not  forget  to  make 
backups! 

5.5.4,3  Locatii^  the  Master  Species  Database 

The  GUI  follows  a  pre-defined  search  path  when  trying  to  locate  the  Master 
Species  Database.  First  the  environment  variable 
XMALLSPD_SPECIES_LIST  is  examined.  If  this  environment  variable  has 
been  defined,  the  species  database  is  read  from  this  location.  If  not  defined,  the 
search  path  continues  in  the  following  order: 

.xmallspd/species_list 

$HOME/.xmallspd/species_list 

$XMALLSPD_ROOT/species/speciesJist 

The  first  species  database  file  which  is  located  is  then  used  to  read  the  master 
species  list  into  the  database. 

When  updating  the  master  list,  the  file  search  order  is  as  follows: 
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$XMALLSPD^SPECIESJJST 

.xmallspd/speciesjist 

$HOME/.xmallspd/speciesJist 

If  the  file  is  not  located  in  any  of  these  locations,  this  means  that  the  original 
distribution  species  list  has  been  altered.  To  prevent  overwriting  the  original, 
the  .xmallspd  directory  is  created  in  the  current  working  directory,  and  the  spe¬ 
cies  list  is  written  there  (.xmallspd/speciesjist). 

The  GUI  uses  the  chemical  name  for  searching  the  database,  so  this  must 
always  be  unique.  Also,  at  the  time  of  writing  this  manual,  the  GUI  is  ‘case 
sensitive.’  Most  species  in  the  database  have  been  entered  with  uppercase  let¬ 
ters.  Some  however  require  lowercase  letters  (such  as  iC4H10  for  iso-Butane). 

5  J.4.4  Locating  the  Master  ‘Default  Chemistry  Package’  Database 

The  GUI  follows  a  pre-defined  search  path  when  trying  to  locate  the  Master 
‘Default  Chemistry  Package’  Database.  First  the  environment  variable 
XMALLSPD_CHEMISTRY_LIST  is  examined.  If  this  environment  variable 
has  been  defined,  the  DefaultChemList  (and  any  selected  packages)  are  read 
from  this  location.  If  hot  defined,  the  search  path  continues  in  the  following  or¬ 
der; 

jcnudlspd/DefaultChemList 

$HOME/.xmallspd/DefaultChemList 

$XMALLSPD_ROOT/chemistry/D(^aultChemUst 

The  first  DefaultChemList  database  file  which  is  located  is  then  used  to  read 
the  master  ‘Default  Chemistry  Package’  fist  into  the  database. 

When  updating  the  master  ‘Default  Chemistry  Package’  list,  the  file  search 
order  is  as  follows: 

$XMALLSPD_CHEMISTRYJJST.xmallspd/DefaultChemList 

$HOME/.xmallspd/D^aultChemList 
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If  the  file  is  not  located  in  any  of  these  locations,  this  means  that  the  original 
distribution  DefaultChemList  has  been  altered.  To  prevent  overwriting  the  orig¬ 
inal,  the  jcmallspd  directory  is  created  in  the  current  working  directory,  and  the 
DefaultChemList  is  written  there  (.xmallspd/D^aultChemList), 

5^.4.5  Model/Properties/Spray  Initialization 

The  user  cannot  enter  the  Spray  Initialization  window  unless  the  Spray  prob¬ 
lem  option  has  been  set.  Once  this  window  is  activated,  the  user  may  either  enter 
spray  initialization  information  (default),  or  read  this  information  from  a  file 

(spray.dat). 

At  the  time  of  writing  this  manual  the  spray.dat  filename  must  be  entered  via  the 
Model/Properties/Spray  window  if  the  allspd.dat  file  indicates  that  such  a  file  is 
to  be  read.  If  the  ALLSPD-3D  code  expects  the  filename  but  it  has  not  been 
entered  through  the  GUI,  an  error  message  will  be  displayed  when  the  user 
attempts  to  start  the  code. 

If  the  user  opts  to  enter  the  spray  initialization  information,  another  window 
opens  allowing  the  user  to  enter  information  such  as: 

Block  with  nozzle 

Spray  location 

Number  of  angles  and  range 

Particle  size  information 

Spray  source  term  treatment  option 

Iteration  to  start  spray 

Spray  Frequency 

Spray  temperature 

Mass  flow  rate 

Velocity  magnitude  _ 


5^.4.6  Model/Properties/Customize  Liquid  Species 

The  user  is  able  to  either  select  the  Spray  Species  from  a  list,  or  to  customize 
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the  spray  property  information  via  the  GUI.  If  the  user  chooses  to  select  the 
spray  species  from  a  list,  the  list  is  built  from  the  current  species  database,  as 
determined  above. 

If  the  user  opts  to  customize  the  spray  species  information,  this  information  can 
then  be  added  to  the  user’s  species  database  for  future  use.  Properties  to  enter 
for  customizing  include: 


Species  Name 

(for  example,  ‘Oxygen’) 

Chemical  Name 

(for  example,  ‘02’) 

Critical  Pressure 

Critical  Temperature 

Melting  Temperature 

Boiling  Temperature 

Phase  change  enthalpy 

Specific  heat  for  liquid  phase 

Thermal  conductivity  for  liquid  phase 

Density  at  injection  temperature 

Vapor  pressure  coefficients 

5.5.4.7  Model/Properties/Read  Default  Chemistry  Package 

This  section  of  the  GUI  allows  the  user  to  read  a  pre-defined  default  chem¬ 
istry  package.  Several  default  chemistry  packages  have  been  distributed  with  the 
ALLSPD-3D  GUI.  The  following  packages  currently  exist: 
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By  building  the  species  database  and  customizing  the  thermal  properties  and 
kinetics,  the  user  can  continue  to  build  and  add  to  the  default  chemistry  pack¬ 
ages. 

5.5.4.8  ModeiyPropertiesAVrite  Default  Chemistry  Package 

This  section  of  the  GUI  allows  the  user  to  write  a  customized  chemistry 
package  to  the  chemistry  database.  Thermal  and  Kinetics  information  must  be 
specified  before  the  user  is  allowed  to  generate  a  chemistry  package.  Once  the 
new  package  is  written  to  the  database,  the  user  is  able  to  read  the  default  chem¬ 
istry  package  just  like  the  pre-defined  packages  which  were  supplied  with  the 
GUI. 

NOTE:  Currently  if  a  name  is  entered  for  a  new  Default  Chemistry  Package 
which  duplicates  a  pre-existing  package,  the  GUI  posts  an  error  message  and 
expects  the  user  to  select  a  new  name.  More  flexibility  may  be  added  at  a  later 
date. 

5.5.4.9  Model/Properties/Customize  Thermal 

The  ‘Customize  Thermal’  option  allows  the  user  to  either  select  Thermal 
Properties  firom  the  species  database,  or  to  enter  new  species  thermal  properties 
by  hand. 

If  the  user  chooses  to  enter  the  Thermal  property  information,  the  following 
information  is  requested  for  the  species  database: 


Species  name 

(for  example.  Oxygen) 

Chemical  name 

(for  example,  ‘02’) 

Specific  Heat 

Enthalpy  constant 

Entropy  constant 

Viscosity 

Thermal  Conductivity 

Diffusion  parameters 
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There  are  some  restrictions  upon  the  order  in  which  species  should  be  selected. 
The  eddy  breakup  model  assumes  that  the  major  fuel  constituent  is  first  and 
that  oxygen  (02)  is  second.  This  order  is  also  assumed  for  a  global  reaction 
mechanism.  In  general,  if  nitrogen  is  chosen  as  a  component,  it  should  be 
picked  as  the  last  species  in  the  list  so  that  round-off  error  will  be  minimized. 

5.5.4.10  Model/Properties/Customize  Kinetics 

The  user  cannot  enter  the  ‘Customize  Kinetics’  window  unless  the  ‘React¬ 
ing’  problem  option  has  been  set.  Once  this  window  is  activated,  the  user  is  able 
to  enter  reaction  information. 

The  user  may  enter  a  reaction  set  name.  Also,  the  user  must  specify  the  known 
reaction  rates  (forward  reaction  rate  constants  only,  backward  reaction  rate 
constants  only,  or  both)  for  all  reactions. 

The  user  may  then  add  any  number  of  reactions.  For  each  reaction,  the  user  is 
asked  for: 

Reaction  type  (Global,  Westbrook  &  Dryer,  Finite  Rate) _ _ 

Forward  Modified  Arrhenius  rates _ _ 

pre-exponential  _ _ _ - _ 

temperature  exponent 

activation  energy  _ _ _ _ 

Backward  Modified  Arrhenius  rates _ _ 

pre-exponential  _ _ _ _ 

temperature  exponent 

activation  energy _  _ 

And  for  each  species. 

Forward  Stoichiometric  Coefficient _ 

Forward  Mass  Action  Coefficient  (if  appropriate) _ 

Backward  Stoichiometric  Coefficient  _ 
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Backward  Mass  Action  Coefficient  (if  appropriate) 
Third  Body  Efficiency  (if  appropriate) 


5.5.4.11  Ignition 

The  ‘Ignition’  window  allows  the  user  to  enter  various  ignition  properties, 
including: 


Iteration  to  apply  ignition 
Iteration  to  stop  ignition 
Temperature  limit 
Spark  strength 
Ignition  Locations 


5.5.4.12  Initial  Mass  Fractions 

Once  the  user  selects  species  information  (via  ‘Default  Chemistry  Package’ 
or  ‘Customize  Thermal’),  a  window  requesting  mass  fraction  information  is  dis¬ 
played.  The  user  is  asked  to  enter  an  initial  mass  fraction  value  for  each  species 
that  has  been  selected. 

5.5.5  Model/Boundary  Conditions 

The  ‘Boundary  Condition’  window  allows  the  user  to  specify  boundary  con¬ 
ditions  for  each  block  on  each  of  six  faces.  It  also  allows  specifying  Dilution 
Holes  and  Transpiration  Cooling  areas. 

The  boundary  condition  options  include: 

Upstream 
Downstream 
Adiabatic  Wall 
Isothermal  Wall 
Block  Interface 
Symmetry 
Periodic 


89 


I  Singularity  _ 

Far  Field  _ 

If  a  downstream  or  far  field  boundary  condition  is  specified,  the  user  is  asked 
for  an  exit  pressure  value.  If  an  isothermal  wall  is  specified,  the  user  is  asked  to 
enter  a  wall  temperature  value.  If  a  block  interface  is  specified,  the  user  is 
asked  to  enter  the  block  number  of  the  neighboring  block. 

If  a  particular  block  has  one  or  more  dilution  holes,  the  user  can  enter  dUution 
hole  information  by  pressing  the  ‘Dilution  Hole’  pushbutton.  The  user  may 
enter  any  number  of  dilution  holes.  The  user  is  asked  to  enter  the  followmg 
information  for  each  dilution  hole: 


Dilution  hole  location _ 

Mass  flow  rate 

Temperature _ _ 

Turbulent  intensity _ 

Direction  cosine  (x,y,z)  _ 

If  species  have  been  selected,  the  user  is  asked  to  enter  dilution  mass  firactions 
for  each  species. 

If  a  particular  block  has  one  or  more  Transpiration  Cooling  areas,  the  user  can 
enter  transpiration  cooling  information  by  pressing  the  ‘Transpiration  Cooli^’ 
pushbutton.  The  user  may  enter  any  number  of  transpiration  coolmg  areas.  The 
user  is  asked  to  enter  the  foUowing  information  for  each  transpiration  coohng 

area: 

Transpiration  cooling  area  location _ 

Mass  flow  rate _ 

Direction  cosine  (x,y,2) _ 

If  species  have  been  selected,  the  user  is  asked  to  enter  transpiration  mass  frac¬ 
tions  for  each  species. 


5.5.6  Model/Help 

Help  files  are  included  in  the  GUI.  If  a  help  file  does  not  currently  exist, 
pressing  that  ‘Help’  button  will  have  no  affect. 


5.6  Solve 

The  solve  menu  is  illustrated  in  Figure  9 .  Each  of  the  options  from  this  solve 
menu  are  briefly  described  below. 


5.6.1  Solve/File 

The  File  option  from  the  Solve  Module  simply  allows  you  to  close  the  Solve 
Module  window. 

5.6.2  Solve/Control 

The  ‘Control’  menu  option  allows  the  user  to  specify  whether  a  restart  file 
is  to  be  used,  and  to  set  iteration  and  other  miscellaneous  parameters. 

The  ‘Restart’  option  allows  the  user  to  specify  whether  a  restart  file  is  to  be 
used.  If  a  restart  file  is  to  be  used,  a  file  selection  box  is  displayed  so  the  user 
may  enter  the  name  of  the  restart  file.  Note  that  if  you  are  running  the  parallel 
code,  a  restart  file  exists  for  each  process.  These  are  named: 

save.dat-001 

save.dat-002 
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save.daUOOn 

where  n  is  the  number  of  processes  used  in  this  particular  calculation.  When 
entering  the  filename  for  the  restart  file  in  this  situation,  do  not  enter  a  specific 
file  extension.  In  other  words,  enter  save.dat  rather  than  save.dat-OOl. 

Also  note  that  a  restart  file  created  by  the  serial  code  is  NOT  compatible  with 
the  parallel  code,  and  vice  versa. 

At  the  time  of  writing  this  manual  the  restart  filename  must  be  entered  via  the 
Solve/Control  window  if  the  allspd.dat  file  indicates  that  a  restart  file  is  to  be 
used.  If  the  ALLSPD-3D  code  expects  the  filename  but  it  has  not  been  entered 
through  the  GUI,  an  error  message  will  be  displayed  when  the  user  attempts  to 
start  the  code. 

The  ‘Iterations’  window  allows  the  user  to  specify  the  following  parameters: 


Total  number  of  Iterations  _ _ _ _ 

Iteration  frequency  for  writing  a  PLQT3D  file  _ 

Iteration  frequency  for  writing  a  restart  file _ _ _ _ 

Number  of  iterations  between  turbulent  equation  iterations _ _ 

Number  of  iterations  between  solutions  to  the  species  continuity  equations _ 

Number  of  iterations  to  under  relax  the  time  step _ _ _ 

Relaxation  factor  for  time  step  _ 

The  ‘Miscellaneous’  window  allows  the  user  to  specify  the  following  parame¬ 
ters: 


Artificial  \^scosity _ 

CFL  Numbers 

Pre-conditioning  numbers _ 

If  Artificial  Viscosity  is  ‘ON’,  the  user  is  asked  for  2nd  and  4th  order  coeffi¬ 
cients. 
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5.6.3  Solve/Output 

The  ‘Output’  window  allows  the  user  to  specify  PLOT3D  parameters.  The 
user  is  given  a  list  of  variables  which  can  be  written  to  PLOT3D  files.  The  list 
includes: 


Density 

U  Momentum 

V  Momentum 

W  Momentum 

Static  Temperature 

Gage  Pressure 

Static  Pressure 

T\irbulent  CMU 

Static  Enthalpy 

Total  Enthalpy 

Turbulent  Kinetic  Energy 

Dissipation 

Turbulent  Eddy  Viscosity/Laminar  Viscosity 

Normal  Distance 

plus  all  mass  fractions  for  the  selected  species. 


A  maximum  of  15  choices  can  be  selected.  These  are  written  in  the  order 
selected  to  three  PLOT3D  files.  The  first  five  selections  are  written  to  the  first 
PLOT3D  file  (default  =  allspdLPFV).  For  compatibility  with  CFD-VIEW,  this  file 
should  include  any  vector  quantities,  e.g.  momentum,  as  the  2nd,  3rd,  and  4th 
choices.  The  second  five  selections  are  written  to  the  second  PLOT3D  file 
(default  =  cdlspd-PFSl).  And  the  third  set  of  five  are  written  to  the  final  PLOT3D 
file  (default=  allspdPFSZ).  These  files  are  written  in  ASCII  format  for  compati¬ 
bility  across  platforms. 

5.6.4  Solve/Solution 

The  ‘Solution’  option  from  the  Solve  menu  allows  the  user  to  setup  an  ap- 
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plication,  and  start  and  stop  the  application. 

5.6.4.1  Solve/Solution/Setup  Environment 

The  ‘Setup  Environment’  option  allows  the  user  to  specify  whether  the 
ALLSPD-3D  executable  has  been  compiled  for  ‘Serial’  or  ‘Parallel’  mode.  It 
also  allows  the  user  to  specify  the  job  submission  mode  (Interactive  or  Batch), 
and  the  host  selection  mechanism. 

In  interactive  mode,  the  user  can  specify  the  hosts  to  be  used,  or  allow  LSF  to 
select  the  hosts  to  be  used,  subject  to  user  specified  resource  requirements. 
Refer  to  the  LSF  manual  for  details  on  specifying  resource  requirements.  An 
example  resource  string  might  be 

cs  (type  ==  RS6K) 

which  might  mean  that  the  user  wants  an  RS6000  compute  server.  In  batch 
mode,  the  user  can  submit  a  job  to  a  specified  LSF  queue.  Again,  refer  to  the 
LSF  manual  for  details  on  specifying  batch  queues. 

If  the  user  manually  selects  hosts  to  run  the  application,  the  GUI  searches  for  a 
host  list  specifying  existing  valid  hosts.  The  XMALLSPD_HOST_LIST  envi¬ 
ronment  variable  defines  the  location  of  this  host  list.  If  this  environment  vari¬ 
able  has  not  been  set,  the  GUI  follows  a  search  path  through: 

,xmall$pd/host_list 

$HOME/.xmallspd/hostJrist 

If  these  files  cannot  be  found,  the  GUI  creates  a  new  file  in  .xmallspd/hostjist 
using  the  current  host  as  the  initial  entry. 

5.6.42  Solve/Solution/Solve 

The  ‘Solve’  option  allows  the  user  to  specify  whether  the  grid  pre-processor 
should  be  run,  and  actually  submits  the  job  to  the  selected  environment.  Note 
that  at  this  time,  the  ALLSPD-3D  executable  must  be  compiled  to  match  the 


grid  size  and  number  of  species  selected  by  the  GUI.  Eventually  the  GUI  may 
query  the  code  to  verify  that  the  ALLSPD-3D  code  matches  the  GUI  inputs, 
however  this  has  not  been  implemented  at  this  time. 

The  Solve  option  starts  a  ‘run-allspd’  script  which  is  located  in  the 
$XMALLSPD_ROOT/bin  directory.  The  run-allspd  script  uses  ‘alhpdSd’  as  the 
default  name  for  the  ALLSPD-3D  executable.  This  can  be  modified  by  setting 
the  XMALLSPD_EXEC  environment  variable,  'pregrid’  is  the  default  name  for 
the  grid  pre-processor.  Again,  this  can  be  modified  by  setting  the 
XMALLSPD_PREGRID  environment  variable. 

The  GUI  creates  three  ‘temporary’  input  files  when  the  ‘SOLVE’  button  is 
pushed.  These  files  reflect  the  current  values  entered  via  the  GUI.  These  files 
are  named  JOB.allspd3d_input.<pid>,  JOB.allspd3d_prop.<pid> ,  and 
JOB.allspd3dJink.<pid>,  where  <pid>  is  the  UNIX  process  id  of  the  run-allspd 
script.  These  files  are  written  to  the  current  working  directory  and  are  used  as 
input  to  the  ALLSPD-3D  code. 

5.6.43  Solve/Solution/Stop 

At  the  time  of  writing  this  manual,  the  Stop  option  has  not  been  fully  imple¬ 
mented.  Currently  the  user  may  stop  a  serial  or  parallel  job  which  is  started  on 
the  local  host  (i.e.,  the  host  which  is  running  the  GUI),  and  a  parallel  job  which 
is  stented  on  a  remote  host.  Jobs  submitted  to  a  batch  scheduler  should  be 
stopped  by  using  the  batch  scheduler  stop  mechanism.  Serial  jobs  which  are 
started  on  remote  hosts  must  be  killed  manually  using  the  UNIX  ‘kill’  com¬ 
mand. 

5.6.5  Solve/Help 

The  help  files  are  included  in  the  GUI.  If  a  help  file  does  not  currently  exist, 
pressing  that  ‘Help’  button  will  have  no  affect. 
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5.7  View 

The  View  section  of  the  main  menu  gives  the  user  access  to  post-processing 
packages.  Currently  the  GUI  aUows  starting  FAST,  and  the  CFD-VIEW  pack¬ 
age  developed  by  CFD  Research. 

The  user  may  set  the  XMALLSPD_CFDVIEW_EXEC  environment  variable 
to  point  to  the  location  of  the  CFD-VIEW  executable.  If  this  environment  vari¬ 
able  is  not  set  it  is  assumed  that  the  default  name  ‘CFD-VIEW’  is  used,  and 
that  the  executable  is  in  the  user’s  path.  It  is  assumed  that  licenses  to  access 
CFD-VIEW  are  in  place,  and  that  CFD-VIEW  will  run  on  the  same  host  as  the 
GUI. 

The  user  may  set  the  XMALLSPD_FAST_EXEC  environment  variable  to 
point  to  the  location  of  the  FAST  executable.  If  this  environment  variable  is  not 
set  it  is  assumed  that  the  default  name  ‘fast’  is  used,  and  that  the  executable  is 
in  the  user’s  path.  It  is  assumed  that  licenses  to  access  FAST  are  in  place,  and 
that  FAST  will  run  on  the  same  host  as  the  GUI. 


5.8  Help 

Most  of  the  help  files  are  included  in  the  GUI,  however  some  have  not.  If  a 
help  file  does  not  currently  exist,  pressing  that  ‘Help’  button  will  have  no  affect. 


5.9  Known  Bugs 

Several  locations  in  the  GUI  require  a  user  to  hit  an  “Add”  button  before  en¬ 
tering  data.  If  the  user  enters  data  without  hitting  Add,  the  data  appears  to  be  en¬ 
tered,  but  it  is  actuaUy  lost.  This  must  be  modified  to  prevent  a  user  from 
entering  data  if  a  specific  item  has  not  been  added. 
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5.10  Future  Plans 

The  following  items  may  be  modified  in  the  future: 

The  assumptions  regarding  the  GUI  host  and  all  compute  nodes  sharing  a  com¬ 
mon  file  system  should  be  relaxed  to  allow  an  arbitrary  file  system. 

Currently  an  entry  cannot  be  deleted  from  the  species  database  or  the  default 
chemistry  packages.  More  flexibility  to  manipulate  these  databases  should  be 
added. 
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VI :  Code  Structure  and 
Hardware 
Implementations 


6.1  General  Structure 

The  ALLSPD-3D  code  is  not  currently  considered  a  “black  box”  production 
code,  however,  it  is  not  a  pure  “research”  code  either.  During  the  course  of  this 
development,  we  have  tried  very  hard  to  run  test  cases  by  modifying  only  the 
user  interface  files,  not  the  source  code.  The  current  version  of  the  code  can  run 
in  either  serial  or  parallel  mode.  In  other  words,  it  can  run  on  a  Cray-YMP,  a 
single  workstation  (SGI,  IBM,  HP,...),  or  a  cluster  of  workstations.  Also,  the 
Graphical  User  Interface  (GUI)  allows  users  to  build  input  data  files  (grid, 
chemistry,  flow  conditions)  and  initiate  interactive  post-processing.  This  sec¬ 
tion  wiU  discuss  these  features.  The  ALLSPD-3D  code  consists  of  a  main  pro¬ 
gram  and  a  set  of  subroutines.  The  general  structure  illustrated  in  Figure  10 
shows  a  flowchart  encompassing  the  entire  calculation.  Beside  each  box  in  the 
flow  diagram  appears  the  name(s)  of  the  primary  subroutine(s)  responsible  for 
the  associated  task.  Short  conunents  in  each  subroutine  source  listing  also  de¬ 
scribe  its  function. 
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PRIMARY 

SUBROUTINES: 


SUPPORTING 

SUBROUTINES: 


QPLUSl 


FIGURE  10.  ALLSPD-3D  Flow  Chart 
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FIGURE  11.  ALLSPD-3D  Flow  Chart  (cont’d) 
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MULTI-BLOCK  LOOP 


> 


6.2  Working  Units 

The  M-K-S  system  of  units  is  used  for  all  variables  throughout  the  code.  The 
units  for  several  variables  are  listed  as  follows. 

Length;  meter  (m) 

Time:  second  (s) 

Mass:  kilogram  (kg) 

Temperature:  Degrees  Kelvin  (K) 

Energy:  Joules  (J) 

Force;  Newton  (N) 

Velocity:  meter/second  (m/s) 

Pressure:  N/m^ 

Density;  kg/m^ 

Specific  heat:  J/(kg  K) 

Dynamic  viscosity:  kg/(m  s) 

Kinematic  viscosity:  m^/s 

The  dimensional  units  required  for  chemical  kinetic  rate  pre-exponents  Af  and 
Ab  are  in  the  M-K-S  system.  Typically,  the  activation  energy  is  normalized  by 
the  ideal  gas  constant  and  represented  as  an  activation  temperature  in  the  code. 
The  required  dimensional  units  for  these  quantities  are, 


where. 


[A^(*)]or[A^(/:)] 


kmols 


l-rm* 


3nnt-3l 


sec 


(80) 
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N 

N 

or 

1  '’’it 

i  =  1 

i  =  1 

respectively,  and 


(81) 


(82) 


The  ALLSPD-3D  finite  chemistry  modules  have  been  developed  with  a  gener¬ 
alized  chemistry  formulation.  Therefore,  any  finite  rate  and  global  chemistry 
models  that  can  be  fitted  into  the  Arrhenius  form  shown  in  Equation  (25)  can 
be  used  with  the  code.  The  method  of  inserting  various  chemical  reaction 
mechanisms  is  shown  in  the  Input/Output  Files  section  along  with  several 
examples. 


63  Coordinate  System 

The  governing  equations  in  the  code  are  transformed  from  the  original  Car¬ 
tesian  components  to  generalized  coordinates.  Either  right-hand  or  left-hand 
rule  based  coordinates  are  acceptable  for  an  ALLSPD-3D  calculation,  although 
the  right-hand  rule  is  the  standard  convention  and  most  commercial  grid  gener¬ 
ation  tools  will  not  allow  the  use  of  left-handed  grid  system.  Figure  13  shows 
two  typical  examples  of  acceptable  grid  systems.  Here,  physical  (x-y-z)  coordi¬ 
nate  system,  the  block  numbers  and  transformed  coordinates  (^-ti-Q  are  illus¬ 
trated  as  a  reference.  Although  there  is  no  particular  restriction  on  the  primary 
flow  direction  in  the  ALLSPD-3D  code,  the  ^  direction  is  preferred. 
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FIGURE  13.  Coordinates 


FLOW  OVER  A  WEDGE  FLAME  HOLDER 


FIGURE  14-  Two  Dimensional  Planar  and  Axisymmetric  Configurations 


6.4  ALLSPD-3D  For  2-D  Problems 

The  current  ALLSPD-3D  code  can  be  used  to  compute  two-dimensional  or 
axisymmetric  flows.  The  optimal  method  for  calculating  2-D  or  axisymmetric 


flows  is  to  use  the  periodic  boundary  conditions,  which  require  a  minimum  of  2 
grid  points  in  the  periodic  direction.  In  the  current  version  of  the  code,  it  is  sug¬ 
gested  that  the  periodic  boundary  condition  be  used  in  they  or  T|  direction,  if  pos¬ 
sible.  (See  the  I/O  section  for  details  on  defining  periodic  boundaries.)  Users 
must  generate  either  two  parallel  grid  points  for  a  2-D  planar  case  or  2  grid 
points  with  an  angle  for  an  axisymmetric  case  (see  Figure  14 ). 

6.5  ALLSPD-3D  Code  Compilation 

The  code  is  designed  to  run  on  different  computer  platforms  in  either  serial 
or  parallel  modes.  The  compilation  for  different  machines  is  controlled  by  the 
Makefile  using  the  proper  make  commands,  which  are  described  in  this  section. 
But  before  doing  that,  users  must  define  some  parameters  required  to  define  the 
problem  size.  The  parameters  in  the  allspd.h,  which  resides  in  the  subdirectory 
CFD/src/Include,  are  used  to  define  the  problem  size.  These  parameters  are  de¬ 
scribed  as  follows. 

I_MAX:  The  maximum  number  of  grid  points  in  the  i  direction  among  all 
blocks. 

J_MAX:  The  maximum  number  of  grid  points  in  they  direction  among  all 
blocks. 

K_MAX:  The  maximum  number  of  grid  points  in  the  direction  among  all 
blocks. 

UK  TOTAL:  X  '^\imax  +2)x.(jmax  +2)  x(kmax  +2)  where  imax  Jmax  , 

/i=l  fi  Ti  ri  n  n 

and  kmax^  are  the  grid  points  of  the  lath  block  in  i,  j,  k  respectively. 

NUM_ZONES:  Number  of  blocks  defined  in  the  domain. 

NUM_SPECIES:  Number  of  species. 

NUM_REACT:  Number  of  reaction  equations. 

MAX_OBSTACLES:  Maximum  number  of  internal  blockages. 

MAX_DILUTION:  Maximum  number  of  dilution  holes. 
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MAX_TRANSPIR:  Maximum  number  of  transpiration  holes. 

MAX_SPARKS:  Maximum  number  of  ignition  sparks. 

MAX_SPRAY_PARTICLES:  Maximum  number  of  spray  particles. 

MAX_SPRAY_INTEGRATION:  Maximum  number  of  integration  time 
steps  for  the  spray  particles. 

The  make  conunand  is  utilized  for  easy  compilation  on  a  variety  of  computer 
platforms.  But  before  compiling  for  a  new  test  case  on  any  machine,  run  make 
veryclean  to  get  rid  of  any  old  files  leftover  from  previous  compilations.  The 
make  command  targets  for  both  serial  and  parallel  versions  follow. 


TABLE  9.  The  make  Commands  for  Supported  Platforms 


Platform 

Serial  Code 

Parallel  Code 

Cray 

make  cray 

N/A 

HP-UX 

make  hppa 

make  pjtppa 

IBM  RS/6000 

make  rs6k 

make  p_rs6k 

SGl/miXS 

make  sgi5 

make  p^sgiS 

SUN/SunOS4.1 

make  sun4 

make  p_sun4 

SUN/Solaris2 

make  sun4sol2  1 

make  p_sun4sol2 

6.6  Program  Execution  (sans  GUI) 

After  compiling  the  code,  you  can  use  the  make  coimnand  to  install  the  ex¬ 
ecutables  allspdSd  and  pregrid  wherever  you  intend  to  run  your  test  case.  This 
is  done  quite  simply  by  utilizing  the  DEST  definition  in  the  Makefile  with  the 
install  target.  For  example,  if  you  would  like  to  install  the  ALLSPD-3D  execut¬ 
ables  in  the  J.JTests/Swirl-combustion  directory  (relative  to  where  you  built  the 
code),  you  would  run  make  DEST=../../rests/Swirl-combustion  install  after  a 
successful  compilation. 

Once  the  executables  are  installed,  you  need  to  create  the  various  input  files 
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needed  for  your  test  case.  At  a  minimum,  the  allspd.dat,  griddat,  gridlintdat, 
and  prop.dat  files  (or  their  equivalents  named  according  to  your  preferences) 
must  exist  in  the  test  directory.  Other  input  files  may  be  needed  based  upon  the 
details  of  your  particular  test  case.  (See  the  File  I/O  section  for  details  on  file 
descriptions,  formats,  etc.)  After  you  have  all  the  needed  files  in  place,  you  can 
run  the  code. 

If  your  test  case  utilizes  the  turbulence  model,  you  must  first  run  pregrid 
which  will  pre-process  the  grid  and  generate  a  file  called  ynbl.dat.  The  ynbldat 
file  contains  normal  distances  (yn)  to  the  nearest  wall  for  each  point  in  the  grid. 
It  also  contains  a  crude  estimate  of  the  largest  dimension  of  the  grid  (xlmax)  for 
initialization  of  the  turbulence  quantities.  The  options  for  pregrid  are  listed  be¬ 
low  with  the  default  file  names  which  will  be  used  in  their  absence. 


TABLE  10.  The  pregrid  Command  Line  Options 


option 

default  file  name 

example 

-0 

out-gruLdat 

-Oout-griddat 

-f 

gridlink.dat 

•fgridlinLdat 

-g 

griddat 

“ggriddat 

-y 

ynbldat 

~y ynbldat 

During  the  execution  of  pregrid,  various  bits  of  information  will  appear  on 
the  standard  output  to  inform  you  of  the  pre-processor’s  progress  as  it  reads  and 
processes  grid.dat.  When  pregrid  is  finished,  the  out-grid.dat  file  will  contain  a 
synopsis  of  the  data  found  in  gridlink.dat  for  your  verification.  This  pre-process¬ 
ing  step  is  needed  only  once  unless  you  change  the  grid.  If  you  decide  to  restart 
a  test  case,  there  is  no  need  to  re-run  pregrid. 

Now  you  are  ready  to  run  the  flow  solver.  In  the  simplest  case,  all  you  have 
to  do  is  run  ollspdSdJhQ  various  options  for  using  alternate  file  names  are  listed 
in  the  following  table.  (Again,  see  the  File  I/O  section  for  details  on  file  formats 
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and  descriptions.) 

TABLE  11.  The  alkpdSd  Command  Line  Options 


option 

default  file  name 

example 

-i 

allspddat 

-lallspddat 

-0 

outdat 

•‘Ooutdat 

-r 

res.dat 

'■rres.dat 

-k 

reske.dat 

•kreske.dat 

-P 

allspd 

’‘pailspd 

-s 

save.dat 

•ssave.dat 

-F 

face.dat 

•¥face.dat 

-g 

griddat 

-ggriddat 

-f 

gridlink.dat 

-fgridlinLdat 

-u 

uvwtdat 

•nuvwtdat 

-I 

spray.dat 

Aspray.dat 

-c 

prop.dat 

~cprop.dat 

-y 

ynbl.dat 

-yynbldat 

As  allspdSd  starts  up,  it  reads  the  various  input  files  and  reports  a  synopsis 
to  outdat  and  standard  output.  As  the  calculation  progresses,  numerous  messag¬ 
es  will  appear  on  the  standard  output  and  in  various  output  files.  When  the  code 
finishes,  your  directory  will  contain  output  files  and  (if  selected)  restart  files. 

To  continue  the  calculation,  you  need  to  edit  allspddat  accordingly  and  ex¬ 
ecute  allspdSd  again.  The  code  will  then  read  in  the  files  save.dat  andface.dat, 
in  addition  to  all  the  other  inputs.  This  will  include  the  ynbLdat  file  previously 
generated  hy  pregrid  if  a  turbulent  flow  is  being  simulated.  You  should  see  the 
code  report  a  successful  restart  message  after  reading  all  the  necessary  files.  Af¬ 
ter  this  message,  the  code  executes  as  normal. 


6.7  Parallel  Implementation 

6.7.1  Theory 

The  theory  behind  the  parallel  processing  version  of  AIJLSPD-3D  is  quite 
simple:  the  code  is  inherently  divided  in  the  data  domain,  therefore  domain  de¬ 
composition  is  used.  (How’s  that  for  simplicity?)  The  multiple  grid  zone  feature 
provides  natural  dividing  lines  in  the  data  for  decomposing  the  problem  onto 
multiple  processes,  i.e.,  each  grid  zone  is  a  natural  candidate  for  parallel  pro¬ 
cessing.  Choosing  which  data  to  send  between  processes  is  also  quite  simple:  es¬ 
sentially  whatever  variables  are  copied  out  of  the  1-D  arrays  and  into  the  3-D 
arrays  during  CALL  CONV13  in  the  serial  code.  This  data  is  transferred  between 
processes  using  a  message-peissing  library  instead  of  shared  memory  to  maxi¬ 
mize  portability. 

The  transferred  data  supplies  exactly  the  information  a  process  needs  from 
its  neighboring  zones  without  requiring  memory  space  for  all  of  the  data  in  other 
zones.  Because  each  process  needs  data  for  only  its  own  grid  zone  (including 
those  ghost  cells  which  actually  belong  to  neighboring  zones),  each  process  only 
needs  enough  memory  for  the  largest  zone.  (ALLSPD-3D  sizes  all  of  its  3-D 
data  arrays  based  on  the  largest  zone  in  the  problem.)  This  reduced  memory  fea¬ 
ture  of  parallel  processing  can  be  very  beneficial  with  large  problem  sizes.  Also, 
since  each  process  only  calculates  data  on  its  zone,  the  time  needed  to  calculate 
a  single  iteration  is  reduced  to  (more  or  less)  the  time  needed  for  the  most  nu¬ 
merically  intensive  zone.  The  only  cost  for  these  great  benefits  of  parallel  pro¬ 
cessing  is  the  time  it  takes  to  transfer  data  between  neighbors. 

6.7.2  Implementation 

The  PVM  (Parallel  Virtual  Machine)  message-passing  library  was  chosen 
because  of  its  wide  acceptance,  installed  user  base,  and  portability.  PVM  is  used 
in  a  wide  variety  of  applications  on  numerous  architectures  and  has  become  a 
de-facto  standard  for  message-passing  libraries.  The  PVM  library  has  many  fea- 
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tures  including  spawning  of  processes  on  a  virtual  machine  and  the  communi¬ 
cation  of  various  message  types  between  architectures  which  may  have 
inherently  different  data  structures.  These  features  are  used  in  the  parallel  ver¬ 
sion  of  ALLSPD-3D.  (For  more  information  on  PVM,  visit  http:// 
WWW.  epm.  omLgov/pvm/. ) 

When  the  parallel  version  of  ALLSPD-3D  is  executed,  it  automatically 
spawns  a  child  process  for  each  grid  zone  except  the  one  zone  that  the  parent 
runs.  The  parent  process  then  becomes  the  controlling  node  in  the  parallel  cal¬ 
culation  and  collects  standard  output  from  all  the  children  for  display  in  the  con¬ 
troller’s  output.  (This  can  produce  some  ugly  output,  but  it  is  still  the  best  way 
to  gatbp-r  all  the  information.)  Each  process  then  runs  as  a  single  zone  calcula¬ 
tion  with  needed  data  passed  between  neighbors  across  the  zone  interfaces. 

The  data  that  is  transferred  between  grid  zones  (and  therefore  processes)  in¬ 
cludes  both  geometric  data  and  flow  data.  The  geometric  data  is  passed  between 
neighbors  twice  at  startup:  first  to  transfer  X,  Y,  Z  coordinate  information  about 
neighboring  grids,  then  to  update  the  various  derived  quantities.  This  is  neces¬ 
sary  since  each  process  reads  in  the  entire  grid,  but  only  stores  the  data  for  its 
own  zone.  Then  the  flow  information  (pressure,  velocities,  enthalpy,  species 
concentrations,  viscosity,  etc.)  is  transferred  between  neighbors  after  initializa¬ 
tion  and  after  each  iteration.  To  minimize  the  size  of  each  message,  the  data  is 
transferred  in  individual  messages,  e.g.,  one  message  for  X  values,  another  for 
U  velocity  values,  etc.  (All  these  messages  are  sent  using  PVM’s  default  data 
format  which  allows  a  heterogeneous  mix  of  machines  to  transparently  share 
data  without  concern  for  data  representation  details  on  each  architecture.)  All 
this  message  passing  fills  in  the  data  at  the  ghost  cells  for  each  zonal  interface. 

Sending  data  to  neighbors  after  each  iteration  impacts  the  behavior  of  the 
code  both  superficially  and  significantly.  The  transfer  of  data  obviates  any  need 
to  move  data  in  and  out  of  the  3-D  arrays.  It  also  forces  the  ghost  cells  to  contain 
data  which  is  lagged  a  full  iteration,  i.e.,  the  data  one  zone  uses  from  a  neighbor 
is  old  information.  In  the  serial  version  of  ALLSPD-3D,  the  zone  interfaces  can 
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use  data  updated  during  the  current  iteration  depending  upon  the  zone  ordering. 
(This  fully-lagged  boundary  condition  at  the  zone  interfaces  is  not  necessarily 
bad.  Some  test  cases  have  actually  had  unproved  convergence  in  parallel  pro¬ 
cessing,  presumably  from  the  dampening  of  errors  across  lagged  boundaries.) 
This  is  a  fundamental  change  in  the  propagation  of  information  through  the  flow 
field  and  results  in  differing  convergence  histories  for  the  parallel  and  serial  ver¬ 
sions  of  ALLSPD-3D. 

By  sending  and  receiving  data  between  neighbors  and  then  synchronizing 
with  only  neighboring  zones,  a  loose  synchronization  of  the  total  calculation  is 
achieved.  Basically,  the  process  that  takes  the  longest  to  calculate  its  zone  wiU 
force  all  the  others  to  wait  for  it  because  of  the  interdependence  of  the  grid 
zones.  This  also  means  that  load  balancing  will  be  best  achieved  by  dividing  the 
grid  into  more  or  less  equally  sized  zones.  (The  only  caveat  here  is  that  the  spray 
model  requires  additional  resources  and  therefore  fewer  grid  points  should  be 
used  in  the  zone  containing  spray  calculation  to  achieve  load  balance.) 

6.73  Notes  About  Coding  Details 

The  ALLSPD-3D  code  is  primarily  written  in  Fortran,  but  much  of  the  par¬ 
allel  coding  is  in  C.  This  results  in  the  sharing  of  some  data  between  C  and  For¬ 
tran  which  is  not  pretty.  Fortunately,  the  basic  data  transfer  of  geometric  and 
flow  data  is  handled  by  Fortran  code  (to  avoid  the  perverse  practice  of  accessing 
Fortran  COMMON  BLOCK  data  from  C)  and  the  parallel  management  is  per¬ 
formed  by  C  code. 

If  you  decide  to  modify  the  basic  sequencing  of  sends,  receives,  and  syn¬ 
chronizations,  be  very  careful.  The  checking  for  abort  signals  is  critical  to  avoid 
deadlocking  the  calculation  because  one  node  has  decided  to  exit  while  another 
is  waiting  for  some  message  that  will  never  arrive. 

The  amount  of  data  transferred  between  neighbors  is  based  on  the  largest 
possible  size  of  a  face  because  of  limitations  in  Fortran.  This  means  that  mes¬ 
sages  could  be  greatly  oversized  for  particular  problems,  especially  2-D  calcu- 
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lations.  If  you  wish  to  reduce  the  data  passed  between  zones  to  a  more 
reasonable  size  for  a  specific  test  case,  you  can  change  the  definition  of  MAX- 
FACE  in  the  CFD/src/Include/allspd.h  file.  (But  be  very  careful.)  Also,  the 
communication  time  for  data  transfer  can  be  improved  for  particular  test  cases. 
For  generality,  the  parallel  version  of  ALLSPD-3D  uses  the  default  data  transfer 
mode  of  PVM.  If  your  virtual  machine  is  homogeneous,  you  could  switch  to  raw 

mode  or  data  in  place  mode.  Just  look  for  PVM_DATA_XFER_C  and 
PVM_DATA_XFER_F  in  the  CFD/src/Include/allspd.h  file  to  change  the  data 
transfer  mode. 

Currently,  ALLSPD-3D  has  no  provision  for  spawning  children  to  specific 
machines  based  upon  CPU  speed,  grid  zone  size,  etc.  Therefore,  the  gridding 
will  have  a  direct  influence  on  the  load  balance  of  the  overall  calculation.  Also, 
the  PVM  daemons  will  decide  which  host  will  calculate  which  zone,  usually 
without  regard  to  existing  CPU  loads.  This  means  that  you  can  easily  wind  up 
with  processes  spawned  in  the  worst  (or  best)  possible  allocation,  especially  in 
highly  heterogeneous  virtual  machines.  Another  impact  from  PVM  is  the  con¬ 
nection  used  between  machines.  A  non-shared  communication  medium  (e.g., 
ATM,  high-speed  switches,  etc.)  is  best.  Standard  ethemet  is  very  common,  but 
it  is  also  the  slowest,  especially  when  all  the  processes  need  to  communicate  at 
the  same  time. 

6.7.4  Compilation  and  Run-time  Details 

To  compile  and  run  the  parallel  processing  version  of  ALLSPD-3D,  you 
must  have  PVM  3.3  installed  on  all  hosts  that  you  intend  to  use.  The  Makefile 
uses  the  environment  variable  PVM_ROOT  and  PVM_ARCH  to  know  where 
PVM  include  files  and  libraries  are  located,  so  make  sure  PVM  is  set  up  prop¬ 
erly.  Then  compilation  is  just  a  matter  of  issuing  the  proper  command.  (I  recom¬ 
mend  that  you  always  ran  make  veryclean  whenever  starting  a  new  test  case  or 
building  for  a  different  architecture.)  Assuming  you  want  to  ran  ALLSPD-3D 
in  parallel  on  a  cluster  of  IBM  RS/6000  workstations,  you  would  run  make 
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p_rs6k.  The  make  command  will  then  compile  the  various  modules  of  ALLSPD- 
3D  including  those  needed  for  parallel  processing.  The  pre-processing  of  the 
source  code  wiU  automatically  include  those  sections  of  code  needed  for  parallel 
processing  and  leave  out  sections  not  needed.  If  your  cluster  of  workstations  is 
heterogeneous,  just  build  with  the  correct  parallel  target  (e.g.,  make  pjsgi  for 
SGI  machines)  and  install  the  resulting  executables  in  separate  directories  ac¬ 
cording  to  the  needs  of  PVM.  (And  don’t  forget  to  run  make  veryclean  between 
builds.) 

After  the  executable(s)  are  built  and  installed  in  the  proper  locations,  run 
pregrid  as  normal,  then  start  PVM  with  a  hostfile  which  is  configured  for  your 
particular  cluster  and  your  problem.  (This  is  basically  a  PVM  issue  and  will  not 
be  covered  here  except  to  say  that  the  ep  [executable  path]  and  wd  [working  di¬ 
rectory]  options  in  the  hostfile  are  critical.  Note  that  the  GUI  will  handle  these 
details  for  you.)  After  the  pvm  daemons  are  running  on  all  the  hosts,  you  start 
up  ALLSPD-3D  pretty  much  as  normal.  Be  careful  not  to  use  pathnames  on  het¬ 
erogeneous  virtual  machines  since  ALLSPD-3D  will  spawn  processes  on  the 
virtual  machine  using  the  same  command  line  you  issue.  After  this,  each  process 
will  read  all  the  input  files  but  only  use  the  parts  it  needs,  e.g.,  each  process  will 
read  the  entire  grid  but  store  only  its  zone.  Whenever  a  process  writes  out  an  out¬ 
put  file,  it  uses  the  name  supplied  on  the  command  line  (or  the  default  name) 
with  a  grid  zone  extension,  i.e.,  the  process  calculating  zone  1  will  output  its  tur¬ 
bulence  residual  data  for  its  zone  (and  only  its  zone)  to  reske.dat-001  instead  of 
reske.dat.  A  slight  exception  to  this  rule  is  the  global  residual  file  res.dat.  Each 
process  gets  its  own  residual  file  for  the  flow  equations  on  its  zone.  The  control¬ 
ling  node,  however,  also  outputs  a  residual  for  the  entire  flow  field  to  res.dat  af¬ 
ter  gathering  information  from  all  the  processes  and  calculating  a  residual 
spanning  the  entire  grid.  After  all  the  processes  finish  the  calculation,  you 
should  probably  halt  the  PVM  daemons. 

When  the  code  finishes,  you  will  have  numerous  files  because  of  each  pro¬ 
cess  dumping  out  its  own  results.  The  solution  files  are  single  zone  solutions  in 
multi-zone  Plot3D  format,  but  can  be  easily  assembled  into  a  single  file.  The  re- 
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start  ffles  {save.dat-nnn  and  face.dat-nnn)  are  inherently  incompatible  with  the 
serial  code.  These  files  contain  only  the  data  needed  for  a  single  zone  to  restart 
the  calculation  while  the  serial  code  needs  all  the  zones  in  a  single  file.  Assem¬ 
bling  the  numerous  restart  files  of  a  parallel  calculation  into  a  single  file  for  use 
by  the  serial  code  is  probably  possible  but  likely  to  be  a  painful  experience  and 
is  not  recommended. 
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VII :  Examples 


This  section  describes  two  test  cases  computed  using  the  ALLSPD-3D  code. 
These  cases  are  designed  to  help  the  users  obtain  the  experience  needed  in  set¬ 
ting  up  the  ALLSPD-3D  input  data  files.  The  ALLSPD-3D  solutions  are  also  in¬ 
cluded  so  that  user-obtained  solutions  can  be  compared  with  the  other 
computational  results  and  also  with  the  experimental  data.  Re-running  these 
cases  will  help  users  gain  experience  needed  to  create  the  AULSPD-3D  inputs 
and  is  a  highly  recommended  step  for  new  users.  The  two  test  cases  included  in 
this  section  are: 

(1)  3-D  non-reacting  laminar  and  turbulent  transition  duct  flows 

(2)  Axisynunetric  turbulent  swirl-can  combustor  flow  with  spray 


7.1  3-D  Non-Reacting  Transition  Duct  Flow 

The  fully  turbulent  and  three-dimensional  circular  to  rectangular  transition 
duct  problem  of  Davis"^^  was  chosen  as  the  first  test  case  because  of  the  exten¬ 
sive  and  detailed  experimental  data  available.  The  grid  has  52,521  cells 
[41x21x61]  and  three  typical  streamwise  stations  are  shown  in  Figure  15.  The 


tests  conditions  are  shown  in  Table  12. 


Rej) 

195,000 

"^ref 

298K 

Uref 

29  m/s 

TABLE  12.  Reference  Conditions  for  the  Transition  Duct 


FIGURE  IS.  The  Tlransition  Duct  Grid  at  Three  Stream  wise  Stations 

Before  running  the  fully  turbulent  case,  we  investigated  a  laminar  flow  in  the 
transition  duct  by  reducing  the  inlet  flow  speed.  The  Re5molds  number  is  about 
195  for  the  laminar  flow.  The  input  data  files,  gridlintdat,  allspd.dat,  and 
prop.dat,  needed  to  set  up  this  laminar  flow  are  listed  in  Figure  16  to  Figure  18. 
The  velocity  vectors  obtained  from  the  ALLSPD-3D  solution  are  shown  in 
Figure  19,  and  the  convergence  history  is  shown  in  Figure  20. 
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FIGURE  16.  The  gridUnkJat  Input  FOe  for  the  Tl-ansition  Duct 
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FIGURE  17.  The  aUspd.dat  Input  File  for  the  Laminar  It-ansition  Duct 
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FIGURE  18.  The  prop.dat  Input  file  for  the  Ik-ansition  Duct 


FIGURE  19.  Velocity  Vectors  for  the  Laminar  lYansition  Duct 


Transition  Duct  -  3D 


FIGURE  20.  Convergence  History  for  the  Laminar  and  Itarbulent 
lYansition  Duct 

The  fully  turbulent  transition  duct  flow  case  requires  the  input  data  file  all- 
spd.dat  listed  in  Figure  17.  The  gridlintdat  arid  prop.dat  are  identical  to  those 
for  the  previous  laminar  case.  For  the  turbulent  flow,  it  is  noted  that  iuvwt  is 
equal  to  1  and  therefore  the  uvwtdat,  an  initial  condition  file,  has  to  be  provided 
as  a  separate  input  file.  Figure  22  shows  the  velocity  vectors  on  the  two  symme¬ 
try  planes  of  the  transition  duct.  Comparisons  of  the  u  velocities  at  various 
streamwise  stations  along  the  two  major  axis  with  experimental  data  are 
shown  in  Figure  23.  This  figure  shows  that  reasonable  agreement  with  expert- 


mental  data  can  be  obtained  with  the  ALLSPD-3D  code.  The  streamwise  locations 
of  the  velocity  profile  comparisons  shown  in  Figure  23  are  illustrated  in  Figure  22. 
Here,  additional  velocity  vectors  on  the  synunetry  planes  are  shown.  Figure  20  also 
shows  the  convergence  history  for  this  turbulent  flow  ALLSPD-3D  calculation. 
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2000 
100 
100 
10  1 
0 
1 

S.OOOOOOe-Ol  l.OOOOOOe+00 
l.OOOOOOe+10  l.OOOOOOe+10  l.OOOOOOe+10 
l-050000e+00  1.050000e+00  1.050000e+00 
1 

l.OOOOOOe-01 

le-13 

1 

1 

1 

2.995000e+01  O.OOOOOOe+00  O.OOOOOOe+00 

2.970000e+02 

1.013250e+05 

1.021400e-01 


0 

0 

30000 

50000 

0 

1.600000e+03 

O.OOOOOOe+00 

0 

0 

12345  11  12  13  7900000 
2 

'02_T<1000'  2.335700e-01 

'N2_T<1000'  7.664300e-01 

FIGURE  21.  The  allspd.dat  Input  FUe  for  the  Ttarbulent  Transition  Duct 


FIGURE  22.  Velocity  Vectors  for  the  Thrbulent  Transitions  Duct 
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U  valocity  along  seml-maior  axis 


-  MJLSn>-30.  41x21xS1 


u/“o 


T 


U  velocity  along  semi-minor  axis 


FIGURE  23.  Velocity  Profile  Comparison 


7.2  Axisyinmetric  Swirl-Can  Combustor  Flow  with 
Spray 

The  second  sample  test  case  was  used  to  validate  the  current  turbulence 
model  and  the  spray  model  for  swirling  flow.  The  experiment  was  conducted  by 
Roback  and  Johnson^^  ^  non-reacting  turbulent  flow.  However,  we  use  this 
case  to  further  demonstrate  the  ALLSPD-3D  spray  model  by  initiating  combus¬ 
tion  in  the  flow  field.  Since  this  is  an  axisymmetric  flow  (see  Section  6.4  of  the 
ALLSPD-3D  manual),  only  two  grid  points  in  the  circumferential  direction  and 
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the  periodic  boundary  condition  are  used.  Here,  the  combustion  is  initiated  by 
injecting  the  liquid  methanol  fuel  into  the  flow.  The  input  data  files  needed  are 
listed  in  Figure  24  to  Figure  26.  It  should  be  noted  that  iuvwt  is  equal  to  1  and 
therefore  the  uvwtdat  input  file  is  also  needed.  The  reference  conditions  are  pro¬ 
vided  in  Table  13. 


Rep 

61,180 

"^ref 

300  K 

Uief 

16  m/s 

TABLE  13.  Reference  Conditions  for  the  Swirl  Can  Combustor 


0 

6 

49  2  38 
15  2  38 
19  2  38 
49  2  24 
15  2  24 
19  2  24 
4  1  6  6  4  7 
4  4  6  6  4  7 
2  4  6  6  4  7 
4  1  6  6  3  4 

4  4  6  6  3  4 
2  4  6  6  3  4 

2  0  0  0  4  0 

3  1  0  0  5  0 
0  2  0  0  6  0 

5  0  0  0  0  1 

6  4  0  0  0  2 
0  5  0  0  0  3 
0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0  0  0  0  0  0 
0  0  0  0  0  0 
101325  0000 
0  0  0  0  0  0 
0  0  0  0  0  0 
101325  0000 


0 


0 


FIGURE  24,  The  gridlink.dat  Input  File  for  the  Swirl  Can  Combustor 
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1.050000e+00  1.050000e+00  1.050000e+00 
1 

l.OOOOOOe-01 

le-13 


1 

1 

1 

1.600000e+01  O.OOOOOOe+OO  O.OOOOOOe+OO 
3.000000e+02 
1.013250e+05 
6.100000e-02 
0 
0 

300 

50000 

1 

1  25  40  1  2  2  30 
1.600000e+03 
5.000000e+07 
0 
1 
1 

30  30  1  2  5  24 
1  O.OOOOOOe+OO  O.OOOOOOe+OO 
l.OOOOOOe+01  3.700000e+01 
300 
20 
0 

3.000000e+02 

1.907808e”04 

3.000000e+01 

1  2  3  4  5  15  16  17  18  19  11  12  13  14  7 


5 

'CH40'  O.OOOOOOe+OO 
'02'  2.335700e-01 

'C02'  O.OOOOOOe+OO 
'H20'  O.OOOOOOe+OO 
'N2'  7.664300e-01 


FIGURE  25,  The  aUspd,dat  Input  File  for  the  Swirl  Can  Combustor 


ntl«: 


Speclva 

♦OMO' 

'02' 

'C02' 

'H20' 


Special  itetM 
'CHi0‘ 
.9774$0e«01 
•02' 

.314340e-03 

'C02' 

.3143401*03 

'H20' 

.314340e*03 


.314340e*03 
Specie*  nom 
'CH40' 

•02' 

'C02' 

'H20' 

'N2' 

Species  luuM 
'CH40' 

•02' 

•C02' 

•H20' 

•N2' 

Species  none 
'0140' 


'ICO' 


Species  name 
'CK40' 


'02' 

'C02' 

'H20' 


Molecular  Welghc 


3.204200e*01 

3.199680e*01 

4.4009S0e*01 

1.601528**01 

2.801348**01 

CpO (as) 

7.8S3B14e-01 

3.086805**00 

2.655681e*00 

3.882612**00 

3.430237e*00 

EuO<ns) 

-2.233179*-06 

4.470847e-06 

-3.2744140-07 

-2.400723e-06 

4.305119e-06 

CdOlns) 

-1.7932590-02 

7.33l868e-04 

-1.1656600-02 

2.605747e-04 

S.409501e-03 

SiBB(as) 

3.626000**00 

3.4670000*00 

3.941000**00 

2.641000**00 

3.79e000**00 

EcUcIns) 

4.818000**03 

1.067000**02 

1.952000**02 

8.091000**02 

7,140000**01 


Cpl(ns) 

1.6862240-02 

l,731246*-03 

7.327241e-03 

1.320022*'04 

6. 646834*- OS 

Eul(ns) 

4.243590*-0e 

5.9215S4e-08 

5.787S88O-08 

4. 07582 6*' 08 

S.  05079 1*-08 

Cdllasl 

9,7496110-05 

9.096378e-0S 

9. 871455*- 05 

3.825095e-05 

7.014240e-05 


Cp2(ns) 

-9. 249801 e-06 

-8.450126O-07 

-4.5854030-06 

1.5128170-06 

6.908434e-07 

Bu2 (ns) 

-1. 0913660-11 

-2. 007 533 0-11 

-2.163243e-ll 

2.725794e-12 

-1.8178420-11 

Cd3(nsi 

6.49S632e-08 

-2.590399e-08 

-2.4718800-08 

7.966066e-08 

-6.616883e-09 


Cp3<as) 

2. 629736 e-09 

5. 6579770-10 

1. 3731210-09 

-6.7634770-10 

-1.741 6860-10 

Eu3(nsl 

2.1672380-15 

S.556300O-1S 

6. 0118530-15 

-3. 4474040-15 

5.223197e-15 

Cd3(as) 

-3  .108S40e-ll 
6.7982390-12 
4.44S081e-12 
-2.3467700-11 
-6.710531*- 13 


Descrlpclon  of  Reaccloa  0:  vacanc 

Roac  •  Porx_Pr^.S3(p.  Porv_AcCEaor.  Pory_T«a©_coef . 

1  4.600000e*ll  -1.509810**04  0.000000**00 

Reac  *  B*c>;_Pr«_B35>.  Back;>cCEaer.  B*c)e..T«BS)_coef . 

1  0.000000**00  0.000000**00  0.000000**00 
Reac  •  rorvard  Stolcloowtric  Coefficieacs  (Species  — >1 
1  1.000000**00  1.500000**00  0.000000**00  0.000000**00  0.000000**00 
Reac  #  BaOcward  Stoiclo«crlc  Coefficients  {Species  -->) 

1  0.000000**00  0.000000**00  1,000000**00  2.000000e*00  0.000000e*00 
Roae  »  Third  Body  Zfficloacl**  (species  — » 

1  0.0  1.000000e*00  1.000000**00  0.000000**00  0.000000**00  0.000000**00 
1  0.0  1.000000**00  1.000000**00  0.000000«*00  0.000000**00  0.000000«*00 


'0140' 

7.950000**06  5.130000**02  3.377100e*02  1.750000**02 
1.176600**06  2.490000**03  2.000000e-01  7.915000**02 
-7.289360**00  1.536790**00  -3.083670**00  -1.024560**00 


Cp4(nB) 

-3.033B81e-13 

-7.0606560-14 

-1.655748«-13 

9.806783*-14 

3.2058380-14 

£u4(ns) 

-1.7983100-19 

-6.3511630-19 

-6.7507280-19 

S.6S5870*-19 

-6. 0778580-19 

Cd4(ns> 

4.266252e-lS 

-7.1535890-16 

-3.0S42S8*-16 

2.572095*-15 

2,4651250-16 


-2.506e43**04 

-3.611010**00 

-4.733380**04 

-2.908730**04 

-3.474180**00 


FIGURE  26.  The prop.dat  Input  File  for  the  Swirl  Can  Combustor 


The  grid  used  by  the  ALLSPD-3D  code  on  the  swirl  can  combustor  is  shown 
in  Hgure  28.  In  this  figure,  every  other  point  is  plotted  with  each  zone  in  sepa¬ 
rate  colors  for  clarity.  The  grid  is  comprised  of  six  blocks  as  detailed  in 
Table  14.  The  velocity  vectors  obtained  from  the  ALLSPD-3D  solution  at  2000 
iterations  is  shown  in  Figure  29.  This  solution  was  obtained  using  the  ALLSPD- 
3D  code  in  the  parallel  mode  where  the  calculation  was  performed  on  6  work¬ 
station  nodes.  The  velocity  vectors  are  colored  according  to  flow  field  tempera¬ 
ture  to  show  that  the  fuel  has  ignited.  This  ignition  produces  a  flowfield  where 
the  temperature  ranges  from  300  K  to  2100  K.  In  Figure  29,  the  color  blue  is 
used  to  indicate  low  temperature  and  the  color  red  is  used  to  indicate  high  tem¬ 
perature.  The  convergence  history  for  this  calculation  is  shown  in  Figure  27. 
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Block 

Block  Dimensions 

Number  of  Points 

1 

49  X  2  X  38 

3724 

2 

15x2x38 

1140 

3 

19  X  2  X  38 

1444 

4 

49  X  2  X  24 

2352 

5 

15  X  2  X  24 

720 

6 

19  X  2  X  24 

912 

Total 

10292 

TABLE  14.  Swirl  Can  Combustor  Grid  Block  Sizes 


Swirl  Can  -  2D 


FIGURE  27.  Convergence  History  for  the  Swirl  Can  Combustor 
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FIGURE  28.  Grid  for  the  Swirl  Can  Combustor 


FIGURE  29.  Temperature  Colored  Velocity  Vectors 
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Appendix  A 
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Jacobian  Matrices 

The  number  of  species  equations  is  assumed  to  be  five  Qf=5)  in  the  Jacobian 
matrices  given  below 
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p{uU+a^H)  pCvlZ  +  ajH)  p{wU  +  a^H) 

pl/<t> 

.(2)  _  pi7  0 

[o  pU 


(A  2) 


> 


I  (3) 


pt/©,  p[/r,Q2 

pf/02  pf/yj^ 

pI/FjQj  pC/yjQj  ^UYj£l^ 

pl/y^Qj  pUY^^^  p  1/74^3  pl/04 


(A3) 


j/  =  aiM  +  a2v  +  a3W 


(A  4) 
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“l  =  «2  =  “3  = 

(A  5) 

'  Iw,-  wj  C^T 

(A  6) 

> 

0,.  =  i  +  yA  = 

(A  7) 

here  W  is  the  molecular  weight  of  the  gas  mixture  and  is  the  mixture  spe¬ 
cific  heat  Cj,  .  The  Jacobian  matrix  is  obtained  by  letting 
0^=11^  and  Jacobian  matrix  C^'"^  is  obtained  by  letting 

The  viscous  term  Jacobian  has  the  following  form 


(A  8) 


(A  9) 
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In  order  to  reduce  the  stiffness  of  the  chemical  species  source  term,  the 
source  terms  are  treated  implicitly.  The  chemical  species  source,  Equation(24), 
is  a  function  of  temperature  and  species  concentration  and  therefore  a  chain  rule 
is  required  to  evaluate  the  source  term  in  general.  The  linearization  of  the  left 
hand  side  leads  to  the  following  expressions: 
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+  MW,.X  CR,.,  [iT/ 


Y.  Jll 
^  /  =  1 


C,.  lar 


= MW, 2  n<‘-?  s^"-  n<“ 


/=1  ^/=1 


/=1  ^/=1 


‘  =  MW,X  CR,J  Ja,(k) -MUcr ‘^-^Uik) 


(A  14) 


where,  Ej^ ^  and  Kj;  ^  are  function  of  the  index. 


The  turbulence  source  term  shown  in  Equation(35)  can  be  very  stiff  nunier- 
ically.  Therefore,  the  ALLSPD-3D  code  solves  the  equation  in  a  point  implicit 
manner.  The  linearization  of  the  turbulent  sou^  term  required  the  use  of  the 
following  source  Jacobian  term.  This  formulation  has  been  optimized  for  con¬ 
vergence.  This  form  has  been  tested  for  a  number  of  flow  situations  and  model 
formulations 


pit  k 


£  ^£2 


(A  15) 
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Appendix  B 


Thennodynaniic  Database  Information 

The  database  used  by  the  GUI  is  an  assembly  of  gaseous  and  liquid  proper¬ 
ties  for  several  fuels  and  oxidizers.  Most  of  the  thermodynamic  and  transport 
properties  were  fourth  order  polynomial  curve  fits  to  those  provided  by 
McBride  et.al.^^  The  original  data  from  this  report  are  split  into  two  temperature 
ranges.  These  curves  have  been  combined  into  a  single  fit  for  a  temperature 
range  300K  to  SOOOK.  The  diffusion  constants  used  in  the  Chapman-Enskog 
formulation  for  species  diffusion  were  taken  from  Svehla.^^  The  liquid  proper¬ 
ties  came  from  Reid.^^  An  example  of  a  database  entry  is  shown  in  Table  15 
and  described  in  Table  16  (Superscripts  are  used  to  relate  Table  15  to  the  appro¬ 
priate  descriptions  in  Table  16 ). 


TABUE  15.  Sample  GUI  Database  Entry 


CH40® 

Methanol^ 

32.0420*^ 

D 

■■ 

7.8538i449E-01® 

1.68622442E-02 

-924980079E-06 

2.62973642E-09 

-3.03388054E-13 

-2.23317873E-06” 

4.24359048E-08 

-1.09136602E-11 

2.16723754E-15 

-1.79831037E-19 

-1.79325923E-02^ 

9.74961128E-05 

6.49563168E-08 

-3.10853977E-11 

4.26625171E-15 

3.626^ 

481.8^ 

-2.50864348B+04^ 

1,97745991E+01^ 

79.5c5^ 

513.0° 

337.71^ 

175.0^ 

1178.6c3^ 

2.49e3^ 

0.2'^ 

791.5° 

-728936^^ 

1.53679 

-3.08367 

-1.02456 

r 


TABLE  16.  Description  of  GUI  Database  Entries 


Entry  Name 

Entry  Description 

^Entry  Number 

Database  book  keeping 

^Chemical  Name 

The  chemical  name  of  the  species  in  upper  case  letters.The  GUI  is  case 
sensitive.  However,  a  specie  such  as  iso-butane  would  be  iC3H8 

^Common  Name 

A  more  general  specie  name 

^Molecular  Weight 

The  species  molecular  weight 

^Gas  Phase  Flag 

Indicates  to  the  GUI  that  a  gas  phase  exists  in  the  database 

liquid  Phase  Flag 

Indicates  to  the  GUI  that  a  liquid  phase  exists  in  the  database 

^Specific  Heat  Coefficients 

The  coefficients  used  in  the  fourth  order  polynomial  fit  for  specific  heat 

^Viscosity  Coefficients 

The  coefficients  used  in  the  fourth  order  polynomial  fit  for  viscosity 

^Thermal  Conductivity 

Coefficients 

The  coefficients  used  in  the  fourth  order  polynomial  fit  for  thermal 
conductivity 

Collision  cross  section  for  diffusion 

Kt 

Temperature  for  collision  integral  correction 

^Enthalpy  Constant 

Integration  constant  for  enthalpy 

^Entropy  Constant 

Integration  constant  for  entropy 

^Critical  Pressure 

Critical  pressure  for  species 

^Critical  Temperature 

Critical  temperature  for  species 

^Boiling  Temperature 

Boiling  tenq>erature  for  species 

^Melting  Temperature 

Melting  temperature  for  species 

^Latent  Heat 

Latent  heat  of  vaporization  at  the  spray  ten^rature  (293K) 

^Liquid  Specific  Heat 

Specific  heat  of  the  liquid  at  the  spray  temperature  (293K) 

“^Liquid  Thermal  Conductivity 

Thermal  conductivity  of  the  liquid  at  the  spray  ten^rature  (293K) 

^Liquid  Density 

Density  of  the  liquid  at  the  spray  temperature  (293K) 

^V^r  Pressure  Coefficients 

Coefficients  for  the  polynomial  fit  to  fuel  vs^^or  pressure. 

i 
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15 


rholpar  density  (kg/mO 

vpl,vp2|Vp3,vp4  Critical  pressure  curve  fit  coefficients 
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Dynamic  viscosity;  kg/(m  s)  ^-Dir  =  Top  (kmax! 

Kinematic  viscosity:  m^/s 
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